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are  presented  for  three  truss  structures  (10  member  plane  cantilever  truss,  200  member  plane 
truss,  63  member  space  truss),  three  idealized  wing  structures  (18  element  wing  box  beam, 

39  element  rectangular  wing,  150  element  swept  wing),  and  three  framed  structures  (one- 
bay  two-story  plane  frame,  two-bay  six-story  plane  frame,  48  element  space  frame). 

Results  obtained  with  the  substructuring  formulation  are  compared  first  with  results 
obtained  without  substructuring  and  then  with  results  obtained  with  other  methods. 
is  concluded  that  the  algorithm  with  substructuring  is  more  efficient  than  the  one*  - 
without  substructuring,  with  computational  times  reduced  by  up  to  66%.  Finally,  results 
obtained  with  the  present  method  are  compared  with  other  available  results  and  a reduction 
in  computing  time  per  iteration  by  factors  ranging  from  2.4  to  19  is  obtained.  Based 
on  these  comparisons  and  the  range  of  problems  solved,  it  is  shown  that  the  method  is 
"general"  and  computationally  efficient.  Finally,  extensions  of  the  method  to  fail- 
safe design  of  large-practical  structures  are  discussed. 
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This  report  represents  the  Ph.D.  research  of  Dr.  A.h.  Govil,  done 
under  the  supervision  of  Professors  J.S.  Arora  and  E.J.  Haug.  The 
research  was  initiated  and  chapters  one  through  four  completed  under 
U.S.  Army  Research  Contract  No.  DAAA09-76-C-2013,  from  the  Rock  Island 
Arsenal.  The  applications  of  Chapter  five  and  six  were  completed  under 
Research  Contract  No.  DAAK11-77-C-0023,  from  the  Ballistics  Research 
Laboratory  at  Aberdeen  Proving  Ground,  Maryland.  The  complete  results 
are  presented  in  this  report,  for  convenience  and  ease  of  use. 


CHAPTER  1 


INTRODUCTION 

1.1  Purpose  and  Scope  of  the  Study 

The  purpose  of  this  report  Is  to  present  an  optimal  structural 
design  method  that  integrates  the  concept  of  substructuring  into  its 
formulation  in  an  effective  way,  and  to  suggest  extensions  of  the 
method  to  fail-safe  design  of  structures.  While  substructuring  is 
commonly  used  in  the  analysis  of  large  structures,  no  general  method 
of  optimization  of  such  structures  that  incorporates  substructuring 
has  appeared  in  the  literature.  Therefore,  an  efficient  method  for 
design  sensitivity  analysis  with  substructuring  is  developed  using  a 
state  space  formulation  in  which  design  and  state  variables  are  treated 
separately  and  explicitly.  The  displacement  method  of  structural 
analysis,  an  efficient  design  sensitivity  analysis,  and  the  Kuhn-Tucker 
necessary  conditions  of  nonlinear  programming  are  employed  to  develop 
an  iterative  optimal  design  algorithm.  Further,  the  overall  efficiency 
of  the  algorithm  is  enhanced  through  carefully  coordinated  use  of 
design  variable  linking,  the  c-active  constraint  concept,  and  a finite 
element  analysis  method  that  is  organized  with  the  design  optimization 
task  in  mind. 

Incorporation  of  the  substructuring  concept  in  analysis .essentially 
means  partitioning  of  the  structure  into  substructures.  Using 
this  concept,  the  nodal  coordinates  are  partitioned  into  two  subsets: 

1)  coordinates  common  to  various  substructures,  defined  as  boundary 
coordinates,  and  2)  coordinates  that  occur  only  at  interior  nodes  of 
substructures,  defined  as  interior  coordinates.  The  equilibrium 
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equations  for  the  structures  are  partitioned  accordingly.  The  present 
algorithm  exploits  these  partitioned  equations,  both  in  structural 
analysis  and  in  the  design  sensitivity  analysis,  in  an  efficient  way. 

To  illustrate  applicability  of  the  method  and  to  study  its 
efficiency,  several  structures  of  given  configuration  are  optimized. 
These  include  both  two  and  three  dimensional  truss,  idealized  wing, 
and  frame  structures.  Many  of  these  structures  have  been  used  in  the 
literature  to  compare  optimization  techniques  [1-8],  Constraints  on 
stresses  and  displacements  under  multiple  loading  conditions,  on 
natural  frequency,  and  on  design  variables  are  imposed.  First,  results 
obtained  with  the  present  algorithm  for  truss  examples  are  compared 
with  results  obtained  without  substructuring  [1].  It  is  shown 
that  the  algorithm  with  substructuring,  especially  for  moderate  to 
large  scale  structures,  is  more  efficient  than  an  algorithm  without 
substructuring.  It  is  shown  that  the  present  method  of  state  space 
design  sensitivity  analysis  and  optimization  is  considerably  faster 
than  results  presented  in  the  literature  [2-8]. 

In  fail-safe  design  of  structures,  one  includes  projected  damage 
conditions  in  the  design  procedure.  These  damage  conditions  may 
include  damage  to  a set  of  members  of  structures  or  to  some  joints. 

The  substructuring  concepts  are  highly  significant  in  fail-safe  design, 
because  they  would  allow  one  to  treat  a large  number  of  damage 
conditions  very  efficiently.  The  reason  for  this  high  efficiency  is 
that  structural  analysis  and  design  sensitivity  analysis  for  damaged 
structures  will  be  very  efficient  when  substructuring  concepts  are 
employed.  This  efficiency  results  because  if  damage  to  the  structure 
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occurs  only  in  a substructure,  then  structural  matrices  for  only 
that  substructure  are  modified  and  need  to  be  recomputed. 

1.2  Review  of  Literature 

As  pointed  out  earlier,  a general  procedure  for  optimal  structural 
design  with  substructuring  has  not  been  developed  in  the  literature. 

A paper  by  Kirsch,  Reiss,  and  Shamir  [9]  is  the  only  exception  in 
which  the  problem  of  optimum  design  of  a structure  by  partitioning  it 
into  a number  of  substructures  is  considered.  Their  procedure  is  to 
sequentially  optimize  each  substructure,  with  respect  to  its  own  design 
variables,  while  keeping  other  design  variables  constant.  This  itera- 
tive procedure  is  inefficient  and  can  lead  to  non-optimal  solutions, 
as  pointed  out  by  the  authors. 

A procedure  for  calculating  derivatives  of  nodal  displacements  by 
formal  differentiation  of  various  substructure  equations  is  developed 
by  Noor  and  Lowder  [10].  This  procedure  will  be  required  when 
the  optimal  design  problem  is  defined  in  design  space  alone.  How- 
ever, when  the  problem  is  defined  in  a state  space  setting,  advantages 
of  self-adjointness  of  the  structural  problem  can  be  fully  brought  to 
bear  to  develop  efficient  design  sensitivity  analysis  procedures. 

Applications  of  optimization  theory  in  structures  have  been  a 
subject  of  great  Interest  during  the  past  two  decades.  It  has  been 
recognized  that  in  order  for  an  algorithm  to  be  applicable  to  large 
structures,  such  as  aircraft,  multi-story  buildings,  and  transmission 
towers,  it  must  be  computationally  efficient  and  compact.  It  must  be 
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executed  on  existing  computer  hardware,  without  exceeding  available 
memory  capacity,  and  be  computationally  fast,  i.e.,  the  time  spent  in 
each  redesign  cycle  is  kept  to  a minimum.  A very  detailed  review  of 
literature  on  optimal  design  has  appeared  in  papers  by  Wasiutynski 
and  Brandt  [11],  Sheu  and  Prager  [12],  and  Arora  and  Haug  [1].  For 
extensive  bibliographies  covering  various  aspects  of  optimal  design 
and  an  evaluation  of  the  state-of-the-art  until  1976,  the  reader  is 
referred  to  these  articles.  Also  recent  books  by  Schmit  [13],  Pope 
and  Schmit  [14],  and  Gallagher  and  Zienkiewicz  [15]  summarize  the 
state-of-the-art  quite  nicely. 

Two  philosophically  different  approaches  to  optimization  have  been 
advocated:  1)  Optimality  Criteria  Methods,  and  2)  Search  Methods 
based  on  gradient  calculations.  Methods  based  on  optimality  criteria 
are  died  indirect  and  employ  repeated  application  of  some  optimality 
criteria,  or  necessary  condition  for  optimality.  Many  researchers 
[3-7,16]  have  recently  presented  algorithms  and  results  based  on  such 
methods. 

Many  gradient  based  methods  are  also  available  in  the  literature. 
Among  these  are  the  Sequential  Unconstrained  Minimization  Techniques 
(SUMT)  [17]  used  by  many  researchers  [2,18-21],  the  Steepest  Descent 
methods  of  nonlinear  programming  [22-25],  and  the  State  Space  Gradient 
Projection  Method  developed  initially  by  Bryson  and  co-workers  [26] 
for  optimal  control  and  later  adapted  for  a wide  range  of  structural 
and  mechanical  design  problems  [1,27-33].  These  approaches  are 
collectively  referred  to  as  mathematical  programming  methods.  They 
all  rely  on  gradient  calculations,  each  employing  this  gradient  infor- 
mation in  its  own  iterative  algorithm. 
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Most  recently,  Schmit  and  Mlura  [2,18],  Dobbs  and  Nelson  [5], 

Rizzl  [3,4]  and  Berke  and  Khot  [6]  have  stated  the  structural  optimiza- 
tion problem  in  design  space.  On  the  other  hand,  Haug,  Arora,  and 
co-workers  [27-33]  have  pursued  a different  approach  in  defining  the 
problem.  In  their  approach,  the  problem  is  defined  in  a state  space 
setting,  i.e.,  the  state  and  design  variables  are  treated  separately 
and  explicitly  in  the  problem.  With  the  latter  approach  the  design 
sensitivity  analysis  is  performed  very  efficiently. 
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CHAPTER  2 

PROBLEM  DEFINITION  AND  DERIVATION  OF  THE  METHOD 

2.1  Introduction 

In  this  chapter,  a method  for  optimal  structural  design  using 
substructuring  is  developed.  Constraints  are  imposed  on  member 
stresses,  nodal  displacements,  natural  frequency,  and  design  variables. 

\ 

In  the  following  development  first  a general  optimal  design  problem  is 
defined.  Since  substructuring  is  anticipated  in  the  present  work, 
the  general  problem  is  then  reformulated  by  introducing  the  concept  of 
boundary  and  interior  coordinates.  Hereafter,  this  problem  is  called 
the  "Optimal  Design  Problem  with  Substructuring  - ODPS".  Design 
sensitivity  analysis  for  the  ODPS  is  developed  and  then  an  algorithm 
is  presented  in  a convenient  step-by-step  format. 

2.2  State  Space  Definition  of  an  Optimal 
Design  Problem 

A general  optimal  design  problem  for  structures  in  the  state 
space  setting  may  be  defined  as  follows:  find  a design  variable  vector 
b that  minimizes  a cost  function 


1*1  (b , z , £) 


(2.2-1) 


satifylng  the  equilibrium  equations  (state  equations) 
K(b)z  - S(b) 


' -"So*’-':— 


(2.2-2) 

| 


K(b)y  - CM(b)y 


(2.2-3) 


and  subject  to  the  constraints 

i 

♦8(b.«>  * 0 
d>d  (b)  < 0 

and 


(2.2-4) 

(2.2-5) 


4>e(C)  * 0 (2.2-6) 

where 

b ■ vector  of  D design  parameters,  such  as  cross-sectional 
areas,  moment  of  inertia,  thickness,  etc.  , 
z ■ state  variable  vector  of  N nodal  displacements, 

N “ number  of  degrees  of  freedom  of  the  structure, 

K(b)  ■ (N,N)  structural  stiffness  matrix, 

M(b)  - (N,N)  structural  mass  matrix, 

S(b)  ■ vector  of  N nodal  loads  on  the  structure, 

y ■ an  eigenvector  of  nodal  displacements,  corresponding 
to  the  lowest  eigenvalues  C,  , 

C - the  lowest  eigenvalue  of  Eq.  (2.2-3). 

Further  note  that  <j>e(?)GR1,  $S(b,z)6R  1,  and  <|>d(b)ER  2.  Therefore, 
the  total  number  of  constraints  is  c - 1 + + c2>  and  the  inequality 

applies  to  each  component  of  the  vector  functions  defined  in  Eqs. 
(2.2-4)  and  (2.2-5). 
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The  cost  function  of  Eq.  (2.2-1)  Is  quite  general  and  may  depend 
upon  weight  of  the  structure,  displacements  of  critical  points,  certain 
critical  member  forces,  or  perhaps  natural  frequency  of  the  structure. 
The  cost  function  depends  only  on  the  design  variables  If  it  represents 
weight  of  the  structure.  The  Inequality  (2.2-4)  represents  constraints 
that  depend  upon  state  and  design  variables.  These  are  the  member 
stress  and  the  nodal  displacement  constraints.  The  inequality  (2.2-5) 
represents  constraints  that  depend  only  on  design  variables.  These 
may  represent  either  explicit  bounds  on  the  design  variables  or  any 
relationship  between  them.  The  inequality  (2.2-6)  represents  a con- 
straint on  the  lowest  eigenvalue  of  Eq.  (2.2-3)  (4  > 4 ; 4 ■ allowable 

- o o 

lowest  eigenvalue)  which  is  related  to  the  fundamental  frequency  of  the 
structure  (frequency  f ■ /c/2tt  Hertz).  In  the  present  work,  constraints 
on  only  the  lowest  eigenvalue  are  considered,  but  constraints  on 
higher  eigenvalues,  if  desired,  can  also  be  included  [32]. 

It  should  also  be  noted  here  that  the  lowest  eigenvalue  of  Eq. 
(2.2-3)  may  also  be  related  to  the  buckling  load  for  the  structure. 

In  that  case,  the  matrix  M will  be  a geometric  stiffness  matrix  for 
the  structure  [34]  and  Eq.  (2.2-6)  will  represent  a constraint  on  the 
buckling  load  for  the  structure.  This  type  of  a constraint  can  be 
treated  in  the  algorithm  developed  herein. 

It  may  be  pointed  out  here  that  Schmit  and  Miura  [2],  Rizzi  [3], 
Dobbs  and  Nelson  [5],  and  others  [6-8,16,18]  have  stated  the  optimal 
design  problem  in  purely  design  space.  Generally,  the  cost  function 
of  Eq.  (2.2-1)  and  constraints  of  Eq.  (2.2-4)  are  also  explicit 
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functions  of  the  state  variables  (nodal  displacements)  z,  which  renders 

them  as  implicit  functions  of  the  design  variable  vector  b.  This  re- 

3 z 

quires  that  the  matrix  — be  explicitly  calculated  during  sensitivity 
analysis.  For  structural  problems  this  approach  is  highly  inefficient 
and  is  not  used  in  the  present  work. 

In  order  to  incorporate  the  substructuring  concept,  the  optimal 
design  problem  must  be  reformulated  in  terms  of  the  boundary  and 
interior  coordinates.  The  next  section  presents  structural  analysis 
by  substructuring,  before  the  problem  is  redefined  in  terms  of  new 
state  variables. 

2.3  Structural  Analysis  by  Substructuring 


2.3.1  Static  Analysis 

In  structural  analysis  using  the  substructuring  concept,  the 
stiffness  matrix  K,  and  the  vectors  z and  S are  partitioned  and  Eq. 
(2.2-2)  is  written  as  follows: 


kbb 

hT 

ZB 

V 

*IB 

V 

*I_ 

_SI_ 

(2.3-1) 


where  the  subscripts  B and  I refer  to  boundary  and  Interior  quantities 
for  all  substructures,  z and  S_GRn,  z_  and  S_GRm,  n is  the  number  of 

D D II 

boundary  degrees  of  freedom,  and  m - N - n is  the  number  of  interior 

St  , Sg,  etc.,  have 


degrees  of  freedom.  Submatrices  such  as 


compatible  dimensions  and  will  be  understood  to  be  functions  of  the 
design  variable  b. 
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Now,  the  interior  displacements  z are  eliminated  from  Eq.  (2.3-1) 
and  the  following  reduced  equation  is  obtained: 

kbzb  ■ fb  (2.3-2) 


where 


S “ *SiB  + (2.3-3) 

(n,n) 


FB  = SB  + QTSI  (2‘3^) 

(n.l) 


and 


Q - 

(m,n) 


K 


IB 


(2.3-5) 


Here,  ^ is  a boundary  stiffness  matrix  for  the  entire  structure,  and 
FgGR11  is  the  vector  of  effective  boundary  forces.  It  should  be  noted 
that  efficient  decomposition  procedures  are  used  to  decompose  K^j.  and 
then  solve  for  Q in  Eq.  (2.3-5).  For  more  computational  details  Chapter 
3 should  be  consulted.  The  boundary  stiffness  matrix  and  the  effec- 
tive boundary  force  vector  Fg  are  synthesized  by  considering  contribu- 
tions from  all  substructures.  For  this  purpose,  the  equilibrium 
equations  for  a substructure,  which  is  considered  as  an  isolated  free- 
body,  are  also  expressed  in  the  following  partitioned  form: 


r 


:y 


r • 


(2.3-6) 


where  the  superscript  r refers  to  the  rth  substructure  and  subscripts  B 

and  I refer  to  boundary  and  interior  quantities,  respectively.  Let 

n(r)  and  m(r)  represent  the  number  of  boundary  and  interior  coordinates 

of  the  rth  substructure,  respectively.  It  may  be  noted  that 
L 

m ■ E m(r),  where  L is  the  total  number  of  substructures.  Then, 
r-1 

- [n(r),  n(r)  ] , - [m(r) , n(r)],  - [n(r),  ra(r)], 

- Mr),  m(r)  ] , z<r)  and  sjr)GRn(r),  and  z<r)  and  sJr)GRm(r). 
From  the  second  line  of  Eq.  (2.3-6), 

{2.3-7) 

Now,  substituting  Eq.  (2.3-7)  into  the  first  line  of  Eq.  (2.3-6),  one 
obtains 


(2.3-8) 


where 
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[n(r),l] 


- sfr)  + 0(r>  S(r) 


(2.3-10) 


and 


Q(r) 

[m(r) ,n(r) ] 


- En  ] 1 K 


(r) 

IB 


(2.3-11) 


The  boundary  stiffness  matrix  K^r)  and  the  effective  boundary  force 
vector  Fg1^  for  each  substructure  are  computed  from  Eqs.  (2.3-9)  and 
(2*3-10),  respectively.  Finally,  Kg  and  Fg  are  assembled  according  to 
the  following  equations: 


and 


E icir^ 

r-1  ® 


(2.3-12) 


B 


SB  + 


L 

E 

r-1 


B(r)T  Q(r)T  S<*> 


(2.3-13) 


where  8(r)  is  a Boolean  transformation  matrix  of  dimension  [ n (r) ,n] . 

Using  the  reduced  equilibrium  equation  (2.3-2),  the  boundary 
displacements  zg  are  computed  by  a suitable  numerical  procedure.  The 
interior  displacements  are  then  computed  for  each  substructure,  using 
Eq.  (2.3-7).  Lastly,  element  forces  for  the  rth  substructure  are 
computed  from 


p(t)  , K(r)  z(r) 


(2.3-14) 


■ 

1- 


— ~r  - 


' .*  ' 


I 

! 
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where  p^r^  is  a vector  of  all  element  forces,  Is  a stiffness 

matrix,  and  is  a vector  of  nodal  displacements.  This  completes 

static  analysis  of  a structure,  using  the  substructuring  formulation. 

2.3.2  Frequency  Analysis 

The  natural  frequency  of  a structure  is  computed  from  Eq.  (2.2-3). 
A number  of  techniques,  such  as  Subspace  Iteration  [35],  Jacobi 
Iteration  [35],  Householder’s  Method  [36],  and  a method  based  on  Sturn 
Sequence  properties  described  by  Gupta  [37],  Wilkinson  [38]  and  others 
[39,40,41]  are  available  in  the  literature  for  solution  of  the  general 
eigenvalue  problem  defined  by  Eq.  (2.2-3).  However,  these  techniques 
require  computation  of  stiffness  and  mass  matrices  for  the  entire 
structure  and  decomposition  of  the  stiffness  matrix.  For  a large 
structure,  this  is  not  desirable,  since  it  defeats  the  purpose  of  sub- 
structuring. The  technique  of  frequency  analysis  used  should  either 
take  full  advantage  of  substructuring  or  it  should  not  require  calcula- 
tion and  storage  of  matrices  K and  M for  the  entire  structure.  There 
are  many  component  mode  substitution  techniques  available  in  the  liter- 
ature that  may  be  used.  For  a complete  survey  of  such  techniques,  the 
reader  is  referred  to  Arora’s  paper  [42].  These  techniques  take  full 
advantage  of  substructuring.  However,  they  are  not  suitable  for  inte- 
gration into  an  optimal  design  algorithm,  because  they  Introduce 
design  variable  dependent  transformations  that  complicate  the  design 
sensitivity  analysis. 

There  is  one  technique,  however,  that  has  all  the  features  needed 
to  solve  the  present  problem.  It  is  based  on  minimization  of  the 
Rayleigh  Quotient 
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R(y) 


T 

Ox 

T 

y My 


(2.3-15) 


subject  to  appropriate  constraints  on  orthogonality  of  eigenvectors. 
This  approach  of  finding  the  eigenvalues  and  the  corresponding  eigen- 
vectors has  been  discussed  and  used  successfully  by  many  researchers, 
such  as  Fox  and  Kapoor  [19],  Wilkinson  [38],  and  Bradbury  and  Fletcher 
[43].  Further  it  has  been  concluded  by  these  researchers  that  the 
method  of  conjugate  gradients  is  quite  suitable  for  the  minimization 
of  the  function  of  Eq.  (2.3-15),  because  it  is  efficient  and  requires 
minimum  additional  computer  storage.  Also,  the  method  does  not  re- 
quire generation  and  storage  of  the  matrices  K and  M for  the  entire 
structure,  because  all  calculations  can  proceed  elementwise.  For  more 
details  the  reader  is  referred  to  Refs.  19  and  43. 

There  is,  however,  one  difficulty  with  this  procedure  of. computing 
eigenvalues.  Convergence  to  an  eigenvalue  and  the  corresponding 
eigenvector  can  be  quite  slow  if  a good  initial  guess  of  eigenvector 
is  not  known.  Some  methods  of  selecting  initial  eigenvectors  have 
been  suggested  [19,43,44],  but  no  general  procedure  exists  to  alleviate 
this  problem.  Therefore,  in  the  present  work,  an  approximate  eigen- 
value problem  is  solved  to  obtain  a good  starting  point  for  the 
conjugate  gradient  method.  The  approximate  eigenvalue  problem  is 
defined,  using  the  procedure  of  Refs.  34,  44,  45,  and  46,  as  follows: 


*ByB  " CBMByB 


(2.3-16) 


where  is  a (n,n)  condensed  mass  matrix,  is  an  eigenvalue,  and 
yB£Rn  is  a vector  of  generalized  coordinates  associated  with  the 
boundary  degrees  of  freedom.  The  procedure  for  defining  the  approxi- 
mate eigenvalue  problem  is  summarized  in  Appendix  I,  for  easy  reference. 

The  approximate  eigenvalue  problem  of  Eq.  (2.3-16)  may  be  solved 
by  any  of  the  previously  mentioned  methods,  because  and  its  decompo- 
sition are  available  from  static  analysis  of  the  structure.  In  the 
present  work,  the  method  of  Subspace  Iteration  [35],  which  has  been 
used  quite  successfully  [1,30,47],  is  employed.  Once  y is  known,  the 
approximate  eigenvector  for  the  entire  structure  is  computed  from  the 
following  equation  (Appendix  I) : 

VB  (2.3-17) 

Thus,  in  the  proposed  method  of  frequency  analysis,  the  approximate 
problem  of  Eq.  (2.3-16)  is  solved  by  Subspace  Iteration  [35]  only  in 
the  first  iteration.  In  all  subsequent  iterations,  the  method  of 
conjugate  gradients  is  used  to  minimize  the  Rayleigh  Quotient  to  obtain 
the  lowest  eigenvalue  and  the  corresponding  eigenvector. 

2.4  Statement  of  the  Optimal  Design  Problem 
with  SubstructurinR  (ODPS) 

In  this  section,  the  optimal  design  problem  of  Section  2.2  is 

redefined  in  terms  of  new  state  variables  z_  and  z_,  as  follows:  find 

D X 

a design  variable  vector  b that  minimizes  the  cost  function 


r i y— — k - 


■*-y- 


l .A  ** 
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J = Jtb.Zg.Zj.O 


(2.4-1) 


satisfies  the  equilibrium  equations  (state  equations) 


Vb  - fb 


(2.3-2R) 


W • ^r> 


(r)  (t)  r * 1 to  L 

IB  B 


(2.3-7R) 


Ky  * CMy 


(2.2-3R) 


and  satisfies  the  constraints 


♦ (b.Zg.Zj)  < 0 


(2.4-2) 


$a(b)  < 0 


(2.2-5R) 


*e(c)  < o 


(2.2-6R) 


where  ♦'(b.Zj.^ei'1,  AbiG^2  and  /(OGR1.  For  computational 
convenience  and  efficiency,  it  is  desirable  to  divide  the  member 
forces  and  nodal  displacements,  for  which  design  sensitivity  informa- 
tion is  required,  into  two  subsets:  1)  a subset  consisting  of  member 
forces  and  nodal  displacements  that  explicitly  depend  upon  the  design 
parameters  b,  boundary  displacements  Zg,  and  interior  displacements 
z ; and  2)  a subset  consisting  of  member  forces  and  nodal  displace- 
ments that  depend  on  design  parameters  b and  the  boundary  displacements 

z only.  Thus,  the  constraints  of  Eq.  (2.4-2)  are  rewritten  as. 

B 
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, si 


(h’VV  - 0 


(2.4-3) 


<?82(b,zB)  < 0 


(2.4-4) 


where  ^(b.Zj.J^GR  1 and  4is2(b,zg)eR 
Terms  which  are  used  frequently  in 
defined. 


, such  that  c^  = + y^. 

subsequent  chapters  are  now 


Def/.  The  set  of  points  b£R^  that  satisfies  constraints  of  Eqs. 
(2.2-5),  (2.2-6),  (2.4-3),  and  (2.4-4)  is  called  the  constraint  set 
for  ODPS. 


Def : A constraint  4^  < 0 is  said  to  be  active  if  4^  > 0. 

Def : A constraint  4>i  < 0 is  said  to  be  e-active  if  4>  + e > 0, 
where  e is  a small  positive  number. 

It  should  be  noted  here  that  assumptions  on  continuous  differ- 
entiability of  cost  and  constraint  functions,  linear  independence 
of  constraint  gradients,  and  closedness  and  boundedness  of  the  con- 
straint set  are  made  in  the  present  work. 

2.5  Design  Sensitivity  Analysis  with 
Sub structuring 

The  philosophy  of  an  iterative  technique  is  to  start  with  the 
best  engineering  estimate  of  the  design  variable  vector  b and  to 
improve  it  until  an  optimum  is  reached.  One  must  determine  the  effect 
of  a design  change  on  the  cost  and  constraint  functions  before  a 


# 
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design  improvement  fib  can  be  calculated.  This  is  known  as  design 
sensitivity  analysis.  In  this  section,  design  sensitivity  matrices 
for  cost  and  constraint  functions  are  derived.  In  the  derivation 
first  order  Taylor's  expansion  and  certain  adjoint  relationships  are 
employed.  First,  the  design  sensitivity  analysis  of  a general  function 
is  considered.  Then  the  analysis  is  specialized  to  the  ODPS,  defined 
in  Section  2.4. 


2.5.1  General  Approach 

Let  iKb, Zg.z^.,  ?)e  Ry  represents  a general  vector  function.  The 

function  i|»  will  be  scalar  if  it  represents  the  cost  function.  When  the 

design  is  changed  by  a small  amount  fib,  the  displacements  z and  z , 

B I 

and  the  eigenvalue  C will  also  change  by  small  amounts  fiz  , fiz  , and 

B I 

fiC;  due  to  the  well  posed  nature  of  the  structural  analysis  problem. 

Let  6i|>  represent  a first  order  change  in  the  function  i Then,  taking 

b,  z , z , and  c as  independent  variables,  6 \J>  is  given  as 
B 1 


" lb  6b  + It  6zb  + I h 6zi  + It  65 

o 1 


(2.5-1) 


where  the  derivatives  , and  f-j;  are  computed  at  the 

ob  dZg  o Zj  ot, 

previously  known  values  of  b,  z_,  z , and  C-  The  problem  is  now  to 

O I 


express  fizB,  fiz  , and  6c  in  terms  of  fib,  so  that  fiip  in 


9 ip 
3C 


dzn  B’  3z  I’ 

D X j | 

Eq.  (2,5-1)  is  expressed  as  -jr  6b.  Let  us  first  consider  the  terms 

ab 

dll/  dll/ 

6z^  and  6z^.  In  order  to  obtain  these  expressions  in  terms 

B I 

of  fib,  first  order  expansions  of  the  partitioned  state  Eq.  (2.3-1), 


are  written  as: 
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Sb  6zb  + Si  6zi  = ci  6b  (2.5-2) 

KIB  6zB  + KII  6zi  = c2  6b  (2.5-3) 


where 


and 


dS 

C1  " 3b  3b  (KBB  V “ 3b  (KBI  ZIJ 
(n,D) 


(2.5-4) 


C2  ” 3b  " 3b  (KIB  V " 3b  (KII 
(m,D) 


ZI} 


(2.5-5) 


These  derivatives  are  evaluated  at  the  known  values  of 
Using  Eq.  (2.5-3),  one  eliminates  6zj  from  Eq.  (2.5-2) 


b,  2^1  and  Zj. 
to  obtain 


Kb  6zb  * C6b 


(2.5-6) 


where 


C * C + QTC2  (2.5-7) 

(n,D)  1 


Equation  (2.5-6)  is  simply  a first  order  expansion  of  Eq.  .(2.3-2). 
Now,  define  an  adjoint  matrix  AjOn.y)  to  be  a solution  of 


K * - 

II  I 3zx 


(2.5-8) 


f 


20 


i 

Ik 


Taking  the  transpose  of  Eq.  (2.5-8),  postmultiplying  by  fiZj,  and  sub- 
stituting from  Eq.  (2.5-3),  one  obtains 

3i 


VC2  Sb  ' *1.  62B>  ‘ SzI 


(2.5.9) 


The  right  hand  aide  of  Eq.  (2.5-9)  is  exactly  the  term  needed  in  Eq. 
(2.5-1),  so  making  a substitution,  one  obtains 


J 


6ip 


d + C2]  6b  + (ft  - Kib]  6zB  + H 6C  (2-5-10) 


Another  adjoint  matrix  Ag(n,y)  is  defined  to  be  a solution  of 


S' 


T 

. lil . 
B 3zB 


‘Si  xi 


(2.5-11) 


Substituting  for  Aj  from  Eq.  (2.5-8)  and  making  use  of  Eq.  (2.3-5), 
one  reduces  Eq.  (2.5-11)  to 


3z, 


(2.5-12) 


Thus,  Ag  is  a solution  of  Eq.  (2.5-12).  Taking  the  transpose  of  Eq. 
(2.5-11),  postmultiplying  by  6Zg,  and  using  Eq.  (2.5-6),  one  obtains 


XB  C6b  “ [It  ' Xl  KIb]  6J 


(2.5-13) 


The  right  hand  side  of  Eq.  (2.5-13)  is  exactly  the  term  needed  in  Eq. 
(2.5-10),  so  making  a substitution,  one  obtains 


T- — - 


*♦  * . .a*** 

* •* 

■ i i 
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[Jb  + xi  c2  + xb  c]  Sb  + n 


(2.5-14) 


Now  one  must  treat  the  expression  64.  The  design  sensitivity 
analysis  of  the  eigenvalue  C has  been  considered  by  many  researchers 
[19,27,30].  Therefore,  this  development  is  only  summarized  here.  From 
the  first  order  expansion  of  Eq.  (2.2-3)  and,  using  the  fact  that  K(b) 
and  M(b)  are  symmetric,  one  obtains  the  following  expression  for  6s : 

65  “ TyT  {lb  [K(b)yl  ” C |b  [M(b)yl  }~\  (2.5-15) 

Substituting  this  expression  for  6s  in  Eq.  (2.5-14),  one  obtains 


6\p  * G 6b 


(2.5-16) 


where 


G . il_  + CT  x T + cT  xD  + xe 

x 3b  2 1 B 

(D,y) 


(2.5-17) 


T X 

>e  - [lb  {K<Wy}  ‘ }]  y |f-/xTM0»)y  <2-5-18> 

Equation  (2.5-16)  is  the  desired  relationship  between  the  design  change 
and  the  changes  in  member  forces,  nodal  displacements,  cost  function, 
and  eigenvalue.  Columns  of  the  matrix  G are  the  required  design 
sensitivity  vectors. 


- ■v’s 


* 
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In  Refs.  2,  3,  5,  and  16,  the  optimal  design  problem  is  formulated 
in  design  space  only.  Thus,  a design  sensitivity  matrix  G correspond- 
ing to  Eq.  (2.5-17)  must  be  expressed  as 


G 


h , n 9zb  , »»  8zi 

3b  3z  3b  3z  3b 

— D 1 


M li  1 
3S  3b  J 


(2 . 5-17a) 


This  method  of  computing  the  matrix  G requires  explicit  calculations 

of  matrices  -rr~  , rr—  , and  rr  . Computation  of  — f may  be  calculated 
3b  3b  3b  , ..  3b 

3ZB 

as  in  Eq.  (2.5-18).  Calculation  of  — — and  -rr—  , however,  is  accom- 

3b  3b 

plished  by  formally  differentiating  the  state  equations  (2.3-2)  and 
(2.3-7),  respectively. 

It  is  shown  in  Ref.  10  that  the  equation  for  calculation  of  the 

8zb 

matrix  -rr~  is  similar  to  Eq.  (2.5-12),  but  with  a different  right  hand 

9*r  92B 

side.  The  equation  for  calculation  of  ~r—  can  be  derived  using  tt— - , 

do  db 

and  is  similar  to  Eq.  (2.5-8),  but  with  a different  right  hand  side. 

Some  of  the  calculations  of  this  approach  are  similar  to  the  state  space 
approach,  but  the  major  difference  is  in  the  number  of  vectors  ir.  the 
right  hand  side  of  Eqs.  (2.5-8)  and  (2.5-12)  and  in  the  number  of 
calculations  required  to  generate  the  right  hand  side.  In  the  state 
space  formulation,  the  right  hand  sides  of  these  equations  are  generated 
quite  readily.  Also,  the  number  of  vectors  in  the  right  hand  side  of 
Eq.  (2.5-8)  is  equal  to  the  number  of  active  constraints  that  depend 
upon  any  component  of  the  vector  z^.  Similarly,  the  number  of  vectors 
in  the  right  hand  side  of  Eq.  (2.5-12)  is  equal  to  the  total  number  of 
active  constraints.  In  the  purely  design  space  formulation,  the 
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number  of  vectors  in  the  right  hand  side  of  Eqs.  (2.5-8)  and  (2.5-12) 
is  equal  to  the  interior  and  boundary  degrees  of  freedom  for  the  entire 
structure.  This  number  is  generally  much  larger  than  the  number  of 
active  constraints.  Thus,  the  state  space  formulation  should  be  more 
efficient. 


2.5.2  Reduction  of  Design  Sensitivity 
Matrices  for  the  ODPS 

In  the  following  development,  let  a - over  a constraint  function 
represent  inclusion  of  only  the  active  constraints.  Linear  expansions 
of  the  cost  function  of  Eq.  (2.4-1)  and  the  constraints  of  Eqs. 

(2.4-3) , (2.4-4),  (2.2-5),  and  (2.2-6),  with  the  procedure  outlined 
above,  yield  the  following  expressions: 


(2.5-19) 


6$Sl  “ A®1  6b 

6$S2  - A®2  6b 
T 

64>d  = Ad  6b 

and 


(2.5-20) 

(2.5-21) 

(2.5-22) 
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where  the  vector  A and  other  A matrices  are  obtained  from  Eqs. 
(2.5-15),  (2.5-17),  and  (2.5-18),  as  follows: 


AJ  3J  . _T  . J . T . eJ 

A 3b”  + C2  Xl  + C *B  + A 


(2 


si  _ T si  T si 

3b  C2  I C B 


(2 


's2  ■ ♦ <=T 


(2 


.14- 


3b 


(2, 


aC  “ ft-  [lb  {K(b)y}  " *fb  {M(b)y  0 y/yVb)y  (2 


and 


.eJ 


If  [lb  {K(b)y}  - 5 lb  {M(b)y}]  y/yTM(b)y 


(2, 


The  vectors  X!?  and  Xd,  and  the  matrices  X^,  X*^  , and  X?^  are 

ib  IBB 

solutions  of  the  following  equations: 


K X?  - ti- 
ll I 3z, 


(2. 


if  iJ  3jT  j.  nT  3jT 
*B  B "^+Q 


(2. 


.5-24) 


.5-25) 


.5-26) 


5-27) 


5-28) 


5-29) 


5-30) 


5-31) 
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kt7  »■!  . if!i 

II  I 3 z. 


(2.5-32) 


KB  a 


8l 

B 


^ + qTi, 


si 


B 


3z, 


(2.5-33) 


and 


*8  X 


s2 

B 3z 


s2J 


B 


(2.5-34) 


2.6  Optimal  Design  Algorithm 

Now,  certain  restrictions  are  placed  on  the  linearized  constraint 
functions.  It  is  required  that  the  design  change  fib  be  computed  in 
such  a manner  that  it  corrects,  or  at  least  improves,  all  the  violated 
constraints.  These  requirements  on  Eqs.  (2.5-20)  to  (2.5-23),  can  be 
stated  in  the  form  of  the  following  inequalities: 


Asl  fib  < aJ>s1 


(2.6-1) 


A82  fib  < A*®2 


(2.6-2) 


Ad  fib  < A$d 


(2.6-3) 


I 


A®  fib  < A$® 


(2.6-4) 


• , y*  '..  ./■**«• 

* /*  • 
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***S-L  ~ 92  ~ q *»g 

where  A^  , A4>  , A<J>  , and  A$  are  desired  corrections  in  constraint 
violations.  If  a constraint  < 0 is  active,  then  A^  ■ with 

a^  > 0.  The  positive  multiplier  a^O  < a^  < 1.0)  is  chosen  by  the 
designer  and  is  simply  a scale  factor  for  the  amount  of  violation  to  be 
corrected.  In  many  instances,  it  may  be  desirable  to  call  for  only  a 
partial  correction  of  the  constraint  violation.  In  such  a case  a^  < 1. 
If  full  correction  of  the  constraint  violation  is  required  then  a^  = 1 . 

The  constraints  of  Eqs.  (2.6-1)  to  (2.6-4)  have  similar  forms,  so 
they  can  be  written  in  a compact  form  as 

AT  <Sb  £ A$  (2.6-5) 


where 


- |a81  »*2  / A6] 


(2.6-6) 


and 


y 


,T  _T  .T  T 
•si  ~s2  ~d  ~e 
A<p  * jA<p  A(p  A<ji 


]' 


(2.6-7) 


The  reduced  problem  of  computing  an  optimum  design  change  fib  can 
now  be  stated  as  follows:  find  fib  to  minimize  the  cost  function  of 
Eq.  (2.5-19),  subject  to  the  constraint  of  Eq.  (2.6-5)  and  the  step 
size  constraint 


T 2 

fib  Wfib  < C 


(2.6-8) 
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where  W is  a positive  definite  weighting  matrix  and  £ is  a small 
number.  The  matrix  W,  usually  diagonal,  is  used  to  assign  weights  to 
the  various  components  of  6b  and  is  often  necessary  when  components  of 
b represent  different  physical  quantities  of  different  orders  of 
magnitude.  The  reduced  ODPS  defined  above  is  exactly  the  problem 
defined  in  Refs.  1 and  27.  An  application  of  the  Kuhn-Tucker  condi- 
tions of  nonlinear  programming  gives  the  following  solution  [1,27]: 


6b  = neb1  + 6b2 

6b1  - -W_1  [AJ  + Ap1] 

.,2  -1,  2 

6b  = -W  Ap 

Hfp1! p2]  = [(-A1  W_1  AJ)  : -A*] 
T -1 

H = A W A 


(2.6-9) 

(2.6-10) 

(2.6-11) 

(2.6-12) 

(2.6-13) 


1 2 

p - p + (1/n)  P (2.6-14) 

where  n > 0 is  a step  size  to  be  chosen  by  the  designer  and  p > 0 is 
a Lagrange  multiplier  vector.  The  method  of  step  size  selection  is 
the  same  as  used  in  Refs.  27,  33,  and  41. 

The  method  can  now  be  described  by  the  following  step-by-step 
algorithm: 

Step  1.  At  the  jth  design  point  b^\  generate  the  matrices 
and  for  each  substructure.  Decompose  each  and 
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calculate  the  matrix  from  Eq.  (2.3-11).  Store  the  decomposed  part 

of  the  matrix  K^^?  and  the  matrix  Q^r^  for  later  calculations. 

Calculate  the  boundary  stiffness  matrix  Kg  and  the  effective  boundary 
load  vector  Fg  from  Eqs.  (2.3-12)  and  (2.3-13),  respectively.  Decom- 
pose the  matrix  Kg  and  store  it  for  later  use. 

Step  2.  Calculate  the  boundary  displacements  z from  Eq.  (2.3-2) 

D 

and  the  interior  displacements  for  each  substructure  from  Eq. 

(2.3-7). 


Step  3.  Calculate  the  lowest  eigenvalue  and  corresponding  eigen- 
vector from  Eq.  (2.3-15).  If  j-1,  then  an  approximate  eigenvector  for 
the  Rayleigh  Quotient  is  obtained  from  Eqs.  (2.3-16)  and  (2.3-17). 

Step  4.  Compute  the  vectors  XlJ  and  from  Eqs.  (2.5-30)  and 

1 15 

(2.5-31),  respectively.  Assemble  the  matrix  of  Eq.  (2.5-24). 


Step  5.  Check  the  frequency  constraint  of  Eq.  (2.2-6).  If  it  is 


violated,  then  compute  Ac  of  Eq.  (2.5-28),  and  put  A$e  = -a^ 


7 e 


Step  6.  Check  the  constraints  of  Eq.  (2.4-2)  and  form  the  vectors 
“*  sl  s 2 

4>  and  $ of  Eqs.  (2.4-3)  and  (2.4-4),  respectively.  Calculate  the 

3i81  3d>81  3*sl  3is2 

sensitivity  information  , — f — , and  — . Calculate  ■ ■*  , 

_o  ®h  3z  3z  3b 

^2®^  B I ^ . 2 

and  — * — simultaneously.  Also  calculate  A<f>8  and  A$8 

3zb 


Step  7.  Calculate  X8^r^  for  each  substructure  from  Eq.  (2.5-32), 

sl  s2 

Also,  calculate  the  matrices  X and  X from  Eqs.  (2.5-33)  and 

D D 

(2.5-34),  respectively.  The  latter  calculations  are  performed 
simultaneously. 


d. 
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SteP  8-  Calculate  the  matrices  C1  and  from  Eqs.  (2.5-4)  and 
^•5-5),  respectively.  Also,  calculate  the  matrix  C from  Eq. 

(2.5-7). 

SteP  9-  Assemble  the  matrices  Asl  and  As2  of  Eqs.  (2.5-25)  and 
(2.5-26),  respectively. 

■SJLe.P.  10 • Check  the  constraints  of  Eq.  (2.2-5)  and  form  a vector 
$d.  Compute  the  matrix  Ad  of  Eq.  (2.5-27).  Also,  compute  A<j>d . 

Step  11.  Finally,  assemble  the  matrix  A and  A$  of  Eqs.  (2.6-6) 
and  (2.6-7),  respectively. 

1 2 

-Step  12.  Compute  u and  y from  Eq.  (2.6-12).  Choose  a step 
size  n and  compute  the  Lagrange  multiplier  from  Eq.  (2.6-14). 

Step_13.  Check  the  sign  of  each  component  of  u.  If  any  component 
of  \i  is  negative,  remove  the  corresponding  row  from  AT  and  and 
return  to  Step  12. 

1 2 

Step  14.  Compute  6b  , 6b  , and  6b  from  Eqs.  (2.6-10),  (2.6-11), 
and  (2.6-9),  respectively.  Put  b^+1^  - + 6b  and  y^+1^  = 

Step  15.  If  all  constraints  are  satisfied  and 

,,1M  iT  i 1/2 

l|6b  ||2  = [6b  W 6b  ] is  sufficiently  small  [27],  terminate 
the  process.  Otherwise,  return  to  Step  1 with  b^+1^  as  the  best 
available  design. 
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CHAPTER  3 

DISCUSSION  OF  THE  METHOD  AND  COMPUTATIONAL  CONSIDERATIONS 

..  3.1  Introduction 

The  method  for  ODPS  of  Chapter  2 is  quite  general  and  can  be 
applied  to  any  type  of  structure,  because  no  assumption  is  made  re- 
garding the  type  of  structural  elements.  However,  details  of  applica- 
tion of  the  method  to  various  types  of  structures,  such  as  trusses, 
frames,  plates,  aircraft  wing  structures,  and  shells,  will  be  quite 
different.  This  is  due  to  the  fact  that  the  structural  analysis  and 
the  design  sensitivity  analysis  are  quite  different  for  different 
types  of  structures.  The  types  of  design  variables  for  different 
structures  are  also  of  a different  nature.  Therefore,  in  this  chapter, 
only  those  computational  aspects  of  the  method  which  are  common  to  all 
types  of  structures  are  discussed. 

3.2  Design  Variable  Linking 

Due  to  many  practical  design  considerations,  such  as  symmetry  in 
the  structure,  fabrication,  and  joint  connections,  it  Is  desirable  to 
constrain  some  elements  of  a structure  to  the  same  design  variable 
values.  This  is  also  known  as  grouping  of  elements  [30]  or  linking  of 
design  variables.  Elements  of  the  structure  are  divided  into  a number 
of  groups  and  one  design  variable  is  assigned  to  each  group  of  elements. 
Recently,  Pickett,  Rubinstein,  and  Nelson  [48]  have  tried  to  generalize 
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the  idea  of  reducing  the  number  of  design  variables  by  using  the  well 
known  reduced  basis  concept.  A reduced  set  of  design  variables  is 
introduced  by  the  transformation 

b - Ub  (3.2-1) 


where 


b = a vector  of  D original  design  variables, 

b = a vector  of  D reduced  original  variables, 

U = a matrix  (D,D)  whose  columns  are  the  base  vectors, 

and  D >>  D.  A number  of  techniques  have  been  suggested  [48]  for 
selecting  the  base  vectors,  including 

(i)  design  variable  linking, 

(ii)  optimality  criteria  solutions,  and 

(iii)  upper  and  lower  bound  solutions. 

Techniques  for  selection  of  base  vectors  for  a reduced  design  basis, 
however,  need  further  work  [1,48,49]  and  are,  therefore,  not  employed 
with  the  present  algorithm.  When  design  variable  linking  is  used, 
the  form  of  the  matrix  U is  quite  simple.  The  ith  row  of  U contains 
only  a unit  positive  element  at  the  jth  location,  indicating  that  the 
ith  element  is  linked  to  the  jth  design  variable.  In  the  case  where 
design  variable  linking  is  not  required,  U is  simply  the  identity 
matrix.  Design  variable  linking  can  be  readily  incorporated  in  the 
algorithm  of  Chapter  2. 
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3.3  Computational  Considerations  in 
Structural  Analysis 

Structural  analysis  is  an  important  part  of  the  optimal  design 
algorithm.  Therefore,  extra  effort  should  be  devoted  to  make  this 
portion  of  the  algorithm  as  efficient  as  possible.  In  this  section, 
structural  analysis  computations  relative  to  ODPS  are  briefly  described. 

In  static  analysis  of  structures,  the  response  variables  to  be 
determined  are  the  boundary  displacements  z , the  interior  displace- 

D 

ments  for  each  substructure  the  lowest  eigenvalue  £,  and  the 

element  stresses.  Before  these  response  quantities  can  be  computed, 
the  boundary  stiffness  matrix  Kg  and  the  interior  stiffness  matrices 
for  each  substructure  must  be  computed.  The  former  matrix  is 

computed  by  summing  contributions  from  each  substructure  (Eq.  (2.3-12)). 
The  nodes  of  a substructure  are  numbered  in  a manner  such  that  they 
give  smallest  bandwidth  for  K^j\  The  interior  stiffness  matrices 
are  decomposed  and  saved  for  later  use.  The  matrix  Q ' is  also 
computed  for  each  substructure  and  saved.  Contribution  to  the  effec- 
tive boundary  load  matrix  from  the  interior  load  matrix  for  each 
substructure  is  also  computed  here.  The  boundary  nodes  for  the  entire 
structure  are  numbered  in  a way  such  that  Kg  has  the  smallest  bandwidth. 
The  matrix  Kg  is  also  decomposed  and  saved  for  later  use.  The  boundary 
displacements  zfi  for  the  entire  structure,  under  all  loading  conditions, 
are  computed  from  Eq.  (2.3-2).  Then,  each  substructure  is  considered 
separately  and  the  Interior  displacements  are  computed  from  Eq. 

(2.3-7).  Finally  element  forces  p ' for  each  substructure  are 
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computed  from  Eq.  (2.3-14).  It  may  be  noted  here  that  the  element 
forces  are  computed  elementwise,  without  assembling  the  matrix  . 

In  frequency  analysis,  the  fundamental  frequency  of  the  structure 
Is  computed  by  minimizing  the  Rayleigh  Quotient  of  Eq.  (2.3-15).  Thus, 

\ 

the  problem  is  to  find  an  eigenvector  y such  that  the  Rayleigh  Quotient 


of  Eq.  (2.3-15)  is  minimized,  subject  to  one  of  the  following  normali- 
zation constraints: 

T 

(1)  The  Euclidean  norm  constraint,  l.e.,  y y = 1 

T 

(11)  The  M-norm  constraint,  l.e.,  y My  = 1 

(ill)  The  sum  (L-)  norm  constraint,  l.e.,  I |y  | «*  1,  and 

1 1 

(iv)  The  maximum  (LJ  norm  constraint,  i.e.,  max  |y  | - 1. 

i 1 

These  constraints  make  the  minimization  problem  for  Eq.  (2.3-15) 

constrained.  The  last  constraint  being  the  simplest,  is  employed  by 

many  researchers  [19,43].  In  this  approach,  max  J y j is  determined 

i 1 

at  the  first  iteration  of  the  minimization  procedure  and  the  vector 
is  scaled  so  that  this  component  is  unity.  The  corresponding  component 
of  gradient  of  the  Rayleigh  Quotient  is  set  to  zero.  The  minimization 
problem  is  then  treated  as  an  unconstrained  problem. 

The  first  three  approaches  can  also  be  used,  in  which  case  a 
component  of  the  eigenvector  is  found  from  the  equality  constraint  and 
eliminated  from  Eq.  (2.3-15).  The  reduced  minimization  problem  is  then 
unconstrained.  The  later  approaches,  however,  introduce  complications 
in  the  structural  problem,  because  matrices  K and  M will  have  to  be 
computed  and  stored.  Therefore,  in  the  present  approach,  constraint 
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(iv)  is  selected,  since  it  can  easily  be  programmed  as  an  integral 
part  of  the  finite  element  formulation  of  the  ODPS. 

The  Rayleigh  Quotient  and  its  gradient 


VR(y)  - 2(Ky  - R(y)My)/(y  My) 


(3.3-1) 


are  computed  by  considering  an  element  of  the  structure  at  a time. 

Thus,  all  calculations  for  frequency  analysis  proceed  elementwise. 

As  pointed  out  in  Chapter  2,  the  conjugate  gradient  method  is  used 
for  the  minimization  of  Rayleigh  Quotient.  In  this  method,  one  starts 
with  an  initial  estimate  y ^ of  the  eigenvector.  The  conjugate  gradient 
algorithm  is  described  by  the  following  iterative  equation: 


yi+i  ■ yi  + Vi 


i - 0,1,2,... 


(3.3-2) 


where 


si+l  “ " 8i+l  + Vi 


gi+l  = *R(yi+1) 


si+l' 


1 i8ir 


3 » “ 

O O 


go  - VR(yo) 


(3.3-3) 


(3.3-4) 


(3.3-5) 


(3.3-6) 


(3.3-7) 


w . - ••*'*’**■ 
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and  ot^  is  a step  size  for  the  ith  iteration.  In  determining  the  step 
size,  a quadratic  interpolation  formula  is  used.  The  necessary  condi- 
tion, i.e.,  -j^|—  [R(y^  + a^s^)]  **  0 , yields  the  following  quadratic 

expression  for  determining  [43]: 

Vi  + A2°i  + A3  " 0 (3.3-8) 

where 


A1  = (y^MSi)  - (y^KSi)(s^MSl)  (3.3-9) 

A2  = (s^KSi) (y^My±)  - (y^Ky^ (sjfa^  (3.3-10) 

A3  = (y^Ks1) (y^My±)  - (y^Kyf) (y^Msf)  (3.3-11) 

Two  roots  of  Eq.  (3.3-8)  correspond  to  the  maximum  and  minimum  values 
of  Rayleigh  Quotient  in  the  direction  s^.  Further  the  positive  and 
negative  roots  of  Eq.  (3.3-8)  correspond  to  the  minimum  and  maximum 
function  values  [43].  This  forms  a convenient  criterion  for  selecting 
the  value  of  a^.  It  may  be  pointed  out  here  that  most  computations  in 
the  Rayleigh  Quotient,  its  gradient,  and  the  step  size  can  proceed 
s imul t aneous ly . 

As  pointed  out  earlier,  a good  estimate  of  the  eigenvector  must 
be  available  for  fast  convergence  of  the  conjugate  gradient  method. 
Therefore,  the  approximate  eigenvalue  problem  of  Eq.  (2.3-16)  is  solved 
to  obtain  a good  starting  eigenvector.  The  Subspace  Iteration  method 
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[35],  used  for  solving  the  approximate  eigenvalue  problem,  requires 
computation  of  the  boundary  mass  matrix  Mg.  This  matrix  is  assembled 
in  the  same  way  as  the  matrix  Kg.  After  the  first  design  cycle,  the 
approximate  problem  of  Eq.  (2.3-16)  is  not  solved,  since  the  previous 
solution  provides  a good  initial  estimate  of  the  eigenvector  for  the 
subsequent  design  cycle.  Thus,  substructuring  is  not  required  for 
frequency  analysis  after  first  design  cycle.  Finally,  convergence 
criteria,  based  on  the  function  value,  norm  of  the  eigenvector,  and 
the  gradient  of  Rayleigh  Quotient,  are  used  [19,  38].  These  are 
expressed  as 

R(y1+1>  - R(y1) 

R(y±) 

I ly1+1  - yj  I i «• 2 0.3-13) 


S ei 


(3.3-12) 


and 


| | VR(y  ) | | < e (3.3-1/*) 

2 ~ 

where  > 0,  i * 1 to  3,  are  specified  small  numbers. 

3.4  Multiple  Loading  Conditions 
and  Constraint  Checks 

Most  structures  are  designed  to  withstand  multiple  loading 

conditions.  This  situation  is  handled  in  formulation  of  the  ODPS  by 

expanding  the  state  variable  vectors  z and  z to  include  all  states. 

o I 
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The  element  force  vector  Is  expanded  accordingly,  and  the  stress  and 
displacement  constraints  are  checked  under  each  loading  condition. 

The  number  of  constraints  In  ODPS  is  quite  large.  So,  an  effort 
should  be  made  to  logically  delete  constraints  which  will  not  be 
critical  at  the  optimum.  In  the  present  research,  the  idea  of  "worst 
violated  constraint"  is  used,  in  order  to  eliminate  redundant  con- 
straints. In  the  case  of  stress  constraints,  worst  violation  under 
all  loading  conditions  and  for  all  elements  of  a group  is  imposed, 
instead  of  for  each  element  under  each  loading  condition.  Thus,  the 
ith  stress  constraint  is 

4'1  < 0 (3.4-1) 

where 


max  (<t> 

j.k 


ijk 


(b  » Zjj  , z j)  ) , 


J - 1 to  NMG;  and  k - 1 to  NLC 


(3.4-2) 


where  NMG  is  the  number  of  elements  in  the  ith  group  and  NLC  is  the 
number  of  loading  conditions.  Similarly,  a displacement  constraint 
for  the  ith  degree  of  freedom  is  imposed  on  the  worst  violation,  under 
all  loading  conditions 


4>1  - max  (b,tB,zI)},  j - 1 to  NLC 


(3.4-3) 
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In  this  procedure,  the  number  of  violated  constraints  is  reduced 
considerably.  Also,  design  sensitivity  calculations  of  redundant  con- 
straints are  not  performed,  which  improves  efficiency  of  the  algorithm. 

A number  of  other  computational  aspects  are  also  considered  in 
the  present  algorithm.  For  example,  all  constraints  of  the  problem 
are  normalized  with  respect  to  their  limit  values.  Thus  each  con- 
straint is  expressed  as 

<Jii  = R*  - 1.0  < 0 

it 

where  R^  is  a response  ratio, 
straint  is  then  given  by 

(1.0  - R*)  (3.4-5) 

This  procedure  facilitates  comparison  of  amount  of  the  violation  in 
various  constraints. 

Further,  the  e-active  constraint  concept  [1,47,50]  is  imple- 
mented quite  readily  in  the  algorithm.  If  a constraint  is  nearly 
tight,  then  it  is  included  in  the  violated  constraint  set.  This 
procedure  increases  the  number  of  constraint  violations  but  eventually 
helps  in  avoiding  oscillations  in  constraint  violations. 

3.5  Computational  Considerations  in  Design 
Sensitivity  Analysis 

3*J  3J  3J  9 j 

Computation  of  the  sensitivity  vectors  — , - — , - — , and  — 

db  dZj  dZp  35 

and  sensitivity  matrices  ^ » f£-  • and  fj"  are 

BIB 


(3.4-4) 

The  correction  to  any  violated  con- 
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quite  easy,  once  the  form  of  the  cost  function  and  the  constraint  func- 
tion have  been  determined. 

Now  consider  computation  of  the  matrix  A of  Eq.  (2.6-6).  This 

requires  generation  of  matrices  Asl,  As2,  Ad,  and  Ae,  which  further 

require  the  calculation  of  matrices  A®1,  A^1,  A®2,  C„,C,  etc.  The 

1 3 d 2 

sl 

computation  of  matrix  A^  from  Eq.  (2.5—32)  can  proceed  substrucure- 
wise,  since  for  each  substructure  is  completely  uncoupled  and  is 
available  from  static  analysis  in  the  decomposed  form.  Let  y(r)  be 

the  number  of  components  of  that  depend  upon  z^r\  Then 

L - disl  1 

Yi  * E Y(r)*  The  matrix  — f — is  partitioned  as 

1 r**l  9zI 


0 


(3.5-1) 


where  ■ 


3z 


(r) 


is  a [Y(r)»  m(r)]  matrix. 


accordingly  as 


The  matrix  A 


sl 

I 


is  partitioned 


,sKl) 

AI 

X8K2) 


,sl(L) 

AI 


(3.5-2) 
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In  computing  the  matrix  A from  Eq.  (2.5-33),  a matrix  -r* — Q 

B oZ  ^ 

Is  required.  This  matrix  is  also  computed  quite  readily,  since  the 


matrix  Q in  partitioned  form  is  available  from  static  analysis.  There- 


fore, Q 


T 3£ 


islJ 


3z, 


is  given  as 


-si 


3z, 


3z 


(1) 


.(L)T  3^Sl 


dz 


(L) 


(3.5-6) 


. . T -'Si 

where  each  matrix  q''  (r)  is  an  tn(r),  Y(r)l  matrix. 


‘T  3lsl 

^ < 

3Z1 

T 

Computation  of  the  matrix  Q in  Eq.  (2.5-7)  is  also  carried  out 
by  considering  one  substructure  at  a time.  From  Eqs.  (3.5-4)  and 
(2.3-11)  one  obtains 


T 

Q C„ 


Q(1)Tcf> 


q(l)Tc<l> 


(3.5-7) 


In  calculation  of  matrices  C1  and  C_  of  Eqs.  (2.5-4)  and  (2.5-5), 
3Sb  13Si 

respectively,  the  terms  tt—  and  rr—  are  quite  simple,  once  the 

a D o D 

dependence  of  externally  applied  loads  on  the  design  variable  is  known. 
Usually,  S and  S are  taken  to  be  independent  of  the  design  variables. 

D I 


3SB  3SI 

Therefore,  the  terms  and  drop  out  of  these  equations.  Other 

matrices,  such  as  fg  zj , fg  (K^  Zg) . and 

3t>  ^KII  ZI^’  are  a8sembbed  quite  conveniently,  elementwise.  In  this 

case,  there  are  a few  nonzero  elements  ar.d  only  these  are  considered  in 

calculations.  Thus,  knowing  the  matrices  and  C^t  matrix  C of 

Eq.  (2.5-7)  is  easily  computed.  Also,  it  is  noted  that  computations 
s X s2 

for  Xg  and  Ag  proceed  simultaneously.  Further,  computation  of 
matrices  (Ky)  and  (My)  of  Eqs.  (2.5-28)  and  (2.5-29)  are 
performed  elementwise. 

In  case  design  variable  linking  of  elements  across  substructure 
boundaries  does  occur,  the  matrix  C ^ of  Eq-  (2.5-5)  is  not  uncoupled, 
as  given  by  Eq.  (3.5-4).  In  this  case,  one  needs  to  define  a global 
numbering  system  for  the  design  variables.  It  is  easily  seen  that  all 
computations  still  proceed  substructurewlse. 

Calculation  of  the  matrix  is  similar  to  that  of  the  matrix  A. 

3.6  Lagrange  Multiplier  Calculations 

The  Lagrange  multiplier  vector  p that  is  required  for  computation 

of  the  optimum  design  change  6b  is  computed  using  Eq.  (2.6-14).  In 

order  to  enhance  computational  efficiency,  special  properties  of  the 

matrix  A of  Eq.  (2.6-6)  can  be  exploited.  It  may  be  noted  that 
8X  s2  6 

matrices  A , A , and  A do  not  possess  any  special  properties,  where- 
as A^  does  when  only  fixed  bounds  on  design  variables  are  considered. 
Thus,  the  matrix  A of  Eq.  (2.6-6)  can  be  rewritten  in  partitioned  form 


Similarly,  of  Eq.  (2.6-7)  and  the  coefficient  matrix  H of  Eq. 

* 

(2.6-13)  can  be  partitioned  as 


44 


Further,  when  only  bounds  on  design  variables  are  considered,  the 
matrix  H ^ is  diagonal.  This  property  of  the  matrix  H is  effectively 
utilized  in  computation  of  Lagrange  multipliers.  By  expanding  Eq. 
(2.6-12)  and  writing  it  for  in  partitioned  form,  one  obtains 


H 

H 

X 

1 

H12 

1 

•H  *— 1 

L_2 

« 

m 

5i 

_H21 

1 

CM 

CM 

X 

1 

|_^2j 

62 

and 


-AT  Vf1  AJ 


From  Eq.  (3.6-10),  one  obtains 


* 

H 


* 


D 


(3.6-10) 


(3.6-11) 


(3.6-12) 


(3.6-13) 


(3.6-14) 
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■ 


(3.6-15) 

(3.6-16) 

(3.6-17) 

and 

A - - H12  (3.6-18) 

Thus  p*  and  p*  are  determined  from  Eqs.  (3.6-14)  and  (3.6-15),  respec- 

2 

tively.  The  vector  p is  computed  using  a similar  approach.  In  the 
present  work,  a Gaussian  elimination  procedure  with  "total  pivoting" 
is  used  for  solving  Eq.  (3.6-14).  Further,  this  method  is  adapted  to 
handle  dependent  set  of  equations  in  Eq.  (2.6-6).  Other  techniques, 
such  as  the  Gram-Schmidt  orthogonalization  method,  can  also  be  used 
for  removing  dependent  equations  but  they  are  usually  computationally 
inefficient. 

Lastly,  it  should  be  pointed  out  that  in  multi-dimensional  struc- 
tural problems,  the  number  of  violated  design  variable  constraints  can 
be  quite  large.  Hence,  the  above  approach  provides  considerable  savings 
in  computing  time,  as  well  as  in  computer  core  storage  requirements. 
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3.7  Step  Size  Determination 

Determination  of  a proper  step  size  n is  quite  critical  for  rapid 
convergence  of  the  optimization  algorithm.  There  are  many  techniques 
presented  in  the  literature  [31]  for  determining  a step  to  be  taken  in 
the  constrained  steepest  descent  direction.  A simple  technique  pre- 
sented in  Ref.  31  and  used  successfully  by  many  researchers  [27,33,41] 
is  also  used  here.  Briefly,  it  consists  of  choosing  n to  yield  a few 
percent  reduction  in  cost  function,  when  A<j>  « 0.  This  implies  that  for 
a desired  reduction  ratio  ? in  the  cost  function,  the  step  size  n can 
be  expressed  as 


n 


r x j 


(3.7-1) 


It  may  be  pointed  out  here  that  choosing  q is  still  an  art.  In 
the  present  work,  step  size  is  specified  based  on  the  amount  of  con- 
straint violations.  It  is  chosen  to  be  quite  large  if  the  design 
point  lies  in  the  interior  of  the  constraint  set,  whereas  a moderate 
step  size  is  selected  if  there  are  violations.  Further,  a provision 
can  be  made  to  change  the  step  size,  if  necessary,  during  computations. 

3.8  Convergence  Criterion  and  Computational 
Checks 

In  order  to  check  the  calculations  for  6b,  computational  checks 
of  Theorem  4-16  [31],  are  implemented  at  each  iteration  of  the 
algorithm.  For  easy  reference,  these  checks  are: 


. MrS? 


■ 


■*.y  . ■ " 

.* ' t 
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W 6b1 


0 


(3.8-1) 


T 1 

A 6b  » 0 


(3.8-2) 


T 2 
A 6b 


A4> 


(3.8-3) 


and 


< 0 


(3.8-4) 


Finally,  since  the  cost  function  of  Eq.  (2.4-1),  state  equations 
(2.2-3)  and  (2.3-1),  and  constraints  of  Eqs.  (2.2-5),  (2.2-6),  and 
(2.4-2)  meet  the  hypotheses  of  Theorem  4-17  [31],  the  | | 6b^ | | ^ is 
monitored  as  a convergence  check  for  the  ODPS.  It  is  shown  in  Ref.  [31] 
that  as  the  algorithm  approaches  the  optimum,  (jAb1)^  -*•  0. 
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CHAPTER  4 

OPTIMAL  DESIGN  OF  TRUSS  STRUCTURES 

4.1  Introduction 

In  this  chapter,  the  method  developed  for  general  ODPS  Is  applied 
to  plane  and  space  trusses.  These  types  of  structures  are  encountered 
quite  frequently  in  practical  situations.  Most  common  among  these  are 
industrial  buildings,  transmission  towers,  and  spacecraft  structures 
such  as  dispenser  trusses  and  antennas.  In  all  these  cases,  it  is 
desirable  that  the  structure  should  simultaneously  meet  strength, 
deflection,  and  natural  frequency  requirements  and  be  of  minimum  weight. 

The  class  of  trusses  considered  herein  is  assumed  to  have  a 
specified  geometry  and  the  loads  are  applied  only  at  the  joints.  The 
objective  function  is  taken  as  the  total  weight  of  the  structure  and 
the  design  variable  for  each  member  is  taken  as  its  cross-sectional 
area.  Finally,  it  is  assumed  that  the  trusses  are  loaded  so  as  to 
deform  only  elastically.  The  displacement  method  of  structural 
analysis  is  used,  considering  nodal  displacements  of  the  truss  as  the 
basic  state  variables.  Equations  developed  for  structural  analysis 
in  Chapter  2 are  used  here.  The  stiffness  and  mass  matrices  used 
herein  for  a truss  element  are  well  documented  in  the  literature  [34]. 
For  a complete  derivation,  the  reader  is  referred  to  Ref.  34. 
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4.2  Reduction  of  the  Method  to  Trusses 
The  ODPS  for  trusses  can  be  defined  as  follows:  find  the  cross- 
sectional  area  of  each  member  of  the  truss,  such  that  its  weight  is 
minimized  and  state  equations,  and  constrains  on  stress,  buckling,  dis- 
placement, frequency  and  member  size  remain  satisfied. 

The  cost  function  of  Eq.  (2.4-1)  is  a linear  function  of  design 
variables  and  it  is  given  as 


J(b)  = 


NG(r)  NM(i) 

l l 


i=l  j=l 


piVi 


(4.2-]) 


where 


NG(r) 

NM(i) 


= material  density  of  members  of  the  ith  group 
= cross-sectional  area  of  members  of  the  ith  group 
= length  of  the  jth  member  of  the  ith  group 
= number  of  groups  in  the  rth  substructure 
* number  of  members  in  ith  group. 


Since  the  cost  function  depends  only  on  the  design  variables,  some  of 

the  calculations  of  Chapter  3 are  eliminated  and  the  vector  of  Eq. 
3JT 

(2.5-24)  is  simply  tt—  , which  is  quite  easy  to  calculate. 

o D 

Next,  the  stress  or  buckling  constraint  of  either  Eq.  (2.4-3)  or 
(2.4-4),  for  a typical  member  is  written  as: 


- 1.0  < 0 


(4.2-2) 


f 
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C 3 

where  o and  c are  the  calculated  and  the  allowable  stresses  for  the 

member,  respectively.  In  order  to  simultaneously  implement  stress  and 

buckling  constraints,  the  allowable  stress  a is  chosen  as  follows: 

(i)  for  members  in  tension,  oa  ■ oa+,  where  oa+  is  the  allowable 

tensile  stress  for  the  member,  and 

(ii)  for  members  in  compression,  a3  * min(o3  ,ab),  where  o3  >0 

and  ob  > 0 are  allowable  compressive  and  critical  buckling 

stresses  for  the  member,  respectively. 

Sit"  3 — 

The  stresses  a and  a are  specified  by  the  designer,  whereas 
a*5  depends  on  the  Euler  buckling  load  and  is  given  as 


v2  E I 

I2  b 


(4.2-3) 


where  E,  £,  b and  I are  modulus  of  elasticity,  equivalent  length  (for 
truss  members  l = i) , design  variable,  and  moment  of  inertia  of  the 
typical  member,  respectively.  In  this  research,  it  is  assumed  that  the 
members  in  a design  group  have  the  same  cross-sectional  geometry  and 
material  properties.  Thus,  the  moment  of  inertia  can  be  expressed  as 

I =*  ci  b2  (4.2-4) 

where  a is  a positive  constant  that  depends  only  on  the  cross-sectional 
geometry  of  the  member.  This  constant  is  specified  by  the  designer 
quite  readily  [31,41].  Thus,  Eq.  (4.2-3)  is  rewritten  as 


ob  * 0 b 


(4.2-5) 
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where 


TT2  E « 

I2 


If  the  constraint  of  Eq,  (4.2-2)  is  violated,  then 


(4.2-6) 


A<jS8  ■*  -a 


■"  oC 


- 1.0 


(4.2-7) 


In  the  example  problems  considered  herein,  the  value  of  the  constant 
'a'  is  taken  to  be  unity. 

The  displacement  constraint,  of  either  Eq.  (2.4-3)  or  (2.4-4),  for 
a typical  degree  of  freedom  is  expressed  as: 


= 


- 1.0  < 0 


(4.2-8) 


where  z and  z are  the  calculated  and  the  allowable  displacements, 
respectively.  If  this  constraint  for  the  typical  displacement  com- 
ponent is  violated,  then 


A$8  - - 


- 1.0 


(4.2-9) 


It  may  be  pointed  out  here  that  static  analysis  and  constraint 
checks  on  stress,  buckling,  and  displacement  proceed  substructurewise. 
The  sensitivity  analysis  proceeds  as  explained  in  Chapter  3 and 
matrices  Afll  and  A82  are  assembled  at  this  stage. 


J 
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The  constraint  of  Eq.  (2.2-6)  is  Imposed  only  on  the  lowest 

2 

eigenvalue  of  the  truss  U - (2irf)  ).  Using  the  method  stated  in 
Section  2.3.2,  the  lowest  eigenvalue  x,  and  the  associated  eigenvector 
y are  obtained.  Thus,  the  eigenvalue  constraint  is  written  as 


$e(c)  = 1.0  - < 0 


(4.2-10) 


where  i,  is  related  to  a resonant  frequency  of  the  structure.  If  this 
constraint  is  violated,  then 


A*e  - - (1.0  - |~) 

‘’o 


and 


_ JL 

* Co 


(4.2-11) 


Finally,  the  design  variable  constraint  4>d(b)  of  Eq.  (2.2-5)  is 
considered.  For  a typical  design  variable  it  is  expressed  as 

bi  1 bi  1 b±  (4.2-12) 

where  bb  and  b^  are  the  lower  and  upper  bounds  on  the  ith  design 
variable,  respectively.  The  inequality  of  Eq.  (4.2-12)  may  be  split 
into  two  parts  as  follows: 

(i)  Lower  bound  design  variable  constraint 

d bi 

d*  (b)  5 1.0  - ~Y  <0  (4.2-13) 

bi  " 

and 


(ii)  Upper  bound  design  variable  constraint 
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*d(b)  = -J  - 1.0  < 0 

bi 

If  a constraint  of  Eq.  (4.2-13)  is  violated,  then 


and 


(4.2-14) 


(4.2-15) 


(4.2-16) 


The  upper  bound  design  variable  constraint  is  treated  in  a similar 
way. 


4.3  Computer  Program 

In  order  to  test  the  efficiency  of  the  formulation  of  ODPS,  a 
computer  program  that  implements  the  algorithm  of  Chapter  2 is 
developed  for  plane  and  space  trusses.  A general  flow  diagram  for 
the  program  is  given  in  Fig.  4.1. 

All  the  computational  aspects  of  Chapter  3,  i.e.  multiple  loading 
conditions,  design  variable  linking,  worst  constraint  violation  concept, 
e-active  constraint  concept,  normalization  of  constraints,  and  compu- 
tational checks  have  been  incorporated  in  this  program.  Further,  the 
element  stiffness  and  mass  matrices  for  truss  elements  are  generate*4 
in  each  design  cycle.  In  frequency  analysis,  provision  is  made  in  the 
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If  any  consCraint  is  not  imposed,  then  skip  the  corresponding  portion 
of  the  program. 


Figure  4.1.  Flow  Diagram  for  the  Computer  Program 


program  for  a restart  of  the  method  for  minimization  of  Rayleigh's 
Quotient  after  a certain  number  of  iterations.  This  approach  speeds 
up  convergence  of  the  conjugate  gradient  method. 

In  order  to  improve  the  efficiency  of  the  method,  provision  for 
stress-ratio  design  is  made  in  the  computer  code.  In  this  concept, 
design  variables  (member  areas  in  the  present  case)  are  computed  from 
the  condition  that  stress  in  each  member  be  at  its  limiting  value.  Tt 
may  be  pointed  out  here  that  this  does  not  yield  an  optimum  design 
under  even  the  stress  constraints  for  indeterminate  structures,  but 
gives  a good  starting  point  for  the  optimal  design  algorithm.  A 
parameter  is  defined  in  the  computer  program  for  controlling  the  number 
of  stress-ratio  design  cycles. 

Also,  provision  is  made  in  the  computer  program  for  assigning  a 
predetermined  value  to  any  design  variable  at  the  start  of  the  program. 
Thus,  the  number  of  design  variables  can  be  less  than  the  number  of 
groups  for  the  truss.  This  is  a valuable  feature  in  the  program,  since 

it  allows  a designer  to  fix  some  member  areas  of  the  truss. 

1 2 

The  6b  and  6b  components  corresponding  to  some  design  variables 
attain  zero  values  after  a few  design  Iterations.  This  indicates  that 
these  design  variables  are  at  their  optimum  value  at  the  current  design 
iteration.  It  Is  assumed  that  they  will  remain  at  these  values  in 
subsequent  Iterations  and  are,  therefore,  kept  fixed.  However,  at  the 
final  solution  all  design  variables  are  released  and  the  convergence 
criterion  is  checked.  If  the  criterion  is  violated,  then  the  problem 
is  resolved.  Also,  stress  constraints  on  all  members  of  the  structure 
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are  Imposed,  irrespective  of  whether  the  design  variable  associated 
with  the  member  is  fixed  or  free. 


4.4  Example  Problems 

The  following  three  truss  examples  are  solved  using  the  computer 
program  developed  for  trusses  with  substructuring  formulation: 

(i)  Example  4.1:  10  Member  Cantilever  Truss 

(ii)  Example  4.2:  200  Member  Plane  Truss 

(iii)  Example  4.3:  63  Member  Wing-Carry-Through  Structure  (WCTS) 

These  examples  have  been  studied  by  several  researchers  [1-8]  and  thus 
serve  as  standard  problems  for  testing  a new  algorithm.  In  subsequent 
subsections  each  problem  is  discussed  separately.  The  results  obtained 
are  first  compared  with  a similar  formulation  developed  without  using 
the  sub structuring  concept  [1],  and  then  with  completely  different 
methods  [2-8].  The  computational  times  (in  double  precision)  reported 
herein  are  for  an  IBM  360/63(H)  computer. 

In  each  example,  the  constant  a used  to  compute  the  moment  of 
inertia  of  each  member,  for  imposing  buckling  constraints,  is  taken 
as  unity  (see  Eq.  (4.2-4))  and  the  value  Of  e is  taken  as  0.1%,  for 
imposing  e-active  constraints.  Further,  the  weighting  matrix  W (see 
Eq.  (2.6-8))  is  taken  as  an  identity  matrix. 

4.4.1  Example  4.1:  10  Member 

Cantilever  Truss 

Figure  4.2  shows  geometry  and  dimensions  of  the  10  member  canti- 
lever truss.  In  the  present  formulation,  the  truss  is  divided  into 
two  substructures,  by  partitioning  it  at  joints  3 and  4.  Thus,  the 


A 
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first  substructure  consists  of  members  2,  4,  6,  9,  and  10,  with  four 
boundary  and  four  interior  degrees  of  freedom.  The  second  substructure 
consists  of  members  1,  3,  5,  7,  and  8,  with  four  boundary  and  no 
interior  degrees  of  freedom.  The  truss  is  designed  for  the  two  loading 
cases  given  in  Table  4.1. 

For  Case  I,  the  optimum  design  is  obtained  for  five  separate 
problems,  designated  1(IA1),  1(IA2),  1(IA3),  1(IA4),  and  1(IB).  The 
first  four  problems  are  solved  by  imposing  only  the  direct  stress  con- 
straints. The  stress  limit  for  all  members,  except  member  9,  is  kept 
at  ±25  ksi.  The  stress  limit  for  member  9 is  varied  as  ±25,  ±30,  ±50, 
and  ±70  ksi  (see  Table  4.1).  The  fifth  Problem  1(IB)  is  solved  by 
imposing  both  stress  and  displacement  constraints. 

For  Case  II,  optimum  designs  are  obtained  for  three  problems, 
designated  l(IIA),  l(IIB),  and  l(IIC).  The  first  two  problems  are 
similar  to  Problems  1(IA1)  and  1(IB),  respectively,  whereas  the  last 
Problem  l(IIC)  considers  all  constraints,  i.e.,  stress,  buckling, 
displacement,  and  frequency  constraints.  Bounds  on  design  variables 
are  imposed  in  all  the  problems.  Design  data  for  all  the  problems  of 
this  example  are  given  in  Table  4.1.  No  design  variable  linking  is 
used  in  this  example. 

Table  4.2  gives  the  final  design,  number  of  active  constraints, 
maximum  constraint  violation,  and  | | 6b^ | at  the  optimum  design  for 
each  problem.  It  also  gives  the  maximum  value  of  | | 6b^ | \^,  and  a 
comparison  of  CPU  times  per  design  iteration  for  the  present  fomulation 
and  a similar  forumulation  without  substructuring  [1).  A reduction  of 
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TABLE  4.1.  DESIGN  DATA  FOR  10  MEMBER 
CANTILEVER  TRUSS 


Modulus  of  elasticity  - 10  ksi 

Material  density 
Initial  value  of  A 
Lower  limit  on  A 
Upper  limit  on  A 
Displacement  limit 
Resonant  frequency 
Stress  limit 

Number  of  loading  conditions  * 1 
Load  data: 


- 0.10  lb/ in. 

2 

*1.0  in.  .where  A is  cross-sectional  area 

* 0.10  in.2 

* None 

* ±2.0  in. 

* 22.0  Hz. 

■ ±25  ksi 


Loading 

Case 

Number 

Loading 

Condition 

No.  of 
Loaded 
Joints 

Joint 

No. 



Load  Component  (kips) 

— 

in  direction ; 

X1 

X2 

x3 

mm 

0.0 

-100.0 

j 

0.0 

I 

1 

2 

fll 

0.0 

-100.0 

i 

0.0 

50.0 

0.0 

2 

0.0 

-150.0 

0.0 

II 

1 

4 

3 

0.0 

50.0 

0.0 

4 

0.0 

-150.0 

0.0 

Starting  reduction  ratio,  r: 

Problem  No.  - 1(IA1)  1(IA2)  1(IA3)  1(IA4)  1(IB)  l(IIA)  l(IIB)  l(IIC) 
r - 0.018  0.091  0.090  0.090  0.081  0.001  0.081  0.040 


For  Problems  1(IA2),  1(IA3),  and  1(IA4)  the  stress  limit  for  member 


61 

approximately  1.5  to  23%  in  CPU  times  per  design  iteration  ts  achieved 
with  the  present  formulation. 

For  Problem  1(IA1),  the  set  of  active  constraints  at  the  optimum 
includes  stress  constraints  in  members  1,  3,  4,  and  7 to  9,  as  well,  as 
minimum  size  constraints  on  members  2,  5,  6,  and  10.  The  total  number 
of  active  constraints  is  equal  to  the  number  of  independent  design 
variables.  For  Problem  1(IA2),  the  set  of  active  stress  constraints  is 
the  same  as  for  Problem  1(IA1).  For  Problems  1(IA3)  and  1(TA4),  the 
optimal  solutions  are  identical.  This  is  due  to  the  fact  that  for 
these  two  problems,  the  stress  constraint  for  member  9 at  optimum  is 
not  critical  (stress  in  member  9 = -37.50  ksi).  Thus  increasing  the 
stress  limit  in  this  member  from  ±50  ksi  to  ±70  ksi  does  not  change  the 
optimum  solution.  The  set  of  active  constraints,  for  the  final  design 
given  in  Table  4.2,  includes  stress  constraints  for  members  1 to  4,  6 
to  8,  and  10,  as  well  as  the  minimum  size  constraint  on  member  5.  For 
Problem  l(IIA),  the  set  of  active  constraints  at  optimum  includes  stress 
constraints  for  members  1,  3,  4 and  6 to  9,  along  with  minimum  size 
constraints  on  members  2,  5,  and  10.  The  set  of  active  constraints  at 
the  optimum  for  Problem  1(IB)  includes  only  the  downward  vertical  de- 
flection constraint  at  nodes  1 and  2.  For  Problem  l(IIB),  active  con- 
straints are  the  downward  deflection  at  node  2,  stress  in  members 
5 and  6,  and  minimum  size  in  member  2. 

The  last  Problem  l(IIC)  includes  all  the  constraints.  As  shown 
in  Table  4.2,  six  constraints  are  active:  frequency  bound,  stress  in 
members  2,  5,  and  6,  vertical  downward  deflection  at  node  2,  and 
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minimum  & ize  in  member  10.  It  may  be  pointed  out  here  that  the  CPU 
time  for  each  design  iteration  for  this  problem  depends  upon  the  number 
of  conjugate  gradient  iterations  performed.  For  this  problem,  an 
average  of  six  conjugate  gradient  iterations  are  performed  for  the 
frequency  analysis. 

Tables  4.3  and  4.4  give  the  details  of  the  computational  times 

and  the  cost  function  history  for  all  the  problems,  respectively.  It 

1 2 

is  noted  here  that  components  of  6b  and  fib  associated  with  design 

variables  for  member  5 of  Problem  1(IA1);  members  3,  5,  6,  and  10, 

of  Problems  1(IB);  and  members  5 and  10  of  Problem  l(IIB)are  equal  to 

zero  after  the  1st  iteration.  Hence  in  the  subsequent  iterations, 

these  design  variables  are  kept  fixed.  In  the  last  iteration,  these 

design  variables  are  released  in  order  to  check  the  final  convergence 

(see  Section  4.3).  Also,  in  each  problem  the  step  size  parameter  was 

monitored  and  adjusted  at  appropriate  iterations,  to  obtain  the  final 

solution,  The  rate  of  convergence  to  the  optimal  solution  (number  of 

design  cycles)  is  dependent  on  the  step  size.  A very  small  step  size 

can  slow  the  rate  of  convergence,  whereas  a very  large  step  may  cause 

oscillations  In  the  algorithm.  For  the  Problem  1(IA4),  an  additional 

computer  run  was  made  by  keeping  the  step  size  constant  (r  = 0.09). 

2 

The  initial  design  was  1.0  in.  for  all  members  and  the  stress-ratio 
design  was  performed  only  once.  The  final  solution  was  obtained  in 
23  iterations,  with  the  cost  function  history  being  1582.03,  1631.02, 
1656.86,  1636.97,  1612.25,  1591.74,  1570.42,  1517.97,  1533.06,  1518.48, 
1446.24,  1458.16,  1516.78,  1505.48,  1517.48,  1515.83,  1513.56,  1511.27, 


► 
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TABLE  4.3.  COMPUTATIONAL  TIMES  IN  SECONDS  FOR 
10  MEMBER  CANTILEVER  TRUSS 


Problem 

Total  CPU  Time  For 

Total  CPU 

Problem  Setup 

Stress-Ratio  Design 

★ 

Optimal  Design 

i ime 

(2)+(3)+(4) 

1 

2 

3 

4 

5 

1(IA1) 

0.547 

0.593  (10)+ 

2.046  (9)* 

3.186 

1 (IA2) 

0.525 

0.833  (14)+ 

1.447  (5)* 

2.805 

1(IA3) 

0.570 

0.067  (1)+ 

4.486  (15)* 

5.123 

1(IB) 

0.600 

0.832  (14)+ 

3.529  (15)* 

4.961 

l(IIA) 

0.570 

0.667  (10)+ 

— 

1.237 

KIIB) 

0.570 

0.665  (10)+ 

4.052  (16)* 

5.287 

Includes  CPU  time  required  for  analysis,  sensitivity  analysis,  and 
optimization. 

^Number  of  stress-ratio  design  iterations  performed. 

^Number  of  design  iterations  performed. 


64 


TABLE  4.4.  COST  FUNCTION  HISTORY  FOR 
10  MEMBER  CANTILEVER  TRUSS 


DI 

Number 


Weight  in  lbs. 


1(IA1) 


1(IA2) 

1 (IA3) 

1 ( IB) 

l(IIA) 

1593.17 

1606.66 

1592.26 

— 

1664.53* 

1572.26 

1650.65 

2486.86 

1551.30* 

1651.29 

3735.76 

1542.07 

1629.26 

4734.86 

1545.13* 

1603.79 

5133.65* 

1583.33 

5123.27 

1561.59 

5115.11 

1508.96 

5106.08 

1526.74 

5099.30* 

1515.62 

5091.32 

1524  06* 

5086.62* 

1509.12* 

5080.94 

1502.82* 

5078.92* 

1499.68* 

5076.84 

1498.11* 

5076.76* 

! 

1497.60* 

5076.56* 

— J 

1(IIB) 


1594.52 

1599.83 

1605.87* 

1587.77 

1597.11 

1595.97 

1594.61 

1593.99 

1593.46 

1593. 18* 


1664. 

2189. 

3236. 

4262. 

4711. 

4773. 

4751. 

4725. 

4716. 

4725. 

4722. 

4715. 

4705. 

4692. 

4679. 

4674. 

4676. 


53 

59 
84 
16 

79 

03 
52 

* 

45 

38 

86 

* 

60 

* 

80 

* 

48 

* 

04 

* 

51 

* 

74 

13^ 


"^Underlined  iteration  is  not  counted  as  it  is  used  for  checking 
solution  only. 

* 

Means  change  in  step  size. 

^Optimal  design  iteration. 
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1508.99,  1506.76,  1504.44,  1502.16,  1499.88,  and  1497.60.  The  final 
solution  is  the  same  as  in  Problem  1(IA3),  but  the  number  of  iterations 
in  the  latter  case  is  greater.  This  shows  the  importance  of  a proper 
step  size  for  faster  convergence  of  the  algorithm. 

Finally  Table  4.5  gives  a comparison  of  final  weights,  number  of 
analyses,  average  CPU  times  per  design  iteration  and  total  times  with 
other  available  results.  For  all  the  problems,  the  cost  function  is 
the  same  as  obtained  in  Refs.  2 and  3,  but  there  is  some  reduction  in 
comparison  with  Ref.  5.  Further,  there  is  a significant  reduction  in 
computational  times.  With  the  present  formulation,  the  CPU  time  per 
design  iteration  is  reduced  by  factors  of  up  to  5.6  and  2.4,  in 
comparison  to  CPU  times  reported  in  Refs.  2 and  4,  respectively.  Note 
that  the  solutions  are  obtained  using  different  computers.  In  order 
to  convert  the  times  from  CDC  7600  to  IBM  360/91  one  needs  to  multiply 
the  computing  times  with  CDC  7600  by  5 [3].  Also  in  a private  communi- 
cation with  IBM,  it  was  found  that  the  IBM  360/91  is  about  8 times 
faster  than  IBM  360/65.  Hence,  in  order  to  compare  the  times  with  the 
present  work,  one  needs  to  multiply  the  times  of  Ref.  2 (IBM  360/91) 
and  Refs.  3,  and  4 (CDC  7600)  by  factors  of  8 and  40,  respectively. 

4.4.2  Example  4.2:  200  Member 

Plane  Truss 

The  second  structure  (Fig.  4.3)  optimized  by  the  present  technique 
Is  a moderate  size  plane  truss  with  200  members  and  150  degrees  of 
freedom.  For  the  purpose  of  obtaining  computational  times  with  dif- 
nuafcers  of  substructures,  the  following  three  cases  are  con- 
rod  for  this  oxsapls: 
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TABLE  4.5.  COMPARISON  OF  RESULTS  FOR  10  MEMBER 
CANTILEVER  TRUSS 


Problem 

KIAl) 

1 ( 1A2) 

1 ( TA3) 

I (IB) 

l(IIA) 

HUB) 

Present 

1593.18 

V/i 

a* 

v/> 

1497.60 

5076.56 

1665.53 

4676.13 

Schmit-Miura  1 2] 

1593.23 

1545.17 

1497.65 

5076.85 

1664.55 

4676.93 

Weigh 
In  lb« 

Rizzl  [3] 

1593.18 

1545.13 

1497.60 

5076.66 

1664.53 

4676.92 

Dobbs-Nclson  [5] 

1622.00 

1582.00 

1509.00 

5080.00 

— 

5059.70 

Present 

9 

5 

15 

15 

0 

16 

o 

Schmlt-Mlura  (2] 

15 

15 

15 

12 

10 

10 

o 

Rizzl  [3] 

15 

15 

15 

11 

12 

12 

z 

Dobbs-Nelson  [5] 

10 

14 

114 

14 

-- 

11 

| 

+- 

a 

Present^ 

0.2273 

0.2895 

0.2991 

0.2353 

0.0 

0.2532 

V 

2 S 

H co 

Schmlt-MJura*[ 2] 

0.1714 

0.1767 

0.1811 

0.1697 

0.1982 

0.1755 

C 

fc" 

Rlzzl**tA) 

0.0166 

0.0161 

0.0187 

0.0141 

0.0169 

0.0142 

i 

Dobbs-Nelson  [5] 

-- 

-- 

— 

" 

Present* 

3.1860 

2.8050 

5.1239 

4.9613 

1.2370 

5.2868 

s s 

u 

c 

Schmlt-Mlura*  [2] 

2.8e04 

2.9296 

3.0028 

2.3537 

2.2783 

2.0383 

* 

4j  4 

O E 

Rizzl** [A ] 

0.3160 

0.3080 

0.3490 

0.2220 

0.2710 

0.2410 

H 

Dobbs-Nelson  [5j 

— 

” 

— 

— 

— 

Design  Iteration. 

♦ IBM  360/&S(H). 

IBM  360/91.  For  comparison,  time  should  be  multiplied  by  8. 

** 

CDC  7600.  For  rompar luoti,  time  should  be  multiplied  by  40. 
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Case  I:  The  truss  is  divided  into  two  substructures  by  partitioning 

it  at  joints  34  to  42.  Members  1 to  93  are  in  substructure 
1 and  94  to  200  are  in  substructure  2. 

Case  II:  The  truss  is  divided  into  three  substructures  by  partition- 

ing it  at  joints  29  to  33  and  57  to  61.  Members  1 to  76, 

77  to  152,  and  153  to  200  are  in  substructures  1,  2,  and  3, 
respectively. 

Case  III:  The  truss  is  divided  into  five  substructures  by  partition- 

ing it  at  joints  15  to  19,  29  to  33,  43  to  47,  and  57  to 
61.  Members  1 to  38,  39  to  76,  77  to  114,  115  to  152,  and 

153  to  200  are  in  substructure  1,  2,  3,  4,  and  5,  respec- 

tively. 

In  Cases  II  and  III,  nodes  1 to  5 are  also  considered  as  boundary 

nodes.  In  all  these  cases,  the  nodes  are  renumbered  in  order  to  pro- 
vide a minimum  upper  bandwidth  for  the  matrices  K,  , and 

In  Ref.  1 members  18,  25,  56,  63,  94,  101,  132,  and  139  are  linked  to 
the  same  design  variable  (no.  10),  whereas  in  the  present  formulation 
different  design  variables  are  associated  with  these  members.  The 
design  data  for  all  cases  of  this  example  are  given  in  Table  4.6. 

The  truss  is  optimized  in  each  case,  under  two  separate  constraint 
conditions:  Problem  2(A),  stress  constraints  only,  and  Problem  2(B), 
stress  and  displacement  constraints.  In  Case  III,  the  truss  is  opti- 
mized for  all  constraints  (Problem  2(C)),  i.e.,  stress,  displacement, 
buckling,  and  frequency  constraints.  Table  4.7  gives  the  final  solu- 
tion, number  of  active  constraints,  maximum  constraint  violation,  and 
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Loading  Condition  1.  One  kip  acting  in  positive  x^  direction  at  node 

points  1,  6,  15,  20,  29,  34,  43,  48,  57,  62,  71. 

Loading  Condition  2.  10  kips  acting  in  negative  direction  at  node 

points  1,  2,  3,  4,  5,  6,  8,  10,  12,  14,  15,  16, 
17,  18,  19,  20,  22,  24 71,  72,  73,  74,  75. 

Loading  Condition  3.  Loading  Conditions  1 and  3 acting  together. 

0.006  for  Problem  2(A) 

0.030  for  Problem  2(B) 

0.040  for  Problem  2(C) 


TABLE  4.7.  RESULTS  FOR  200  MEMBER  PLANE  TRUSS 


Group 

Member 

Optimum  Cross-Sectional  Area  in  in.^ 

No. 

Numbers 

Problem  2(A) 

Problem  2(B) 

Problem  2(C) 

1 

1,4 

0.1000 

0.1878 

2.4571 

2 

2,3 

0.1000 

0.1000 

2.6271 

3 

5,17 

0.2469 

4.7832 

4.6970 

4 

6,16 

0.1185 

0.1703 

1.1171 

5 

7,15 

0.1000 

0.1000 

0.32/0 

6 

8,14 

0.3675 

2.3462 

1.9650 

7 

9,13 

0.1000 

0.1876 

0.2549 

8 

10,12 

0.1000 

0.1000 

0.9190 

9 

11 

0.2454 

2.8809 

1.7910 

10 

18,25 

0.1000 

0.1000 

0.2210 

11 

19,20,23,24 

0.1000 

0.1000 

0.1000 

12 

21,22 

0.1000 

0.1000 

0.1000 

13 

26,38 

0.5815 

6.7767 

6.7390 

14 

27,37 

0.1000 

0.1000 

0.1170 

15 

28,36 

0.1736 

0.2361 

1.1600 

16 

29,35 

0.7008 

3.3133 

2.8450 

17 

30,34 

0.1000 

0.1732 

0.9100 

18 

31,33 

0.1000 

0.2227 

0.2447 

19 

32 

0.5792 

4.1473 

2.9270 

20 

39,42 

0.1035 

0.1000 

0.4541 

21 

40,41 

0.1000 

0.1000 

0.3257 

22 

43,55 

0.6760 

8.1292 

7.4171 

23 

44,54 

0.2489 

0.2476 

0.8670 

24 

45,53 

0.1000 

0.1000 

0.1572 

25 

46,52 

1.2622 

4.4206 

5.6820 

26 

47,51 

0.1000 

0.2802 

0.6880 

27 

48,50 

0.1000 

0.2673 

0.5190 

28 

49 

0.8402 

4.7929 

3.7401 

29 

57,58,61,62 

0.1000 

0.1000 

0.1408 

30 

59,60 

0.1000 

0.1002 

0.1190 

31 

64,76 

1.0080 

9.3889 

8.4891 

32 

65,75 

0.1000 

0.1000 

0.2398 

33 

66,74 

0.3103 

0.3362 

0.9540 

34 

67,73 

1.5955 

5.0733 

6.4100 

35 

68,72 

0.1000 

0.3008 

0.5181 

36 

69,71 

0.1000 

0.3096 

0.7290 

37 

70 

1.1731 

5.5744 

4.4271 

38 

77,80 

0.1892 

0.4967 

0.9281 

39 

78,79 

0.1000 

0.3865 

0.2628 

40 

81,93 

1.0402 

9.5196 

8.3681 

41 

82,92 

0.4108 

0.9366 

1.7621 

42 

83,91 

0.1000 

0.1000 

0.1490 

43 

84,90 

2.2576 

6.2617 

7.8100 

TABLE  4.7.  (cont'd.) 


Group 

Member 

Optimum  Cross-Sectional  Area  in  in.^ 

No. 

Numbers 

Problem  2(A) 

Problem  2(B) 

Problem  2(C) 

44 

85,89 

0.1000 

0.3508 

0.8980 

45 

86,88 

0.1000 

0.4835 

0.5779 

46 

87 

1.4567 

5.8679 

5.0921 

47 

95,96,99,100 

0.1000 

0.1000 

0.1012 

48 

97,98 

0.1000 

0.1000 

0.1241 

49 

102,114 

1.3481 

10.4800 

9.1651 

50 

103,113 

0.1111 

0.1108 

0.2005 

51 

104,112 

0.4795 

1.0313 

1.8621 

52 

105,111 

2.5909 

6.8203 

8.3780 

53 

106,110 

0.1007 

0.5012 

0.5439 

54 

107,109 

0.1000 

0.3754 

0.8989 

55 

108 

1.7922 

6.4768 

5.6151 

56 

115,118 

0.2932 

1.9807 

1.4791 

57 

116,117 

0.1000 

1.4784 

0.2829 

58 

119,131 

1.2625 

9.1546 

8.3571 

59 

120,130 

0.6050 

3.1979 

2.5540 

60 

121,129 

0.1000 

0.1000 

0.1477 

61 

122,128 

3.3587 

9.0271 

10.0360 

62 

123,127 

0.1000 

0.2074 

0.8499 

63 

124,126 

0.1069 

0.9717 

0.6821 

64 

125 

2.0771 

6.5338 

6.1791 

65 

133,134,137,138 

0.1000 

0.1000 

0.1000 

66 

135,136 

0.1002 

0.1219 

0.1210 

67 

140,152 

1.5748 

9.9624 

9.0011 

68 

141,151 

0.1325 

0.1341 

0.1983 

69 

142,150 

0.6817 

3.3000 

2.6930 

70 

143,149 

3.6920 

9.5771 

10.5110 

71 

144,148 

0.1092 

0.9814 

0.6610 

72 

145,147 

0.1000 

0.2269 

0.8559 

73 

146 

2.4116 

7.0561 

6.6231 

74 

153,156 

0.4250 

2.5500 

2.2599 

75 

154,155 

0.1000 

0.6074 

0.2332 

76 

157,169 

1.3251 

7.5376 

7.1521 

77 

158,168 

0.8462 

4.1216 

3.8200 

78 

159,167 

0.1000 

0.1000 

0.1020 

79 

160,166 

4.5880 

13.3290 

12.5160 

80 

161,165 

0.1000 

1.8691 

1.8370 

81 

162,164 

0.1276 

0.3045 

0.8090 

82 

163 

2.6772 

7.4246 

7.1681 

83 

171,172,175,176 

0.1000 

0.1000 

0.6620 

84 

173,174 

0.1002 

0.1000 

0.6611 

85  ' 

178,190 

1.6597 

8.2183 

7.1961 

7»—  'at 


Group 

No. 

Member 

Numbers 

2 

Optimum  Cross-Sectional  Area  in  in. 

Problem  2(A) 

Problem  2(B) 

Problem  2(C)1 

86 

179,189 

0.1156 

0.1000 

0.3654 

87 

180,188 

0.9181 

4.1916 

4.9040 

88 

181,187 

4.9213 

13.8330 

12.5350 

89 

182,186 

0.1339 

0.3354 

1.1069 

90 

183,185 

0.1002 

1.9082 

2.1110 

91 

184 

3.0100 

7.8840 

7.1751 

92 

191,194 

1.2392 

5.8649 

5.3130 

93 

192,193 

0.8521 

3.4248 

3.1651 

94 

195,200 

2.3257 

10.6560 

9.9201 

95 

196,199 

5.9232 

17.7770 

18.8650 

96 

197,198 

2.5718 

7.7140 

8.9630 

97 

98 

99 
100 

56,63 

94,101 

132,139 

170,177 

0. 1000 

0. 1000 
0.1000 
0.1000 

- -------  -r  . 

0.1000 

0.1000 

0.1000 

0.1000 

0.2210 

0.2210 

0.2210 

0.1000 

1 

♦j 

a 

o 

u 

< 

Wt.  in  lbs. 

■ H 

7,488 

28,963 

29,707 

No.  of  Active 
Constraints 

90 

28 

25 

Max.  Constraint 
Violation 

0.71xl0~4 

0. 97x10  3 

0. 14*10~2 

54.95 

53.01 

264.30 

Max.  I |6bX | |2 

273.60 

979.90 

— 

728.90 

1 
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| | 6b  j | 2 at  the  optimum  design  for  each  problem,  and  the  maximum  value 
of  | | 6b^|| 2*  Also,  final  solutions  could  be  obtained  in  about  12-15 
iterations  by  monitoring  the  step  size.  Table  4.8  gives  active  con- 
straints for  the  optimum  solution  reported  in  Table  4.7. 

Table  4.9  shows  a comparison  of  computational  times  between  the 
present  formulation  and  a similar  formulation  without  substructuring 
[1].  The  CPU  time  per  design  iteration  depends  upon  many  factors, 
such  as  the  total  number  of  constraint  violations,  the  number  of  load- 
ing conditions,  the  number  of  substructures  under  which  violations 
occur,  etc.  The  CPU  times  per  design  iteration  given  in  Table  4.9  are 
for  the  number  of  active  constraints  reported  in  Table  4.7.  For  the 
Problem  2(A),  violations  occurred  in  all  substructures,  in  all  loading 
conditions,  and  for  all  three  cases.  For  the  stress  and  displacement 
constraint  case  i.e..  Problem  2(B),  violations  occurred  in  all  sub- 
structures and  in  all  loading  conditions  in  some  initial  iterations. 
Near  the  optimum  point,  however,  only  displacement  violations  occurred 
under  loading  conditions  2 and  3,  and  in  substructure  1 for  all  cases. 

A comparison  of  CPU  times  for  structural  analysis  indicates  that 
as  the  number  of  substructures  is  increased  from  2 to  5,  the  CPU  time 

for  analysis  is  reduced  by  a factor  of  more  than  2.  This  reduction  in 

> 

CPU  time  for  analysis  is  due  to  the  reduced  bandwidth  of  the  matrix 
(r) 

. In  comparison,  the  CPU  time  for  structural  analysis  is  50% 
greater,  with  two  substructures,  and  is  33%  lower  with  five  substruc- 
tures. The  CPU  time  per  design  iteration  is  reduced  in  all  cases, 
with  the  greatest  reduction  occurring  in  Case  III.  This  indicates  that 


active  constraints  at  the  optimum 

DESIGN  FOR  200  MEMBER  PLANE  TRUSS 


Problem  2(A):  Case  III 
(i) 


Active  Stress  Constraints  (50)* 

Members  U.  17,  26,  32,  39,  43,  49,  64,  70,  77,  87,  108,  115, 
125,  146,  163,  and  184  under  loading  condition  2; 

Members  14,  16,  28,  35,  44,  52,  66,  73,  82,  90,  104,  110,  111, 
120,  126,  128,  142,  148,  149,  153,  158,  164,  166,  169,  179,  180, 
190,  193,  194,  198,  199,  and  200  under  loading  condition  3. 
Active  Minimum  Size  Constraints  (40)* 

D«lg„  Variables  1,  2,  5,  7,  8,  10  to  12,  14,  17,  18,  21,  24, 

26,  27.  29,  30,  32,  35,  36,  39,  42,  44,  45,  47,  48,  54,  57,  60, 
62,  65,  72,  75,  78,  80,  83,  and  97  to  100. 

Problem  2(B);  Case  I 


(ii) 


(i) 


(ii) 


Active  Displacement  Constraints  (5)* 

At  nodes  1 and  6 In  positive  ^ direction  under  loading  condi- 
tion 1;  and  at  nodes  3 to  5 in  negative  *2  direction  under 
loading  condition  3. 

Active  Minimum  Size  Constraints  (23)* 

Design  Variable.  2,  5,  8,  10  to  12,  14,  20,  21,  24,  29,  32,  42, 
47,  48,  60,  65,  78,  83,  84,  86,  97,  and  98. 


Problem  2(0:  Case  III 

(i)  Active  Stress  Constraints  (18)* 

Members  37,  40,  53,  79,  91,  and  167  under  loading  condition  1; 
Member  179  under  loading  condition  2; 

Members  7.  9,  18,  50,  56,  72,  88,  94,  110,  129,  and  132  under 
loading  condition  3. 
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(ii) 


(iii) 


TABLE  4.8.  (cont’d.) 


Active  Displacement  Constraints  (3) 

At  node  3 under  loading  condition  2,  and  at  nodes  4 and  5 under 
loading  condition  3 in  negative  direction. 

* 

Active  Minimum  size  constraints  (4) 

Design  Variables  11,  12,  65,  and  100. 


* 


Number  of  active  constraints 
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there  is  an  optimum  number  of  substructures,  for  a minimum  computing 
time.  When  the  number  of  substructures  is  below  or  above  this  optimum 
number,  efficiency  of  the  substructuring  technique  is  reduced. 

It  may  be  observed  from  Table  A. 7 that  many  design  variables  are 
at  their  lower  bounds  at  the  optimum  design.  A side  computation  was 
made  by  keeping  these  variables  constant  for  Case  I,  with  only  stress 
constraints,  in  order  to  determine  the  reduction  in  computing  time. 

For  this  run,  computing  time  per  design  iteration  was  12.4  seconds,  as 
compared  with  20.57  seconds,  a reduction  of  approximately  40%.  In  many 
practical  situations,  it  is  known  a priori  that  certain  members  of  the 
structure  must  have  assigned  minimum  areas  or  maximum  areas.  The 
above  result  indicates  that  the  designer  should  take  full  advantage 
of  such  knowledge  to  achieve  maximum  efficiency. 

For  Problem  2(C),  the  number  of  conjugate  gradient  iterations  for 
minimizing  the  Rayleigh  Quotient  is  restricted  to  five.  In  the  initial 
few  iterations,  the  design  changes  significantly.  Therefore,  it  is  not 
advantageous  to  obtain  a precise  solution  of  the  eigenvalue  problem. 

The  eigenvector  obtained  from  the  previous  design  cycle  serves  as  a 
good  approximation  for  the  current  design  cycle.  It  may  also  be  noted 
that,  near  the  optimum,  the  conjugate  gradient  method  takes  only  one 
iteration  for  convergence.  For  this  problem  a reduction  of  about  40% 
in  computing  time  is  achieved,  as  compared  to  results  in  Ref.  1. 


The  results  of  this  example  are  not  compared  with  Refs.  51  and 
52,  since  the  comparison  was  already  made  in  Ref.  1. 
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4.4.3  Example  4.3:  63  Member  Wing-Carry-Through 

Structure  (WCTS) 

Figure  4.4  shows  the  truss  idealization  of  the  63  member  wing- 
carry-through  structure  (WCTS) . This  problem  was  originally  optimized 
by  Berke  and  Khot  [6]  and  later  by  Schmit  and  Kiura  [2],  and  Dobbs  and 
Nelson  [5].  For  details  regarding  idealization  of  this  structure,  the 
reader  is  referred  to  Ref.  6. 

In  the  present  formulation,  the  structure  is  divided  into  two 
substructures  by  partitioning  it  at  nodes  7 through  10.  The  first 
substructure  consists  of  31  members:  1 to  8,  17  to  20,  23  to  27, 

30  to  33,  42  to  49,  58  and  59.  The  remaining  32  members  are  in  the 
second  substructure.  Tables  4.10  and  4.11  define  the  nodal  coordinates 
and  member  locations,  respectively.  Design  data  for  the  structure  is 
given  in  Table  4.12.  The  loading  conditions  given  in  Table  4.12  are, 
according  to  Ref.  6,  "derived  from  a 2g  condition  of  a half  million 
pound  aircraft  and  from  two  assumed  wing  positions".  The  structure  is 
designed  for  two  cases: 

(i)  with  stress  constraints  only  (Problem  3(A)) 

(ii)  with  stress  and  torsional  stiffness  constraints (Problem  3(B)). 
The  torsional  stiffness  constraint  of  Ref.  6 is  imposed  by  limiting 
the  relative  displacement  of  nodes  1 and  2 in  the  x^  direction.  This 
constraint  is  incorporated  by  a minor  modification  in  the  present 
computer  program.  Finally,  no  design  variable  linking  is  used  in  this 
example. 

Table  4.13  gives  the  final  cross-sectional  areas,  value  of  the 
cost  function,  and  number  of  active  constraints,  maximum  constraint 
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Member 

End 

Nodes 

Member 

End 

Nodes 

Member 

End 

Nodes 

Member 

End 

Nodes 

1 

1,3 

17 

3,5 

33 

8,6 

49 

7,4 

2 

2,4 

18 

4,6 

34 

7,13 

50 

9,14 

3 

1,5 

19 

7,9 

35 

8,14 

51 

7,12 

4 

2,6 

20 

8,10 

36 

11,9 

52 

13,10 

5 

7,3 

21 

11,13 

37 

12,10 

53 

11,8 

6 

8,4 

22 

12,14 

38 

11,17 

54 

13,18 

7 

9,5 

23 

1,2 

39 

12,18 

55 

11,16 

8 

10,6 

24 

3,4 

40 

15,13 

56 

17,14 

9 

11,7 

25 

5,6 

41 

16,14 

57 

15,12 

10 

12,8 

26 

7,8 

42 

1,6 

58 

5,4 

11 

13,9 

27 

9,10 

43 

1,4 

59 

3,6 

12 

14,10 

28 

11,12 

44 

5,2 

60 

9,8 

13 

15,11 

29 

13,14 

45 

3,2 

61 

7,10 

14 

16,12 

30 

3,9 

46 

5,10 

62 

13,12 

15 

17,13 

31 

4,10 

47 

3,8 

63 

11,14 

16 

18,14 

32 

7,5 

48 

9,6 

TABLE  4.12.  DESIGN  DATA  FOR  63  MEMBER  WCTS 


3 

Modulus  of  elasticity  = 16  * 10  ksi 

(Titanium  alloy) 

Material  density  = 0.16  lb/in. 3 

Lower  limit  on  cross-sectional  areas  = 0.01  in.^ 

Upper  limit  on  cross-sectional  areas  - None 

Initial  value  of  cross-sectional  areas  = 20.0  in.^ 

Limit  on  relative  displacement  of 

nodes  1 and  2 in  x^-direction  =1.0  in. 

Number  of  loading  conditions  = 2 

Load  data: 


Loading 

No.  of  loaded 

Node 

Load  components  (kips)  in  direction 

Condition 

Nodes 

Number 

X1 

X2 

X3 

1 

2500.0 

-5000.0 

250.0 

1 

2 

2 

-2500.0 

5000.0 

250.0 

1 

5000.0 

-2500.0 

250.0 

2 

2 

2 

-5000.0 

2500.0 

250.0 

Starting  reduction  ratio,  r 


0.010  for  Problem  3(A) 
0.070  for  Problem  3(B) 


TABLE  4.13.  RESULTS  FOR  63  MEMBER  WCTS 


Member 

Number 


2 

Optimum  Cross-Sectional  Area  in  in. 

Problem  3(A) 

Problem  3(B) 

38.0780 

37.4460 

36.1020 

37.5610 

51.9230 

52.7160 

54.2430 

53.6230 

25.0690 

23.8050 

28.2220 

28.8300 

17.4610 

17.1540 

20.7660 

21.4020 

24.8100 

26.2760 

26.9370 

25.0950 

7.5854 

8.8686 

9.0049 

8.9586 

24.0470 

23.4150 

20.4880 

19.5970 

4.2966 

5.1817 

2.9538 

3.0207 

37.0600 

37.09*0 

37.1710 

37.3540 

0.0100 

0.0100 

0.0100 

0.0100 

0.0470 

0.0921 

0.0100 

0.0100 

0.1048 

0.0100 

1.1187 

0.0100 

0.0312 

0.0100 

Member 

2 

Optimum  Cross-Sectional  Area  in  in. 

Number 

Problem  3(A) 

Problem  3(B) 

26 

2.9000 

4.1985 

27 

1.3034 

0.9148 

28 

4.7738 

3.2278 

29 

0.8632 

0.0100 

30 

2.6809 

8.0127 

31 

2.7707 

8.0586 

32 

5.8034 

9.1997 

33 

5.6645 

9.0427 

34 

6.0951 

9.6634 

35 

5.6651 

9.0430 

36 

2.8156 

8.4160 

37 

2.7665 

8.0605 

38 

2.6920 

7.9197 

39 

2.7716 

8.0652 

40 

5.7531 

9.2645 

41 

5.6651 

9.0417 

42 

16.0290 

24.3540 

43 

18.3970 

26.0850 

44 

11.4900 

19.3100 

45 

13.8640 

21.0440 

46 

11.7350 

17.0110 

47 

6.1777 

12.8530 

48 

11.7930 

18.4540 

49 

13.9880 

19.8470 

50 

7.3508 

6.5775 

TABLE  4. 13.  (cont'd.) 


Member 

Number 


^ No.  of  Active 
^ Constraints 

4J  

a 

° Max.  Constraint 
£ Violation 

IlSb1! 


Max.  I I fib  II 


o 

Optimum  Cross-Sectional  Area  in  in. 

Problem  3(A) 

Problem  3(B) 

7.7067 

5.6535 

5.3699 

8.4259 

0.0285 

3.8917 

3.5992 

6.1081 

9.7775 

12.1060 

4.2407 

5.7593 

0.1179 

1.5657 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

4975.06 

6117.38 

55 

23 

-4 

-3 

0.52x10 

0.34x10 

3.68 

0.91 

21.64 

39.46 

V 
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violation,  and  | 1 6b^"  | | ^ at  the  optimum  for  both  the  problems  and  the 
maximum  values  of  | Jfib1 1 1 2 * For  Probl-ein  3(A),  the  set  of  active  con- 
straints at  the  optimum  includes  constraints  for  members  1 to  15,  17, 
18,  21,  23  to  25,  28,  50,  and  52,  under  loading  condition  I,  and  for 
members  26,  30,  32  to  36,  38  to  49,  and  54  to  56,  under  loading  condi- 
tion 2.  Minimum  size  constraints  on  members  19,  20,  22,  and  58  to  63 
are  also  active.  For  Problem  3(B),  the  set  of  active  constraints  at 
the  optimum  includes  stress  constraints  for  members  1 to  8,  10  to  18, 
21,  28,  29,  50,  and  51,  under  loading  conition  1,  and  the  torsional 
stiffness  constraint  under  loading  condition  2. 

Table  4.14  gives  the  cost  function  history  for  these  problems. 
Problem  3(A)  took  16  iterations  to  converge  and  Problem  3(B)  took  22 
iterations.  Details  of  monitoring  the  step  size  are  marked  in  the 
table.  A side  run  for  Problem  3(A)  was  made  with  a constant  step  size 

(r  * 0.09)  and  the  same  optimum  was  achieved  in  34  iterations.  Further 

1 2 

it  is  noted  here  that  for  Problem  3(B)  components  of  6b  and  6b 
associated  with  design  variables  for  members  19,  20,  22  to  25,  29, 
and  58  to  63  are  equal  to  zero  after  the  9th  iteration.  Hence,  in 
the  subsequent  iterations  for  this  problem,  these  design  variables 
are  kept  fixed. 

Table  4.15  gives  the  computational  times  for  these  problems. 
Finally,  Table  4.16  gives  a comparison  of  final  weights  achieved, 
number  of  analyses,  average  CPU  time  per  design  iteration,  and  total 
CPU  times  with  other  available  results.  For  both  the  problems,  the 
cost  function  is  practically  the  same  as  reported  in  Ref.  2,  whereas  in 
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TABLE  4.15.  COMPUTATIONAL  TIMES  IN  SECONDS  FOR 
63  MEMBER  WCTS 


Problem 

j.<  tal  CPU  Time  For 

Total  CPU 

Problem  Setup 

Stress-Ratio  Design 

Optimal  Design 

i lme 

(2)+(3)+(4) 

1 

2 

3 

4 

5 

3(A) 

1.987 

8.955  (10)+ 

97.242  (16)+ 

108.182 

3(B) 

1.988 

L— . ..  — ...  - 

11.641  (13)+ 

— 

62.221  (22)+ 

75.850 

* 

Includes  CPU  time  required  for  analysis, 
optimization. 

^Number  of  stress-ratio  design  iterations 
^Number  of  design  iterations  performed. 


sensitivity  analysis, 


performed. 


and 
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comparison  with  Ref.  6,  it  is  about  1%  lower.  Further,  in  comparison 
with  Ref.  5,  the  final  weight  achieved  by  the  present  method  is  about 
1%  lower  for  Problem  3(A)  and  it  is  9%  lower  for  Problem  3(B).  Finally, 
comparing  the  CPU  time  per  design  iteration  with  Ref.  2,  it  is  ob- 
served that  with  the  present  formulation  there  is  a reduction  by 
factors  of  8 and  19  for  Problems  3(A)  and  3(B),  respectively  (Table 
4.16). 

4.5  Discussion  and  Conclusions  Based  On 
Results  of  Chapter  4 

In  this  section,  conclusions  based  on  the  results  of  this  chapter 
are  presented.  Results  of  Examples  4.1  and  4.2  are  first  discussed 
and  compared  with  results  obtained  without  the  substructuring  formula- 
tion [1].  Then  all  solutions  are  discussed  and  compared  with  other 
solutions  available  in  the  literature  [2-8]. 

The  first  example,  a 10  member  cantilever  truss,  is  a small  scale 
structure.  The  CPU  time  per  design  iteration  for  problems  of  this 
example  is  reduced  by  up  to  23%,  as  compared  with  CPU  times  without 
substructuring  (Table  4.2).  Example  4.2  is  a moderate  size,  200  member 
plane  truss.  For  this  example,  the  reduction  in  computing  time  is 
more  significant,  6.5  to  66%  (Table  4.13). 

There  are  three  major  phases  of  any  optimal  design  algorithm: 

1)  structural  analysis,  2)  design  sensitivity  analysis,  and  3)  calcula- 
tion of  Iterative  design  changes.  Thus,  the  efficiency  of  any  method 
can  be  improved  by  increasing  the  efficiency  of  calculations  in  any  of 
the  three  phases.  The  present  examples  indicate  that  with  substructur- 
ing there  can  be  some  reduction  in  computing  time  of  the  structural 
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analysis  phase.  This  fact  is  also  observed  by  Kirsch,  Reiss,  and 
Shamir  [9],  and  Noor  and  Lowder  [10].  However,  in  the  present  formula- 
tion most  of  the  efficiency  is  achieved  in  the  design  sensitivity 

analysis  phase.  For  example,  if  all  the  constraint  violations  take 

sl 

place  in  one  substructure,  the  adjoint  matrix  X is  computed  for  only 

that  substructure.  For  Case  I of  the  200  member  plane  truss  example, 

this  means  that  in  design  sensitivity  analysis  one  has  to  solve  only 

a 66  * 66  system  of  equations,  rather  than  a 150  x 150  system  of 

equations,  each  having  a half  band-width  of  20  (including  the  diagonal). 

12  3 

Further,  an  operation  count  shows  that  only  (y  nm  + y nm)  operations 
are  required  to  decompose  a symmetric  positive  definite  banded  matrix 
(n,  m)  and  n(2m  + 1)  operations  are  required  for  the  solution  of  system 
of  n equations  [35].  The  parameter  m is  defined  here  as  the  half  band- 
width, excluding  the  diagonal.  In  the  operation  count,  an  operation 
is  defined  to  consist  of  one  multiplication  which  is  near]y  always 
followed  by  an  addition.  Thus,  it  can  be  easily  seen  that  the  number 
of  calculations  is  reduced  considerably,  even  with  just  two  substruc- 
tures, as  compared  to  a formulation  without  substructuring.  The  system 
of  equations  is  further  reduced  when  the  number  of  substructures  is 
increased  (Case  II  and  III).  It  is  noted  here  that  Case  III  of  this 
example  has  the  optimum  number  of  substructures,  which  is  seen  from 
the  operation  count  as  well  as  from  the  results  obtained  (Table  4.13). 

Based  on  the  above  discussion  of  results,  structural  analysis  and 
design  sensitivity  analysis  with  substructuring  is  more  efficient, 
because  it  involves  manipulation  of  matrices  having  smaller  dimensions. 
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This  contributes  to  a reduction  in  computation  time,  as  well  as 
computer  storage  requirements.  For  optimal  design  of  very  large 
structures  with  the  formulation  of  the  ODPS,  the  reduction  in  computing 
time  is  expected  to  be  even  greater  than  that  achieved  for  the  present 
examples. 

Results  obtained  by  the  present  method  for  Examples  4.1  and  4.3 
are  compared  with  the  results  of  Refs.  2-6.  In  all  cases,  the  optimum 
weights  obtained  with  the  present  method  are  the  same  as  obtained  in 
Refs.  2 and  3,  whereas  they  are  better  than  those  reported  in  Refs.  5 
and  6 (Tables  4.5  and  4.16).  Computing  times  with  the  present  method 
are  significantly  better  than  those  given  in  Refs.  2 and  3.  Computing 
times  in  Refs.  5 and  6 are  not  available  and  thus  cannot  be  compared. 

With  the  present  formulation,  there  is  a reduction  in  CPU  time 
per  design  iteration  by  factors  of  up  to  5.6,  8.3,  and  19  for  Example 
4.1;  Problem  3(A);  and  Problem  3(B),  respectively,  in  comparison  to 
CPU  times  reported  by  Schmit  and  Muira  [2],  In  comparison  with  Rizzl 
[4],  there  is  a reduction  by  factor  of  2.4  in  CPU  time  per  design 
iteration,  for  Example  4.1. 

It  should  be  noted  that  in  the  present  algorithm,  selection  of  a 
proper  step  size  requires  some  experience.  A very  small  step  size  will 
slow  down  the  rate  of  convergence,  whereas  a large  step  size  may  cause 
oscillation  in  the  algorithm.  In  the  present  example  problems,  the 
step  size  was  monitored  and  adjusted  to  obtain  final  solutions.  More 
work  is  needed  in  the  area  of  step  size  selection. 


94 


CHAPTER  5 

OPTIMAL  DESIGN  OF  IDEALIZED  WING  STRUCTURES 

5 . 1 Introduc  t ion 

In  the  previous  chapter  it  was  shown  that  the  present  formulation 
offers  a reliable  and  efficient  means  of  obtaining  minimum-weight 
designs  for  general  trusses  subjected  to  static  loads  with  stress, 
displacement,  frequency,  and  minimum  size  constraints.  In  this  chapter, 
the  algorithm  developed  in  Chapter  2 for  general  ODPS  is  applied  to 
idealized  wing  structures,  with  the  intention  of  demonstrating  the 
feasibility  of  developing  a large  scale  general  purpose  finite  element 
structural  synthesis  capability  that  is  practical  and  efficient. 

Results  are  compared  with  the  results  available  in  the  literature 

[2, 3, 4, 7, 8]. 

The  class  of  wing  structures  considered  herein  is  assumed  to  have 
a preassigned  geometric  configuration  and  fixed  structural  material. 
Further,  it  is  assumed  that  these  structures  can  be  idealized  using 
the  following  three  types  of  finite  elements:  (i)  truss  elements  of 
uniform  cross-sectional  area,  (ii)  constant  strain  triangular  (CST) 
membrane  elements  of  uniform  thickness,  and  (iii)  symmetric  shear  panels 
(SSP)  or  symmetric  pure  shear  panels  (SPSP) . It  may  be  pointed  out  here 
that  the  SSP  element  differs  from  the  usual  SPSP  element  in  that  it 
carries  bending  and  shear  loads.  The  element  stiffness  matrices  for 


' - - ' • 


•J? 


- t-  . 


95 


truss  and  CST  elements  are  well  documented  In  the  literature  and  will 
not  be  derived  here.  For  a complete  derivation  the  reader  is  referred 
to  Ref.  34.  However,  element  stiffness  matrices  for  the  SSP  and  the 
SPSP  elements  cannot  be  easily  found,  so  their  derivation  is  briefly 
summarized  in  Appendix  2. 

5.2  Reduction  of  the  Method  to  Idealized 
Wing  Structures 

The  ODPS  for  idealized  wing  structures  is  defined  as  follows: 
find  the  design  variable  associated  with  each  element  of  the  wing  struc- 
ture so  that  its  weight  is  minimized  and  the  state  equations,  and  the 
constraints  on  stress,  displacement,  and  member  size  are  satisfied. 

The  design  variables  for  the  structure  are  taken  as  cross- 
sectional  areas  of  the  truss  elements  and  thicknesses  of  the  CST  and 
SSP/SPSP  elements.  Design  variable  linking  is  used,  limited  of  course 
to  pre-specified  groups  of  finite  elements  that  are  of  the  same  type. 
Thus,  the  cost  function  of  Eq.  (2.4-1)  is 


J(b) 


TP(r) 

l 

k-1 


T 

i-1 


NM^i) 

J-l 


piYi 


(5.2-1) 


where 

TP(r)  ■ type  of  element  in  the  rth  substructure, 

NG(k)  “ number  of  groups  in  the  kth  type  of  element, 
NM(i)  ■ number  of  members  in  the  ith  group, 

■ material  density  of  elements  of  the  ith  group. 
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■ design  variable  for  elements  of  the  ith  group, 

■ length  of  the  jth  truss  element  or  surface  area 
of  the  jth  CST/SSP/SPSP  element,  of  the  ith  group. 


The  cost  function  is  an  explicit  linear  function  of  design 
variables,  so  its  design  derivatives  are  easily  obtained. 

In  the  present  work,  CST/SSP/SPSP  elements  are  required  to  satisfy 
a design  criterion  based  on  the  Von  Mises  equivalent  stress  (also 
known  as  Maximum  Energy  of  Distortion  Theory).  This  theory,  proposed 
by  M.  T.  Huber  in  1904  and  later  developed  by  R.  Von  Mises  (1913)  and 
H.  Hencky  (1925),  is  in  excellent  agreement  with  experiments  on  ductile 
materials  that  are  subjected  to  biaxial  or  triaxlal  state  of  stress,  as 
compared  to  other  theories  of  failure.  For  a complete  development  of 
the  Von  Mises  equivalent  stress  criterion,  the  reader  is  referred  to 
Ref.  53.  According  to  this  criterion,  the  equivalent  stress  (o  ) for 
an  element  in  a general  state  of  stress  is  given  by 


_2  {(°11  °22)  + (°22  ' °33)  + (°33  " °11) 


22  33 


'33  ir 


2 2 

+ 6(o12  + °23  + 


■4»] 


(5.2-2) 


where  o^  (i,j*  * 1,2,3)  are  stress  components  at  a point  of  interest 
(x^,  x^,  x^)  in  the  domain  Ji  of  the  element.  For  CST  or  SSP  (SPSP) 
elements,  Eq.  (5.2-2)  reduces  to 
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It  should  be  noted  that  the  SSP  elements  Include  a uniform  shear 
stress  and  a linear  normal  stress  distribution  (see  Appendix  II). 

Thus,  the  maximum  value  of  the  Von  Mlses  equivalent  stress  in  Is 
easily  calculated  (see  Appendix  II). 

Next,  the  stress  constraint  of  either  Eq.  (2.4-3)  or  (2.4-4),  for 
a typical  element,  is 


1.0  < 0 


(5.2-4) 


c 

where  a Is  the  direct  stress  for  truss  elements,  or  the  maximum  Von 
Mlses  equivalent  stress  calculated  from  Eq.  (5.2-3)  for  CST/SSP/SPSP 
elements,  and  o Is  an  allowable  stress. 

The  design  sensitivity  analysis  for  CST/SSP/SPSP  elements  Is  quite 
straight  forward  and  proceeds  as  explained  in  Chapters  3 and  4,  once 
the  stress-displacement  relationship  is  known.  As  an  easy  reference, 
the  stress-displacement  relationships  for  SSP  and  SPSP  elements  are 
given  in  Appendix  II,  whereas  for  CST  and  truss  elements  the  reader  is 
referred  to  Ref.  34. 

Remaining  constraints,  i.e.,  displacement  and  design  variable  con- 
straints are  the  same  as  In  Chapter  4 and  are  treated  In  the  same 
manner . 


5.3  Example  Problems 

In  order  to  test  the  structural  synthesis  capabilities  of  the 
formulation  of  ODPS,  a computer  program  based  on  the  algorithm  of 
Chapter  2 and  flow  diagram  of  Fig.  4.1  was  developed  and  the  following 
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three  examples  of  idealized  wing  structures  are  solved: 

(i)  Example  5.1:  18  Element  Wing  Box  Beam 

(ii)  Example  5.2:  33  Design  Variable  Rectangular  Wing 

(iii)  Example  5.3:  150  (130)  Element  Swept  Wing. 

These  examples  have  been  studied  by  several  researchers  [2, 3, 4, 7,8]. 

In  the  subsequent  subsections,  each  example  is  discussed  separately. 
The  computational  times  (in  double  precision)  reported,  herein,  are 
for  an  IBM  360/65(H)  computer.  The  value  of  e for  imposing  e-active 
constraints  varies  from  0.1%  to  5%  in  all  example  problems.  The 
weighting  matrix  in  these  examples  is  taken  to  be  diagonal,  as  in 
Chapter  4,  but  different  from  the  identity,  because  the  design  vari- 
ables are  of  different  dimensions.  In  all  examples,  an  element  of  the 
weighting  matrix  associated  with  the  ith  design  variable  is  taken  as 


W 


i 


3 J - 
3bi  Wi 


(5.3-1) 


where  w^  > 0 is  a weighting  multiplier,  to  be  chosen  by  the  designer. 
In  all  examples,  the  value  of  w^  for  truss  elements  varies  from  0.50 
to  200,  whereas  for  CST,  SSP,  and  SPSP  elements  it  is  taken  as  unity. 

5.3.1  Example  5.1:  18  Element  Wing 

Box  Beam 

Figure  5.1  shows  geometry  and  dimensions  of  the  18  element  wing 
box  beam.  The  structure  is  symmetric  with  respect  to  the  x^  - x^ 
plane.  Thus,  only  the  upper  half  of  the  box  beam  is  idealized.  Two 
cases  for  this  example  problem  are  considered.  In  Case  I,  5 truss,  5 
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CST,  and  8 SSP  elements  are  used.  For  Case  II,  SPSP  elements  are  used, 
instead  of  SSP  elements.  In  the  formulation  of  the  OOPS,  the  idealized 
wing  structure  is  divided  into  two  substructures  by  partitioning  it  at 
nodes  3 and  4.  The  first  substructure  consists  of  2 truss  elements 
(1  and  3),  2 CST  elements  (1  and  2),  and  3 SSP/SPSP  elements  (1,  3, 
and  6).  The  remaining  elements  are  in  substructure  2.  Both  substruc- 
tures have  6 boundary  degrees  of  freedom.  The  first  substructure  has 
no  interior  degrees  of  freedom,  whereas  the  second  substructure  has 
9.  Design  data  and  load  data  for  this  example  are  given  in  Table  5.1. 
Further,  linking  for  CST  elements  is  used  (see  Table  5.2). 

For  both  cases,  Table  5.2  gives  the  final  design,  value  of  cost 
function,  number  of  active  constraints,  maximum  constraint  violation, 

| | fib'*' | | 2 at  the  optimum,  and  the  maximum  value  of  | J fib^  I I 2 * The  set 
of  active  constraints  at  the  optimum  includes: 

Case  I:  deflection  at  node  5 in  the  x^-direction,  under  loading 

condition  1;  Von  Mises  equivalent  stress  for  CST  element  5 
and  SSP  elements  3 to  5 and  8 under  loading  condition  1 and 
for  SSP  elements  1 and  2 under  loading  condition  2;  and 
minimum  element  size  for  truss  elements  2 to  5 and  for  SSP 
elements  6 and  7. 

Case  II:  deflection  at  node  5 in  the  x^-direction  under  loading 

condition  2;  Von  Mises  equivalent  stress  for  CST  element  2 
and  SPSP  elements  1 and  2 under  loading  condition  2;  and 
minimum  element  size  for  truss  elements  2,  3,  and  5,  and  for 
SPSP  elements  6 and  7. 
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TABLE  5.1.  DESIGN  DATA  FOR  18  ELEMENT 
WING  BOX  BEAM 


A. 


Data  Common  to  All  Type  of  Elements: 


Modulus  of  elasticity 

Material  density 

Poisson's  ratio 

Displacement  limit  at  nodes 
3 to  7 In  x^-directlon 

Number  of  loading  conditions 

Load  data: 


- 10A  ksl 

- 0.10  lb/in.3 

- 0.30 

- ±2.0  In. 

- 2 


Loading 

No.  of  Loaded 

Node 

Load  Component  (kips) 

in  direction 

Condition 

Nodes 

Number 

X1 

x2 

X3 

1 

1 

7 

0.0 

0.0 

5.0 

2 

1 

5 

0.0 

0.0 

10.0 

Starting  reduction  ratio,  r ■ |0.027  for  Case  I 

0.081  for  Case  II 


B.  Data  for  Individual  Type  of  Elements: 


Starting  weighting  multiplier,  w - |1 00.0  for  Case  I 

(for  trues  elements)  200.0  for  Case  II 


Type  of 
Element 

Stress 
Limit  ksl 

Design  Variable 

Upper  Limit 

TRUSS 

±10.0 

0.10  in.2 

None 

0.98  in.2 

CST 

10.0 

0.02  in. 

1.0  in. 

0.196  in. 

SSP 

10.0 

0.02  in. 

1.0  in. 

0.196  in. 

SPSP 

10.0 

0.02  in. 

1.0  in. 

0.196  in. 

TABLE  5.2.  RESULTS  FOR  18  ELEMENT 
WING  BOX  BEAM 


Type  of 
Element 


(In.2) 


Group 

Number 


Element 

Number 


Optimum  Design 


Case  I I Case  I 


SSP  For 
Case  I 
And 

SPSP  For 
Case  II 
(In.) 


Weight  in  lbs. 

B 

Jj  No.  of  Active  Constraints 

4J  — 

o. 

o Max.  Constraint  Violation 

4J  — 

< I I ,l1 


Max.  <5b  2 


3.8682 

1.39290 

0.1000 

0.10000 

0.1000 

0.10000 

0.1000 

0.54522 

0.1000 

0.10000 

0.084383 

0.122730 

0.053070 

0.053830 

0.037844 

0.038444 

0.370500 

0.088141 

0.224470 

0.089840 

0.128540 

0.074624 

0.120750 

0.073084 

0.091634 

0.075286 

0.020000 

0.020000 

0.020000 

0.020000 

0.030954 

0.024886 

402.73 

357.09 

14 

9 

0. 53xl0-3 

0.48xl0~ 3 

17.13 

44.18 

1513.00 

1204.00 
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Table  5.3  gives  the  cost  function  history.  It  is  noted  here  that 
1 2 

components  of  6b  and  6b  associated  with  design  variables  for  elements 
2 to  5 (truss),  14  and  15  (SSP)  of  Case  I,  and  2,  3,  and  5 (truss),  14 
and  15  (SPSP)  of  Case  II  are  equal  to  zero  after  the  1st  iteration. 

Hence  in  the  subsequent  Iterations,  these  design  variables  are  kept 
fixed.  The  details  of  computational  times  are  given  in  Table  5.4. 

Finally  for  both  cases.  Table  5.5  gives  a comparison  of  results. 
Compared  to  Refs.  2 and  3,  the  cost  function  is  practically  the  same, 
whereas  for  the  Case  II  there  is  a reduction  of  about  9%,  as  compared 
to  Refs.  7 and  8.  With  the  present  formulation,  the  total  CPU  time  is 
less  by  factors  of  2.2  and  3.2  for  Case  I,  as  compared  to  Refs.  2 and 
4,  respectively.  For  the  Case  II,  it  is  less  by  a factor  of  3.5,  as 
compared  to  Ref.  2.  Finally,  a comparison  of  CPU  time  per  design 
iteration  indicates  that  there  are  reductions  by  factors  of  4.3  and 
3.4,  as  compared  to  Refs.  2 and  4,  for  Case  I,  respectively  and  by  a 
factor  of  4.6,  as  compared  to  Ref.  2,  for  Case  II. 

5.3.2  Example  5.2:  33-Design  Variable 
Rectangular  Wing 

Figure  5.2  shows  the  geometry  and  dimensions  for  a wing  structure, 
which  is  symmetric  with  respect  to  the  xj”x2  P^ane  and  corresponds  to 
the  wing  middle  surface.  Thus,  only  the  upper  half  of  the  rectangular 
wing  is  modeled,  using  12  truss  elements  to  represent  spar  caps,  12 
CST  elements  for  the  skin,  and  15  SSP  elements  for  the  vertical  webs. 

In  the  present  formulation,  the  idealized  wing  structure  is  di- 
vided into  three  substructures  by  partitioning  it  at  nodes  5-6,  and  9-10. 
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TABLE  5.4.  COMPUTATIONAL  TIMES  IN  SECONDS  FOR 
18  ELEMENT  WING  BOX  BEAM 


Problem 

Total  CPU  Time  for 

Total  CPU 
Time 

(2)+(3)+(4) 

Problem  Setup 

Stress-Ratio  Design 

Optimal  Design* 

1 

2 

3 

4 

5 

Case  I 

0.847 

1.505  (4)+ 

9.80  (15)* 

12.152 

Case  II 

1.015 

1.316  (4)+ 

5.593  (9)* 

7.924 

Includes  CPU  time  required  for  analysis, 
and  optimization. 

^Number  of  stress-ratio  design  iterations 
^Number  of  design  iterations  performed. 


sensitivity  analysis. 


performed. 
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Truss  elements  1 to  4,  CST  elements  1 to  4,  and  SSP  elements  1 to  4,  | 

and  13  are  in  substructure  1 and  truss  elements  5 to  8,  CST  elements 
5 to  8,  and  SSP  elements  5 to  8 and  14  in  substructure  2.  The  remain- 
ing elements  are  in  substructure  3.  All  substructures  have  six 
interior  degrees  of  freedom.  The  boundary  degrees  of  freedom  for  the 
first,  second,  and  third  substructures  are  6,  12,  and  12,  respectively. 

Design  and  load  data  are  given  in  Tables  5.6  and  5.7,  respectively.  As 
in  the  previous  problem,  design  variable  linking  is  used  only  for  CST 
elements  (see  Table  5.8). 

Table  5.8  gives  the  final  design,  value  of  cost  function,  number 
of  active  constraints,  maximum  constraint  violations,  | | 6b^ | j ^ at  the 
optimum,  and  the  maximum  value  of  [jfib^jl^-  It  should  be  noted  that 
Table  5.8  gives  two  solutions,  one  at  the  15th  and  the  other  at 
the  42nd  iteration.  Both  designs  are  usable,  since  all  constraints 
are  satisfied.  Also,  the  weights  at  the  two  iterations  are  practically 
the  same.  However,  the  final  designs  are  somewhat  different  and  the 
values  of  | | 6b'1'  | j ^ are  quite  different.  The  value  of  | | 6b^  | j ^ at  the 
42nd  iteration  is  much  smaller  than  at  the  15th  iteration.  Further, 
the  cost  function  is  very  near  its  optimum  value  after  15  iterations, 
but  it  takes  27  more  iterations  to  significantly  reduce  the  value  of 
j 1 6b^  ! I 2 * This  behavior  is  observed  in  many  structural  problems,  in 
which  the  cost  function  is  linear. 

The  set  of  active  constraints  at  the  optimum  includes:  deflection 
constraint  at  node  14  in  the  x^-direction  under  loading  condition  2; 

t 

I 

I 

J 

2 * V *•  x . — : " - 1 * . .•  j *'■*  * ”■  ' ■ 


I 
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TABLE  5.6.  DESIGN  DATA  FOR  33  DESIGN  VARIABLE 
RECTANGULAR  WING 


Data  Common  to  All  Type  of  Elements: 


Modulus  of  elasticity  * 10.5  * 10^  ksi 

Material  density  - 0.1  lb/in. ^ 

Poisson's  ratio  =0.3 

Displacement  limit  at  nodes 

13  and  14  in  xj-direction  = ±11.0  in. 

Number  of  loading  conditions  » 2 (see  Table  5.7) 

Starting  reduction  ratio,  r *■  0.021 

B.  Data  for  Individual  Type  of  Elements: 

Starting  weighting  multiplier,  w±  «=  0.667  (for  truss  elements) 


Type  of 
Element 

TRUSS 

CST 


Stress 
Limit  in 
ksi 

±25.0 

25.0 


Design  Variable 

Lower  Limit 

Upper  Limit 

0.10  in.2 

1.50  in.2 

0.01  in. 

None 

0.01  in. 

None 

For  elements  1 to  12 


For  elements  13  to  15 


Initial  Design 


2.0  in.' 

j 0.08  in. 
0.04  in. 

0.04  in. 


r . 


4- 


tr;.-  ' r 


V 


! 

» 
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TABLK  5.7.  LOAD  DATA  FOR  33  DESIGN  VARIABLE 
RECTANGULAR  WINC 


j For  all  nodus,  x^  and  components  of  load  * 0.0 


— 

Node 

Number 

Component  of  Load  in 

xj  direction  (kips) 

Load  Condition  1 

Load  Condition  2 

3 

0.1658 

0.2996 

4 

0.4144 

0.3745 

5 

0.1426 

0.2576 

6 

0.3564 

0.3221 

! 7 

0.1311 

0.2367 

8 

0.3274 

0.2958 

9 

0.1177 

0.2127 

10 

0.2943 

0.2659 

11 

0.0978 

0.1768 

12 

0.2445 

0.2209 

13 

0.0381 

0.0689 

“ 

0.0953 

0.0861 

Ill 


TABLE  5.8.  RESULTS  FOR  33  DESIGN  VARIABLE 
RECTANGULAR  WING 


Type  of 

Group 

Element 

Optimum  Design 

Element 

Number 

Numbers 

At  15th  Dll' 

At  42nd  DI' 

1 

1 

1.32640 

1.3063 

2 

2 

1.37700 

1.3713 

3 

3 

0.83578 

0.82372 

4 

4 

0.91279 

0.92373 

Truss 

5 

5 

0.49935 

0.48875 

fin 

6 

6 

0.52112 

0.52810 

V>n.  ; 

7 

7 

0.24142 

0.22519 

8 

8 

0.24934 

0.26079 

9 

9 

0.10129 

0.10 

10 

10 

0.11074 

0.12404 

11 

11 

0.10 

0.10 

12 

12 

0.10 

0.10 

13 

1,2 

0.123930 

0.127470 

14 

3,4 

0.099328 

0.098920 

CST 

15 

5,6 

0.073211 

0.072230 

(in. ) 

16 

7,8 

0.047007 

0.046674 

17 

9,10 

0.020752 

0.020828 

18 

11,12 

0.010 

0.010 

19 

1 

0.059261 

0.059019 

20 

2 

0.080678 

0.072716 

21 

3 

0.041459 

0.041699 

22 

4 

0.064171 

0.063083 

23 

5 

0.031355 

0.029265 

24 

6 

0.047776 

0.050518 

SSP 

25 

7 

0.026171 

0.022464 

(in.) 

26 

8 

0.038431 

0.040775 

27 

9 

0.018753 

0.014846 

28 

10 

0.028262 

0.031025 

29 

11 

0.010 

0.010 

30 

12 

0.025174 

0.013188 

31 

13 

0.010 

0.010 

32 

14 

0.010 

0.010 

33 

15 

0.010 

0.010 

- 

Weight  in  lbs. 

100.344 

100.212 

No.  of  Active  Constraints 

10 

15 

■U 

a 

o 

Max. 

Constraint 

Violation 

0.88x10"^ 

0.72xl0-7 

4J 

< 

115b1 

Hi 

17.48 

1.94 

Max.  116b1 

ll2 

102.70 

102.70 

t 
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Von  Mlses  equivalent  stress  in  SSP  elements  1 and  3,  under  loading 
condition  2;  and  minimum  element  size  for  truss  elements  11  and  12, 

CST  elements  11  and  12  (design  variable  16),  and  SSP  elements  11,  13, 
14,  and  15.  Table  5.9  gives  a comparison  of  final  design  weight, 
number  of  design  iterations,  CPU  time  per  design  iteration,  and  total 
computing  time  which  includes  initial  set-up  time  (2:086  sec),  4 
iterations  of  stress-ratio  design  time  (3:528  sec),  and  15  iterations 
of  design  time  (17.689  sec).  The  cost  function  is  practically  the 
same  as  obtained  by  Rizzi  [3,4].  Comparing  the  CPU  time  per  design 
iteration  and  total  CPU  time  with  Ref.  4,  it  is  observed  that  with 
the  present  formulation  there  is  a reduction  by  factors  of  6.4  and  5.4, 
respectively. 

Lastly,  the  cost  function  history  (in  lbs.)  for  this  problem  is: 

64.476,  82.551,  99.790,  104.314,  104.157*,  102.071,  101.095,  100.797*, 

100.279,  100.185,  100.121*,  100.359,  100.360,  100.359,  100.365, 

100. 344  (*  means  change  in  step  size).  It  is  noted  that  for  this 

1 2 

problem  components  of  6b  and  6b  associated  with  design  variables  for 
elements  11  and  12  (truss),  and  31  to  33  (SSP)  are  equal  tD  zero  after 
the  4th  iteration.  Hence,  in  the  subsequent  iterations  for  the  prob- 
lem, these  design  variables  are  kept  fixed. 

5.3.3  Example  5.3:  150(130)  Element 

Swept  Wing 

Figure  5.3  shows  geometry,  dimensions,  and  element  numbering  of 
the  idealized  swept  wing  structure.  The  wing  is  symmetric  with  respect 
to  the  x^-x2  plane-  Thus,  only  'the  upper  half  of  the  swept  wing  is 


TABLE  5.9.  COMPARISON  OF  RESULTS  FOR  33  DESIGN 
VARIABLE  RECTANGULAR  WING 


Weight  in  lbs. 

No.  of  DI+ 

Av.  CPU  Time/DI^  in  sec 


Present 

100.344 

15 

1.179 


Rizzl  [3,4] 

100.542 

15 

0.189 


Total  CPU  Time  in  sec 


23.303 


3.142 


modeled,  using  20  truss  elements  to  represent  the  spar  caps,  60  CST 
elements  to  represent  the  skin,  and  70  SSP  elements  to  represent  the 
vertical  web9.  The  wing  without  spar  caps  Is  referred  to  as  Case  1, 
and  with  spar  caps  (truss  elements)  as  Case  II.  In  this  example, 
extensive  design  variables  linking  is  employed  (see  Table  5.16).  The 
number  of  independent  design  variables  describing  the  spar  cap,  skin, 
and  web  are  14,  11,  and  7,  respectively.  For  Case  I,  the  problem 
involves  a total  of  eighteen  design  variables  and  for  Case  II,  it 
increases  to  32. 

In  the  present  formulation,  the  idealized  swept  wing  structure 
is  divided  into  two  substructures  by  partitioning  at  nodes  17  to  20. 
The  first  substructure  consists  of  eight  truss  elements  (1-4,  11-14), 
twenty  four  CST  elements  (1-24),  and  twenty-five  SSP  elements  (1-4, 
11-14,  21-24,  31-34,  41-49).  The  remaining  elements  are  in  substruc- 
ture two.  Both  substructures  have  12  boundary  degrees  of  freedom. 

The  first  substructure  has  36  interior  degrees  of  freedom,  whereas 
the  second  substructure  has  72. 

Design  and  load  data  are  given  in  Table  5.10  and  5.11,  respec- 
tively. The  nodal  coordinates  defining  the  layout  of  the  swept  wing 
are  given  in  Table  5.12.  Tables  5.13,  5.14,  and  5.15  define  locations 
of  truss,  CST,  and  SSP  elements,  respectively. 

Table  5.16  gives  the  final  designs,  values  of  cost  function, 
number  of  active  constraints,  maximum  constraint  violations,  | 1 6b^ | | ^ 
at  the  optimum,  and  the  maximum  value  of  | | 6b^"  I 1 2 ' set  actiye 

constraints  at  optimum  includes: 
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TABLE  5.10.  DESIGN  DATA  FOR  150  (130)  ELEMENT 
SWEPT  WING 


A. 


Data  Common  to  All  Type  of  Elements: 


Modulus  of  elasticity 
(Aluminum  alloy) 

Material  density 

Poisson’s  ratio 

Displacement  limit  at  nodes 
41  and  44  in  xj-direction 

Number  of  loading  conditions 

Starting  reduction  ratio,  r 


10.6  x 10  ksi 

0.096  lb/in.3 
0.30 

±60.0  in. 

2 (see  Table  5.11) 

0.054  for  both  Case  I and  II 


B.  Data  for  Individual  Type  of  Elements: 

Starting  weighting  multiplier,  w^  ■ 20.0  for  Case  II  (for  truss 
elements) 


Type  of 

Stress 

Design  Variable 

Initial  Design 

Element 

Limit  in 
ksi 

Lower  Limit 

Upper  Limit 

With 

Spar  Caps 

Without 
Spar  Caps 

TRUSS 

±25.0 

0.01  in.2 

1.50  in.2 

0.02  in.2 

CST 

25.0 

0.02  in. 

None 

* * ** 
f 0. 20  in. 

‘ 4>4> 

0.10  in. 

0.30  in. 

SSP 

25.0 

0.02  in. 

None 

0.20  in. 



0.15  in. 

- d 

*Foe  elements  1 to  24. 

**For  elements  25  to  60. 


f 
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TABLE  5.11.  LOAD  DATA  FOR  150  (130)  ELEMENT 
SWEPT  WING 

I 

V 


For  all  nodes,  xj^  and  x£  components  of  load  ■ 0.0 


Node 

Number 



Load  in  X3 
direction 
kips 

— 

Node 

Number 

Load  in  x^ 
direction 
kips 

Node 

Number 

Load  in  X3 
Direction 
kips 

Load  Condition  1 

5 

1.282 

19 

0.206 

6 

2.581 

20 

0.431 

7 

3.398 

21 

35 

0.563 

8 

2.380 

22 

36 

0.383 

9 

0.978 

23 

37 

0.144 

10 

2.013 

24 

0.852 

38 

0.302 

11 

2.593 

25 

0.362 

39 

0.395 

12 

1.764 

26 

0.756 

40 

0.269 

13 

0.727 

27 

0.986 

41 

0.620 

14 

1.386 

28 

0.671 

42 

0.129 

15 

1.906 

29 

0.282 

43 

0.169 

16 

1.297 

30 

0.589 

44 

0.116 

17 

0.570 

31 

0.768 

18 

1.190 

32 

0.522 

Load  Condition  2 


5 

2.361 

19 

if 

0.402 

6 

3.876 

20 

sms 

0.646 

7 

2.308 

21 

AJL jJ 

35 

0.398 

8 

0.793 

22 

1.374 

36 

0.154 

9 

1.772 

23 

0.825 

37 

0.311 

10 

2.895 

24 

0.284 

38 

0.482 

11 

1.705 

25 

0.665 

39 

0.306 

12 

0.582 

26 

1.092 

40 

0.135 

13 

1.310 

27 

0.651 

41 

0.133 

14 

2.135 

28 

0.224 

42 

0.206 

15 

1.258 

29 

0.518 

43 

0.131 

16 

0.433 

30 

0.851 

44 

0.58C 

17 

1.047 

31 

0.508 

18 

1.719 

32 

0.175 

Node 

Number 


Coordinate  in  Inches 


dumber 

X1 

X2 

x3 

1 

0.0 

300.00 

10.00 

2 

0.0 

250.00 

15.00 

3 

0.0 

185.00 

13.00 

4 

0.0 

100.00 

5.00 

5 

100.0 

258.3332 

8.5833 

6 

100.0 

214.1667 

12.8333 

7 

100.0 

157.1667 

11.0833 

8 

100.0 

83.3333 

4.3333 

9 

190.0 

220.8333 

7.3083 

10 

190.0 

181.9167 

10.8833 

11 

190.0 

132.1167 

9.3583 

12 

190.0 

68.3333 

3.7333 

13 

260.0 

191.6667 

6.3167 

14 

260.0 

156.8333 

9.3667 

15 

260.0 

112.6333 

8.0167 

16 

260.0 

56.6667 

3.2667 

17 

325.0 

164.5833 

5.3958 

18 

325.0 

133.5416 

7.9583 

19 

325.0 

94.5417 

6.7708 

20 

325.0 

45.8333 

2.8333 

21 

385.0 

139.5833 

4.5458 

22 

385.0 

112.0416 

6.6583 

23 

385.0 

77.8417 

5.6208 

24 

385.0 

35.8333 

2.4333 

25 

440.0 

116.6666 

3.7667 

26 

440.0 

92.3333 

5.4667 

27 

440.0 

62.5333 

4.5667 

28 

440.0 

26.6667 

2.0667 

29 

490.0 

95.8333 

3.0583 

30 

490.0 

74.4167 

4.3833 

31 

490.0 

48.6166 

3.6083 

32 

490.0 

18.3333 

1.7333 

33 

535.0 

77.0833 

2.4208 

34 

535.0 

58.2917 

3.4083 

35 

535.0 

36.0917 

2.7458 

36 

535.0 

10.8333 

1.4333 

37 

570.0 

62.50 

1.9250 

38 

570.0 

45.75 

2.6500 

39 

570.0 

26.35 

2.0750 

40 

570.0 

5.00 

1.20 

41 

600.0 

50.00 

1.50 

42 

600.0 

35.00 

2.00 

43 

600.0 

18.00 

1.50 

44 

600.0 

0.00 

1.00 

2 


TABLE  5.13.  TRUSS  ELEMENT  LOCATIONS  FOR 
150  (130)  ELEMENT  SWEPT  WING 


Member 

End 

Nodes 

Member 

End 

Nodes 

1 

1,5 

11 

3,7 

2 

5,9 

12 

7,11 

3 

9,13 

13 

11,15 

4 

13,17 

14 

15,19 

5 

17,21 

15 

19,23 

6 

21,25 

16 

23,27 

7 

25,29 

17 

27,31 

8 

29,33 

18 

31,35 

9 

33,37 

19 

35,39 

10 

37,41 

20 

39,43 

mm 


TABLE  5.14.  CST  ELEMENT  LOCATIONS  FOR 
150  (130)  ELEMENT  SWEPT  WING 


Member 


Member 


Member 


1.2.5 

6.5.2 

2.3.6 

7.6.3 

3.4.7 

8.7.4 

5.6.9 
10,9,6 

6.7.10 

11.10.7 

7.8.11 

12.11.8 

9,10,13 

14.13.10 

10.11.14 

15.14.11 

11.12.15 

16.15.12 
13,14,17 
18,17,14 


14.15.18 

19.18.15 

15.16.19 

20.19.16 

17.18.21 
22,21,18 

18.19.22 

23.22.19 

19.20.23 

24.23.20 

21.22.25 

26.25.22 

22.23.26 

27.26.23 

23.24.27 

28.27.24 

25.26.29 

30.29.26 

26.27.30 

31.30.27 


27,28,31 

32,31,28 

29.30.33 

34.33.30 

30.31.34 

35.34.31 

31.32.35 

36.35.32 

33.34.37 

38.37.34 

34.35.38 

39.38.35 

35.36.39 

40.39.36 

37.38.41 

42.41.38 

38.39.42 

43.42.39 

39.40.43 

44.43.44 
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TABLE  5.15.  SSP  ELEMENT  LOCATIONS  FOR 
150  (130)  ELEMENT  SWEPT  WING 


Member 

End 

Nodes 

Member 

End 

Nodes 

Member 

End 

Nodes 

1 

1.5 

25 

19,23 

49 

15,16 

2 

5,9 

26 

23,27 

50 

17,18 

3 

9,13 

27 

27,31 

51 

18,19 

4 

13,17 

28 

31,35 

52 

19,20 

5 

17,21 

29 

35,39 

53 

21,22 

6 

21,25 

30 

39,43 

54 

22,23 

7 

25,29 

31 

4,8 

55 

23,24 

8 

29,33 

32 

8,12 

56 

25,26 

9 

33,37 

33 

12,16 

57 

26,27 

10 

37,41 

34 

16,20 

58 

27,28 

11 

2,6 

35 

20,24 

59 

29,30 

12 

6,10 

36 

24,28 

60 

30,31 

13 

10,14 

37 

28,32 

61 

31,32 

14 

14,18 

38 

32,36 

62 

33,34 

15 

18,22 

39 

36,40 

63 

34,35 

16 

22,26 

40 

40,44 

64 

35,36 

17 

26,30 

41 

5,6 

65 

37,38 

18 

30,34 

42 

6,7 

66 

38,39 

19 

34,38 

43 

7,8 

67 

39,40 

20 

38,42 

44 

9,10 

68 

41,42 

21 

3,7 

45 

10,11 

69 

42,43 

22 

7,11 

46 

11,12 

70 

43,44 

23 

11,15 

47 

13,14 

24 

15,19 

48 

14,15 

' 


— 


— 


a ' 
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TABLE  5.16.  RESULTS  FOR  150  (130)  ELEMENT 
SWEPT  WING 


Type  of 
Element 

Group 

Number 

Element 

Numbers 

Optimum  Design 

Case  I 

Case  II 

Truss 

(In.2) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1 

2 

3 

4 

5,6 

7,8 

9,10 

11 

12 

13 

14 

15,16 

17,18 

19,20 

0.030000 

0.010000 

0.010000 

0.010000 

0.010000 

0.010000 

0.010000 

1.500000 

0.637270 

0.339730 

0.716110 

1.002900 

0.389730 

0.010000 

15 

1-6 

0.203340 

0.191260 

16 

7-12 

0.177250 

0.171290 

CST 

17 

13-18 

0.156250 

0.151740 

(in 

• ) 

18 

19-24 

0.129290 

0.123110 

19 

25-36 

0.110810 

0.099201 

20 

37-48 

0.094851 

0.087910 

21 

49-60 

0.020000 

0.020000 

22 

1-4 

0.032631 

0.032648 

23 

5-10 

0.020000 

0.020000 

24 

11-14 

0.033091 

0.031825 

25 

15-20 

0.046713 

0.047229 

SSP 

26 

21-24 

0.221400 

0.228710 

(In.) 

27 

25-30 

0.097980 

0.115860 

28 

31-34 

0.088411 

0.084393 

29 

35-40 

0.060872 

0.055557 

30 

41-49 

0.035310 

0.035857 

31 

50-58 

0.061103 

0.068786 

32 

59-70 

0.105600 

0.099396 

§ 

Weight  In  lbs. 

2462.34 

2442.06 

4-1 

No.  of  Active  Constraints 

13 

22 

a. 

o 

Max. Constraint  Violations 

0.28  x 10“ ^ 

0.59  xlO-* 

u 

< 

116b1 

ll2 

33.12 

35.11 

Max 

• 1 l^b1! 1 2 

2647.00 

3263.00 

i 


Case  I : 


Von  Mises  equivalent  stress  /or  CST  element  8,  14,  20  under 


loading  condition  1,  for  SSP  elements  20,  21,  30,58,  and  61 
under  loading  condition  1,  and  3,  11,  and  42  under  loading 
condition  2;  and  minimum  element  size  for  CST  elements  49  to 
60  and  SSP  elements  5 to  10. 

Case  II:  For  the  final  design  with  spar  caps,  active  constraints  at 

the  optimum  are  the  same  as  for  Case  I,  with  the  addition  of 
minimum  area  constraints  for  truss  elements  1 to  10,  19  and 
20,  and  maximum  area  constraint  for  truss  element  11. 

The  set  of  active  constraints  is  found  to  be  the  same  as  in  Refs. 

2 and  3,  with  the  exception  of  stress  constraint  for  SSP  element  11 
under  loading  condition  2.  It  is  noted  that  none  of  the  tip  deflection 
constraints  are  active.  It  is  noted  that  for  Case  II  the  spar  cap,  web, 
and  skin  materials  account  for  3.57%,  12.61%,  and  83.82%  of  the  final 
wing  weight,  respectively.  For  Case  I,  the  web  and  skin  materials 
account  for  13.11%  and  86.89%  of  the  final  wing  weight,  respectively. 

In  Ref.  2,  spar  cap  material  accounts  for  less  than  0.5%  of  the  final 
wing  weight,  while  the  web  material  accounts  for  approximately  12.5% 
of  the  final  wing  weight  in  both  cases.  Thus,  the  distribution  of  the 
material  achieved  by  the  present  method  is  approximately  the  same  as 
in  Ref.  2,  except  for  the  spar  caps.  The  optimum  design  reported  in 
Table  5.16  is  obtained  in  15  design  iterations  for  both  the  cases. 

Also,  it  is  noted  that  for  Case  I,  a usable  design  is  found  at  the  10th 
iteration  with  a weight  of  2469.36  lbs.,  which  is  within  0.3%  of  the 
minimum  achieved.  For  Case  II,  a usable  design  is  found  at  the  7th 
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iteration,  with  a weight  of  2459.52  lbs,  which  is  within  0.7%  of  the 
minimum  achieved. 

Table  5.17  gives  the  cost  function  history  for  both  cases.  It  is 

1 2 

noted  that  components  of  fib  and  fib  associated  with  design  variables 
for  elements  1 to  10,  19  to  20  (Truss),  49  to  60  (CST) , and  5 to  10  (SSP) 
for  Case  II;  and  49  to  60  (CST),  and  5 to  10  (SSP)  for  Case  I are  equal 
to  zero  after  the  8th  iteration  for  Case  II  and  after  the  9th  iteration 
for  Case  I.  Hence,  in  the  subsequent  iterations  these  design  variables 
are  kept  fixed.  Computational  times  are  given  in  Table  5.18.  Table 
5.19  gives  a comparison  of  results.  The  cost  function,  as  compared  to 
Ref.  3,  is  practically  the  same  for  both  cases.  It  is  the  same  for 

Case  I,  as  compared  to  Ref.  2,  but  about  0.9%  less  for  Case  II.  With 

the  present  formulation,  the  total  CPU  time  is  less  by  factors  of  2.1 
and  6.4  for  Case  I,  as  compared  to  Refs.  2 and  4,  respectively.  For 
Case  II,  it  is  less  by  factors  of  3.1  and  9.1,  as  compared  to  Refs.  2 
and  4.  A comparison  of  CPU  time  per  design  iteration  indicates  that 
there  is  a reduction  by  factors  of  4.4  and  6.5  for  the  Case  I;  and  by 

factors  of  6.0  and  9.4  for  the  Case  II,  as  compared  to  Refs.  2 and  4, 

respectively. 

Lastly,  the  swept  wing  is  divided  into  five  substructures  by 
partitioning  it  at  nodes  9 to  12,  17  to  20,  25  to  28,  and  33  to  36. 

Nodes  41  to  44  are  also  considered  as  boundary  nodes.  With  this 
partitioning,  there  is  considerable  linking  of  design  variables  across 
substructure  boundaries.  The  optimum  results  obtained  are  the  same  for 
both  Cases  I and  II.  The  total  CPU  time  is  also  the  same.  This  is 
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TABLE  5.18.  COMPUTATIONAL  TIMES  IN  SECONDS  FOR 
150  (130)  ELEMENT  SWEPT  WING 


Problem 

Total  CPU  Time  for 

Total  CPU 
Time 

(2)+(3)+(4) 

Problem  Setup 

Stress-Ratio  Design 

r — 

* 

Optimal  Design 

1 

2 

3 

4 

5 

Case  I 

3.101 

9.800  (4)  + 

69.118  (15)$ 

82.019 

Case  II 

3.717 

. 

10.584  (4)+ 

74.067  (15)* 

88.368 

* 

Includes  CPU  time  required  for  analysis,  sensitivity  analysis, 
and  optimization. 

f 

Number  of  stress-ratio  design  Iterations  performed. 

^Number  of  design  Iterations  performed. 
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due  to  the  fact  that  the  total  number  of  calculations  with  two  and 
five  substructures  is  approximately  the  same  (see  Section  4.5). 

5.4  Discussion  and  Conclusions  Based  on 
Results  of  Chapter  5 

In  this  section  results  of  all  examples  are  discussed  and  compared 
with  other  solutions  available  in  the  literature  [2-8]. 

The  first  example,  the  18  element  wing  box  beam,  is  a small  scale 
structure  and  is  used  for  developing  the  computer  program  for  idealized 
wing  structures.  The  second  example  is  of  medium  scale,  whereas  the 
third  problem  is  fairly  large. 

The  optimum  weights  for  all  examples,  obtained  by  the  present 
method,  are  the  same  as  obtained  in  Refs.  2 to  4,  except  for  Example 
5.3  (Case  II),  in  which  the  final  weight  achieved  is  0.9%  less  than  in 
Ref.  2.  Further,  as  compared  to  Refs.  7 and  8,  there  is  a reduction 
of  about  9%  in  the  cost  function  for  Example  5.1  (Case  I). 

The  computing  times  with  the  present  method  are  significantly 
better  than  given  in  Refs.  2 to  4.  This  fact  is  also  observed  in 
Chapter  4.  Computing  times  in  Refs.  7 and  8 are  not  available  and  thus 
cannot  be  compared.  With  the  present  formulation,  there  is  a reduction 
in  CPU  times  per  design  iteration  by  factors  of  4.3,  4.6,  4.4,  and  6.0 
for  both  the  cases  of  Examples  5.1  and  5.3,  respectively,  as  compared 
with  CPU  times  reported  by  Schmit  and  Miura  [2].  In  comparison  with 
Rizzi  [4],  there  are  reductions  by  factors  of  3.4,  6.4,  6.5,  and  9.4 
for  Case  I of  Example  5.1,  Example  5.2,  and  both  cases  of  Example  5.3, 
respectively. 


r ) 30 

\ 

It  should  be  noted  that  in  the  present  method,  as  applied  to 

I 

k 

idealized  wing  structures,  selection  of  a proper  weighting  multiplier 
w^  associated  with  ith  design  variable  is  essential  for  fast  convergence. 

In  the  present  examples  this  multiplier  is  monitored  and  adjusted  to 
obtain  final  solutions.  As  a test  example,  side  computation  for  the 
Case  II  of  150  (130)  Element  Swept  Wing  was  made  keeping  the  weighting 
multiplier  to  be  unity  for  all  design  variables.  The  same  optimum 
as  reported  in  Table  5.16  was  achieved  in  40  iterations  as  compared 
to  15  iterations  reported  in  Table  5.17.  This  shows  that  proper  selec- 
tion of  a weighting  matrix  is  essential  for  rapid  convergence. 

Numerical  results  presented  in  this  chapter  indicate  that  the 
present  formulation  is  quite  efficient  to  serve  as  a practical  struc- 
tural synthesis  method.  However,  more  research  work  is  needed  in  the 
area  of  selection  of  a proper  weighting  matrix. 
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CHAPTER  6 

OPTIMAL  DESIGN  OF  FRAMED  STRUCTURES 


6.1  Introduction 

In  this  chapter,  the  optimal  design  method  developed  with  sub- 
structuring in  Chapter  2 is  applied  to  framed  structures.  Framed 
structures  are  found  quite  frequently  in  practical  situations, 
including  some  buildings,  transmission  towers,  bridges,  airplane 
hangers,  and  wing  structures.  Members  of  a framed  structure  usually 
consist  of  beams  and  columns  and  the  connections  between  them  are  classi- 
fied according  to  their  rotational  characteristics  as  rigid  (Type  I) , 
flexible  (Type  II),  and  semi-rigid  (Type  III)  [54].  In  the  present 
work,  optimal  design  of  plane  and  general  three-dimensional  framed 
structures  using  hollow  circular  sections  is  considered,  under  the 
assumption  that  the  connections  between  members  are  Type  I . Such 
framed  structures  are  called  "rigid  frames". 

The  use  of  hollow  circular  sections  (hereafter  also  called  tubes) 
as  structural  components  is  not  new,  but  this  usage  has  increased  con- 
siderably in  recent  years.  In  the  past,  use  of  tubes  was  restricted  be- 
cause of  connection  difficulties.  These  difficulties  have  been  overcome 
with  the  development  of  fully  automatic  oxyacetylene  tube-cutting 
machines.  A cost  comparison  of  tubular  and  rolled  sections  indicates 
that  the  increased  unit  cost  for  tubular  members  is  usually  offset  by 
the  Increase  in  weight  required  for  rolled  shapes  [55].  In  many 
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Instances,  steel  tubing  is  unsurpassed  in  its  efficiency,  especially 
where  weight  as  well  as  appearance  becomes  an  important  design  con- 
sideration (arenas,  stadiums,  etc.),  or  where  weight  savings  may  be  .1 
very  substantial  economical  consideration  (bridge,  derrick,  tower 
cranes,  etc.).  For  a complete  discussion  on  hollow  circular  sections, 
the  reader  is  referred  to  Ref.  55. 

The  class  of  framed  structures  considered  herein  is  assumed  to 
have  preassigned  geometric  configurations  and  structural  materials. 

The  mean  radius  and  thickness  of  tubes  are  treated  as  design  variables. 
Quite  often  in  structural  problems,  it  is  possible  to  obtain  a single 
parameter  that  completely  characterizes  a member,  thus  allowing  one  to 
treat  only  one  design  variable  for  each  member.  This  point,  however, 
is  not  central  to  the  present  application  and  two  design  variables 
are  used  for  each  member.  It  should  also  be  noted  here  that  a clear 
distinction  is  made  between  the  elements  and  the  members  of  a frame, 
in  order  to  adequately  represent  the  static  behavior  of  the  structure 
with  long  members.  In  other  words,  a member  of  the  frame  may  be 
divided  into  many  elements  for  the  purpose  of  structural  analysis. 
However,  only  uniform  members  are  treated  in  optimal  design.  There- 
fore, all  elements  of  a member  have  the  same  design  variables. 

Provision  of  design  variable  linking  is  also  incorporated. 

In  modeling  the  structure,  points  of  application  of  concentrated 
loads  are  treated  as  nodal  points.  Constraints  are  imposed  on  member 
stresses,  displacements,  and  design  variables.  Multiple  loading 
conditions  for  the  structure  are  treated  in  the  same  way  as  explained 
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in  Chapter  3.  Equations  developed  for  structural  analysis  in  Chapter 
2 are  also  used  here.  The  stiffness  matrix  for  a frame  element  is 
well  documented  in  the  literature  and  is  not  derived  here.  For  a 
complete  derivation,  the  reader  is  referred  to  Ref.  34. 


6.2  Reduction  of  the  Method  to  Design  of  Fr ame s 
The  ODPS  for  design  of  frames  is  defined  as  follows:  find  design 
variables  associated  with  each  member  of  the  frame  such  that  its 
weight  is  minimized  and  the  state  equations  (2.3-2)  and  (2.3-7),  and 
constraints  on  stress,  displacement,  and  member  size  are  satisfied. 

The  cost  function  of  Eq.  (2.4-1)  can  be  expressed  as: 

L NG(r)  NM(i) 

J (b)  « 2v  l l l p R t l (6.2-1) 

r=l  i=*l  j=l  J 

where 


and 


mean  radius  of  members  of  the  ith  group, 
thickness  of  members  of  the  ith  group, 
length  of  jth  member. 


b - {R-jl*  R2,'**’RNG(1),  C1  C2* ' ' ' CNG(1)  ’ * " * CNG(L)} 


(6.2-2) 


In  the  subsequent  three  sections,  treatment  of  constraints  is 
presented,  and  in  the  last  section  results  for  three  example  problems 
are  presented. 
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6.3  Stress  Constraint 

The  members  of  a framed  structure  are  subjected  to  a state  of 
combined  stress  due  to  axial  force,  bending  moment,  shear  force, 
and/or  torsional  moment  (Fig.  6.1).  Thus,  the  effect  of  combined 
stresses  must  be  considered  in  implementing  stress  constraints  of 
either  Eq.  (2.4-3)  or  Eq.  (2.4-4).  Further,  the  nature  of  the  axial 
force  distinctly  divides  the  design  criterion  for  a member  in  the 
following  two  categories: 

Case  1:  member  subjected  to  direct  tensile  force  with  (or  without) 

lateral  loads. 

Case  2:  member  subjected  to  direct  compressive  force  with  (or 

without)  lateral  loads. 

In  the  present  work,  members  are  required  to  satisfy  a design 
criterion  based  on  Von  Mises  equivalent  stress,  which  is  also  used 
in  Chapter  5.  Using  this  approach,  the  stress  constraint  for  an 
element  can  be  written  as: 

s 2 2 1/2 

<*'  (x^ » X2 » » zg » z j » b)  = ot(o  + 3t  } - 1.0  < 0 , ) 

( (6.3-1) 

for  all  , x2,  Xj  . / 
o 

where  a * l/°a,  °a  = , a = yield  stress,  °a  = allowable  stress, 

r b y 

FS  = factor  of  safety,  and  cr  and  t are  direct  and  shear  stresses, 
respectively.  It  may  be  pointed  out  here  that  o and  t are  assumed  to 
be  continuous  functions  of  spatial  coordinates  x , x^,  and  x^,  state 
variables  z and  z , and  design  variable  b.  Further,  in  this  research 


the  thickness  t of  the  tube  is  assumed  to  be  small  in  comparison  with 
the  mean  radius  R. 

In  Case  (2),  members  of  a framework  behave  as  beam  columns. 
Therefore,  a rational  design  procedure  should  take  into  account  end 
restraints.  There  are  many  approaches  for  the  design  of  restrained 
frame  members.  For  complete  detail  the  reader  is  referred  to  the  Ref. 
16.  In  this  dissertation,  Winter's  approach  [57]  which  is  a current 
design  practice  (for  restrained  beam-columns)  in  the  United  States 
and  Canada  is  used.  This  design  procedure  is  simple,  conservative, 
and  at  the  same  time  considers  the  strength  afforded  by  the  end 
restraints  of  adjacent  members.  Winter's  approach  is  to  check  beam- 
columns,  using  the  member  effective  length.  It  also  requires  one  to 
check  the  maximum  9tress  in  a member,  in  order  to  make  certain  that 
it  does  not  exceed  the  basic  allowable  value.  The  latter  requirement 
is  treated  in  the  same  way  as  for  Case  (1).  For  the  first  requirement, 
Gerard  and  Becker  [58]  present  interaction  formulas  together  with 
experimental  verification  under  various  combined  loading  conditions; 
including  1)  axial  compression  and  bending,  2)  axial  compression  and 
torsion,  3)  axial  compression,  bending  and  torsion,  4)  bending  and 
torsion,  and  5)  bending  and  transverse  shear.  Although  a slightly 
different  interaction  equation  is  proposed  for  each  case.  Schilling 
[59]  suggests  that  a single  formula  can  be  conservatively  applied  to 
all  the  above  mentioned  loading  combinations.  For  details  the  reader 
is  referred  to  Refs.  56  and  59. 

Using  Schilling's  approach,  an  additional  stress  constraint  for 


Case  (2)  is  written  as 


* 
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(x1,X2,X3,ZB,ZI’b)  ; a 


1.0  < 0 , 


for  all 


(6.3-2) 


where  a = normal  stress  due  to  axial  force  and/or  bending  moment 
at  point  (x^,  x2,  x^) 

t * shear  stress  due  to  shear  force  and/or  twisting  moment 
at  a point  (x^  x2>  x3> 

acr  = critical  buckling  stress  for  axial  compression  (or 
bending)  alone,  and 

= critical  shear  buckling  stress  for  torsion  (or  trans- 
verse shear)  alone. 


The  values  of  a and  t for  a member  are  given  in  Appendix  3. 
cr  cr 

As  pointed  out  earlier,  a member  of  the  frame  may  be  divided  into 

many  elements  for  the  purpose  of  structural  analysis.  On  the  other 

hand,  the  critical  compressive  stress  for  all  elements  making  up  a 

member  is  the  same  because  the  axial  loads  ate  assumed  to  be  acting 

only  at  the  end  of  a member.  Hence,  the  effective  length  factor  k, 

required  in  the  computation  of  a , is  to  be  calculated  for  each 

cr 

member  of  the  frame.  In  rigid  frames,  two  cases  of  interest  must  be 
considered:  (a)  frames  without  sidesway,  and  (b)  frames  with  sidesway. 

The  value  of  k for  Cases  (a)  and  (b)  are  obtained  by  solving  the 
following  transcendental  equations,  respectively  [54]: 


r 
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"-(i) 

- [H  + V (ij  * 2]  cos  ( k ) + ^ " 0 

arid 

[Vb  (s)‘  - “ ] Sin(i)-  6<Ga  + Gb)  JcoB(J).°  (6.3-4) 

The  subscripts  A and  B on  G in  Eqs.  (6.3-3)  and  (6.3-4)  refer  to  the 
joints  at  the  ends  of  the  member  under  consideration.  The  parameter 
G is  given  by 


(6.3-5) 


where  £ indicates  a summation  for  all  members  rigidly  connected  to 
that  joint  and  lying  in  the  plane  in  which  buckling  of  the  member  is 


being  considered,  Ic  is  the  moment  of  inertia,  and  is  the  corre- 
sponding unbraced  length  of  the  compressed  members.  Also,  1 Is  the 
moment  of  inert  La,  and  i the  corresponding  unbraced  length  of  the 


restraining  members,  and  I and 

c 


I 

8 


are  taken  about  axes  perpendicular 


to  the  plane  of  buckling.  In  the  present  work,  the  secant  method  of 


nonlinear  algebraic  equations  is  used  in  finding  the  lowest  root  of 
Eqs.  (6.3-3)  and  (6.3-4). 


I 

J 
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6.3.1  Treatment  of  Stress  Constraint 
Functions  in  the  stress  constraints  of  Eqs.  (6.3-1)  and  (6.3-2) 
are  assumed  to  be  continuous  in  the  spatial  coordinates  x^,  and  x^, 
state  variables  z and  z , and  design  variable  b.  It  is  required  that 
maximum  normalized  equivalent  stresses  for  all  members  of  a frame  be 
within  the  specified  limit.  In  general,  the  stress  constraint  for 
Cases  (1)  and  (2)  of  the  previous  section  can  be  expressed  as  (Fig. 

6.1) 

Max  $S(x .,0,z _,z  ,b)  < 0 (6.3-6) 

12  1 B I - 

where 

(2  = {x  ,0  I 0 < Xx  < fc,  0 < e < 2tt}  (6.3-7) 

and 

6 = 0(x2,x3)  (6.3-8) 

For  notational  convenience,  Eq.  (6.3-6)  represents  a stress  constraint 
for  a typical  frame  element.  Therefore,  no  subscript  for  the  con- 

g 

stralnt  function  $ will  be  used. 

The  treatment  of  a max-value  constraint  defined  by  Eq.  (6.3-6) 
in  the  general  minimization  problem  is  mathematically  quite  tedious 
and  computationally  quite  time  consuming.  This  type  of  constraint 
function  may  not  even  be  differentiable.  There  are  two  basically 
different  techniques  [33,47]  that  can  be  used  in  treating  this  type 


of  constraint,  which  are  briefly  reviewed  here. 
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Recently,  Kvak  [47]  has  treated  a very  general  class  of  such 
min-max  problems.  He  has  developed  and  evaluated  several  algorithms 
for  such  problems.  His  First  Order  Algorithm  is  quite  efficient  and 
reliable  [47].  Briefly,  the  First  Order  Algorithm  is  a two-step 
procedure.  For  the  OOPS,  in  the  first  step,  points  that  give  max-value 
to  the  constraint  function  of  Eq.  (6.3-6)  are  solved  for  a given  design 
variable  b(0)  and  state  variables  z<0)  and  z{0).  This  is  called  the 
inner  optimization  problem,  i.e.,  find  x^O  e ft  such  that  *8  i6 
maximized.  Then,  the  constraint  Eq.  (6.3-6)  is  replaced  by  constraints 
at  these  points  only.  In  the  second  step,  points  that  give  max-value 
to  constraint  of  Eq.  (6.3-6)  are  kept  fixed  and  the  general  minimiza- 
tion problem,  known  as  an  outer  problem,  is  solved.  The  entire  pro- 
cecure  is  repeated  until  desired  convergence  is  obtained.  There  are 
also  procedures  which  account  for  changes  in  the  solution  points  of 
the  inner  problem,  known  as  higher  order  or  hybrid  methods  [47]. 

Feng  [33],  on  the  other  hand,  first  converted  the  max-value  con- 
straint of  Eq.  (6.3-6)  into  the  pointwise  constraint 

g 

<p  (x1,Q,zB,zI,b)  < 0 , for  all  x^.0  e ft  (6.3-9) 

and  then  replaced  the  constraint  Eq.  (6.3-9)  by  the  equivalent 
functional  constraint 
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g 

Since  4>  (x^.O.Zg.Zj.b)  is  continuous  with  respect  to  x^  and  6,  the 

integrand  in  Eq.  (6.3-10)  is  also  continuous  with  respect  to  x^  and  G 

and  it  is  always  non-negative.  If  any  constraint  in  the  expression 

(6.3-9)  is  violated,  the  corresponding  functional  constraint  in  Eq. 

(6.3-10)  is  also  violated,  and  vice-versa.  Thus,  constraint  Eqs. 

(6.3-9)  and  (6.3-10)  are  equivalent.  In  Feng's  approach,  the  domain 

U of  the  functional  constraint  is  discretized  by  introducing  grid 

points  in  the  x^  direction  and  the  0 direction  of  the  element  (Fig. 

6.1).  The  integrand  of  Eq.  (6.3-10)  is  calculated  at  each  of  the 

grid  points  and  violations  in  the  constraint  are  noted.  In  order  to 

correct  this  constraint  error,  design  sensitivity  analysis  of  the 

functional  must  first  be  developed.  This  involves  calculation  of  an 

adjoint  matrix  at  each  of  the  grid  points  where  the  constraint  is 

violated.  Finally,  numerical  integration  must  be  done  to  calculate 
d<l>s 

-f—  . This  procedure  must  be  followed  for  each  functional  constraint, 
db 

Selection  of  a technique  to  treat  the  max-value  constraint  of 
the  Eq.  (6.3-6)  appears  to  be  problem  dependent.  In  the  case  of 
frame  problems,  Feng's  approach  will  be  numerically  inefficient, 
because  of  the  large  number  of  grid  points  that  will  have  to  be 
introduced  to  discretize  the  domain  (2.  On  the  other  hand,  Kwak's 
First  Order  Algorithm  may  be  quite  efficient,  since  analytic  solution 
of  the  inner  problem  for  frames  is  possible.  Thus,  Kwak's  approach 
is  adopted  in  this  dissertation.  Details  of  stress  constraint  treat- 
ment for  plane  and  space  frame  members  are  presented  separately  in 
the  following  two  subsections. 
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6.3.2  Stress  Constraint  For  Plane  Frames 
Figure  6.2  shows  a plane  frame  elment  and  its  local  coordinate 
system,  with  the  sign  convention  to  be  used  on  element  forces  P^'6 
and  deformations  u^'s.  The  element  is  assumed  to  be  subjected  to  ln- 
plane  loading  only.  As  pointed  out  earlier,  the  points  of  application 
of  concentrated  loads  are  treated  as  nodes,  so  there  is  no  disconti- 
nuity in  the  shear  force  or  bending  moment  in  an  element. 

Consider  a section  at  x^  and  let  axial  force,  shear  force,  and 
bending  moment  at  this  section  be  F^x^),  F^x^,  and  F^x^), 
respectively.  Since  no  axial  load  is  applied  inside  the  element, 

F^  is  assumed  to  be  constant  over  the  length  of  the  element.  Further, 
the  thickness  t of  the  tube  is  assumed  to  be  small  in  comparison  with 
mean  radius  R.  Thus,  the  shear  stress  distribution  can  be  assumed 
constant  over  the  thickness  t.  Now  the  direct  stress  oCx^.x^)  and 
shear  stress  t(x^,x^)  at  a point  (x^,x2)  can  be  respectively  expressed 
as 


a(x1,x2) 


W 


(6.3-11) 


and 


t(x1>x2) 


F2(x1)Q(x2) 

1 b (x2) 


(6.3-12) 


where 


Section  at  a-a 


Figure  6.2.  A Typical  Plane  Frame  Element 


• .. 
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F = P 
1 1 


(6.3-13) 


F2(X1)  = P2  + '"‘l 


(6.3-14) 


F3^X1^  P3  P2X1  ~ 2 


(6.3-15) 


Here,  A is  the  cross-sectional  area,  I is  the  moment  of  inertia  of  the 
cross-section  about  x^axis,  b(x2)  * 2t  for  thin-walled  tubes,  and  the 
parameter  Q(x„)  associated  with  shear  stress  can  be  expressed  as 


[refer  to  Appendix  3]: 


Q(x2)  - 2R  t jl 


-(¥! 


2\  1/2 


(6.3-16) 


Substituting  the  expressions  for  Q(x_)  and  b(x_)  in  Eq.  (6.3-12),  one 


obtains 


2 i 1/2 


(6.3-17) 


The  expressions  for  direct  and  shear  stresses  at  any  point  (x^,x2)  ate 
now  known.  Using  the  First  Order  Algorithm  of  Kwak  [47],  stress 
constraints  for  Cases  (1)  and  (2)  of  Section  6.3  are  now  considered 
separately. 

6. 3. 2.1  Case  (1) 

The  stress  constraint  Eq.  (6.3-1)  for  a plane  frame  element 
i educes  I n 
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1/2 

*S(x1,x2,zB,zI,b)  = a |a2(xlfx2)  + 3t2(Xi,x2>| 

- 1.0  < 0 , for  all  x^  and  x2  (6.3-18) 

where  o^.x^  and  t^.x^  are  given  by  Eqs.  (6.3-11)  and  (6.3-17), 
respectively.  Equation  (6.3-18)  can  be  rewritten  as  a max-value  con- 
straint (see  Eq.  (6.3-6))  and  the  domain  0 for  the  constraint  is  given 


0 = (x1,x2  | 0 < Xj^  < «.  , ~ R < x2  < R} 


(6.3-19) 


This  constraint  function  has  a continuous  derivative  in  the  neighbor- 
hood of  a point  (x  ,x2)  which  is  in  the  interior  of  domain  fl  . The 
necessary  condition  for  a stationary  point  (x^.j^)  of  the  constraint 
function  defined  by  Eq.  (6.3-18)  is 


’ a{°2<V 


x2)  + 3t2(x1,x2)| 


- 1.0  = 0 


(6.3-20) 


Thus,  from  Eq.  (6.3-20)  one  obtains,  after  simplification, 


FI 

*1>  (F^(x^)  + 3wR2}  x2  + — x2  - 3wR^  «■  0 


(6.3-21) 


A OR^O^)  - F2^)}  x2  - IF^Sj)  = 0 


(6.3-22) 


Equation  (6.3-21)  implies  that  either  F„(x, ) » 0,  i.e.  x,  » if 

* 1 1 w 


w / 0,  or 


I 


l/*6 


2 -2  _ 4 

(F^x^)  + 3wR  } x2  + — — x2  - 3wR  = 0 


Solving  Eq.  (6.3-22),  one  obtains 


- , JsLZiW__ 

2 3R2F2(51)  - FjUj) 


R < x2  < R 


(6.3-23) 


(6. 3-24) 


Substituting  the  expressions  for  x2,  F^x^),  and  F^x^)  from  Eqs. 
(6.3-24),  (6.3-14),  and  (6.3-15),  respectively,  in  Eq.  (6.3-23);  one 
obtains  an  8th  degree  polynomial  in  x^: 


-8  -7  , 

Vl  + a7xl  + 


+ a,x  + a = 0 
1 o 


(6.3-25) 


where  the  coefficients  (a  , i = 0 to  8)  are  given  in  Appendix  3. 

P2  . 

Using  the  permissible  zeroes  of  Eq.  (6.3-25)  and  x,  = if  w # 0, 

1 w 

corresponding  values  of  x2  from  Eq.  (6.3-24)  are  calculated.  Thus, 
one  obtains  a set  of  stationary  points  at  which  the  stress  constraint 
of  Eq.  (6.3-18)  needs  to  be  evaluated. 

In  case  w = 0,  the  8th  order  polynomial  defined  by  Eq.  (6.3-25) 
reduces  to  a linear  equation  (see  Appendix  3 for  coefficient  a^'s) 
and  the  stationary  point  along  the  x^-axis  is  given  by  = (P^/P2), 
which  implies  that  the  point  at  which  the  bending  moment  is  zero  (refer  to 
Eq.  (6.3-15))  is  a stationary  point  and  from  Eq.  (6.3-24)  one  obtains 
x2  - 0.  Thus,  in  this  case,  the  necessary  condition  gives  only  one 


stationary  point.  It  may  be  noted  that  whether  w * 0 or  w # 0,  one 
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Substituting  the  values  of  F^x^,  and  F2(x  ) in  Eq-  (6-3-28)» 
after  simplification,  one  obtains 

(f  “ 6w2^2)xJ  + (?2  + 28wP2x2)x1  + 6^2P2  - P-j  = 0 (6.3-30) 

Now  the  zeroes  of  quadratic  Eq.  (6.3-29)  are  obtained  and  then,  for 
each  permissible  value  of  x2>  the  quadratic  Eq.  (6.3-30)  is  solved 
for  x^.  Thus,  one  obtains  a set  of  stationary  points.  In  case  w = 0, 
then  Eq.  (6.3-29)  yields  x2  = 0 and,  from  Eq.  (6.3-30),  x^  = 

Thus,  the  stationary  point  for  the  stress  constraint  of  Case  (2), 
when  w = 0,  is  the  same  as  obtained  in  Section  6. 3. 2.1.  The  remainder 
of  the  procedure  for  obtaining  the  global  maximum  of  the  stress 
constraint  defined  by  Eq.  (6.3-26)  is  exactly  the  same  as  in  Section 
6. 3. 2.1. 


6.3.3  Stress  Constraint  for  Space 
Frames 

A typical  three  dimensional  frame  element  and  its  local  coordi- 
nate system,  with  the  sign  convention  to  be  used  for  element  forces 
P^'s  and  deformations  u^'s,  is  shown  in  Fig.  6.1.  The  element  is 
assumed  to  be  subjected  to  only  nodal  loads  that  pass  through  the 
shear  center.  Hence,  x^  = 0 and  x^  = i are  the  only  critical  points 
along  the  x^-axis  at  which  one  needs  to  impose  the  stress  constraint. 

Consider  a section  at  either  x^  * 0 or  x^  = t and  let  axial 
force  F^,  shear  forces  F2  and  F^,  twisting  moment  F^,  and  bending 
moments  F,.  and  F^  be  the  force  components  at  that  section.  As  in 
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Section  6.3.2,  the  thickness  t of  the  hollow  circular  section  is 
assumed  to  be  small  in  comparison  to  mean  radius  R.  The  direct  stress 
o(6)  and  shear  stress  t(0)  at  a point  P (see  Fig.  6.1)  are  expressed  as 

F,  FCR  F.R 

a(e)  * — + -5-  cos  9 + — sin  0 (6.3-31) 

A 1 *■ 


and 


t(0) 


2 2 
F.R  F R F R 

+ — - — cos  0 + — ~ — sin  0 


(6.3-32) 


1 


As  in  Section  6.3.2,  stress  constraints  for  Cases  (1)  and  (2)  of 
Section  6.3  are  considered  separately. 

6.3. 3.1  Case  (1) 

The  stress  constraint  of  Eq.  (6.3-1)  for  a space  frame  element 
reduces  to 

2 2 1/2 

4>S(0,z  ,zT,b)  = a { o (0)  + 3t  (0))  - 1*0  < 0 , 

D 1 

for  all  9 . (6.3-33) 

where  o(0)  and  t(0)  are  given  by  Eqs.  (6.3-31)  and  (6.3-32),  respec- 
tively. Equation  (6.3-33)  can  be  rewritten  as  a max-value  constraint 
(see  Eq.  (6.3-6)),  with  the  domain  « given  as 

£2  - {0|  0 < 0 < 2ir}  (6.3-34) 

Following  the  same  approach  as  in  the  previous  section,  after  simpli- 
fication, one  obtains  the  following  necessary  condition  for  the 
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stationary  point9  (0)  of  constraint  Eq.  (6.3-33): 

a.  cot  20  + a_  sec  0 + a,  cosec  0 + a * 0 (6.3-35) 

3 2 1 o 

where  the  coefficients  (a^,  i = 0 to  3)  are  given  in  Appendix  3. 

From  the  transcendental  Eq.  (6.3-35),  one  obtains  a set  of 
stationary  points  at  which  the  stress  constraint  of  Eq.  (6.3-33)  needs 
to  be  evaluated.  Finally,  one  obtains  a global  maximum  of  this  con- 
straint, which  should  be  satisfied. 

6. 3. 3. 2 Case  (2) 

For  a space  frame,  the  buckling  constraint  given  by  Eq.  (6.3-2) 
reduces  to 

2 

4>S (0 , zD , z ,b)  i ^ - 1.0  < 0 , for  all  0.  (6.3-36) 

B 1 °cr  t 

cr 

Using  the  same  approach  as  in  Section  6. 3. 3.1,  one  obtains  the  follow- 
ing transcendental  equation  for  stationary  points  (0)  at  which  the 
above  buckling  constraint  needs  to  be  evaluated: 

a^  cot  20  + a ^ sec  0 + a^  cosec  0 + Sq  = 0 (6.3-37) 

where  the  coefficients  (a^,  i = 0 to  3)  are  given  in  Appendix  3. 

Finally,  one  obtains  a global  maximum  of  this  constraint  which  should 


be  satisfied. 
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6.4  Displacement  Constraints 

In  earlier  work  [28],  displacement  constraints  for  a frame  were 
imposed  at  the  nodal  points  only.  If  the  constraint  was  to  be  imposed 
at  an  interior  pcint  of  a member,  due  to  practical  and/or  design 
considerations,  then  that  point  was  treated  as  a node.  In  the  above 
treatment,  the  displacement  constraint  does  not  depend  upon  spatial 
coordinates  x^,  x^,  or  x^,  so  it  can  be  implemented  in  a way  exactly 
similar  to  the  case  of  trusses  or  idealized  wing  structures  (see 
Chapters  4 or  5) . For  large  framed  structures,  this  approach  greatly 
increases  the  dimensionality  of  the  problem,  as  well  as  computing  time. 

In  the  present  work,  displacement  constraints  are  imposed  on  the 
maximum  value  of  displacement  in  the  x^  and  x^  directions,  in  the 
interior  of  the  element  (see  Fig.  6.1),  along  with  constraints  on 
nodal  displacements.  The  max-value  displacement  constraint  is  treated 
like  the  stress  constraints  of  Section  6.3,  since  both  constraints 
are  of  the  same  form.  For  implementation  of  the  constraints  on  nodal 
displacements  the  reader  is  referred  to  (developments  of)  Chapter  4. 

Neglecting  deformation  due  to  shear,  the  displacement  constraint 
in  a lateral  direction  is  imposed  along  the  x^-axis,  and  for  a typical 
element  is  written  as 


i 


8 - 5<V 

$ (x],,zB,zI,b)  = 6 — - — - 1.0  < 0 , for  all  x^  (6.4-1) 

where  C(x^)  is  the  lateral  displacement  at  a point  x^,  C is  the 

allowable  displacement,  and  8 is  an  amplification  factor  [60].  In 
- P 

general,  6 will  depend  on  the  ratio  - — (where  P is  the  applied 

cr 
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compressive  axial  force  and  is  the  Euler  critical  load)  and  the 

type  of  lateral  load.  In  the  present  work  8 is  taken  to  be  unity  for 

tensile  axial  force,  which  is  somewhat  conservative.  For  compressive 

- P 

axial  force,  a simplified  expression  for  8 - 1/(1-  p — ) is  used  with 

P cr 

good  accuracy,  if  the  ratio  - — is  not  large.  For  complete  details 

cr 

the  reader  is  referred  to  Ref.  60.  Further,  Eq.  (6.4-1)  is  rewritten 
as  a max-value  constraint 

V - S<xl>  i 

Max  } 8 — -’.Of  <0  (6.4-2) 

a I ca  > 

where 

n = {%1  | o < x1  < a}  (6.4-3) 


In  subsequent  sections,  the  displacement  constraint  is  considered  for 
plane  frames  and  space  frames,  separately. 


6.4.1  Displacement  Constraint  for 
Plane  Frame 

Using  the  Euler-Bernoulli  equation,  the  general  expression  for 
elastic  curve  of  a plane  frame  element  (Fig.  6.2)  can  be  written  as 


w 


w x. 


24EI 


. 1 3 

+ 6 C3X1 


. 1 2 
+ 2 C2X1 


C1X1  + C0 


(6.4-4) 


where  is  displacement  in  the  X2~direction  and  c^,,  i = 0 to  3 
are  constants  of  integration  (see  Appendix  3).  Using  the  necessary 
conditions  of  Section  6.3,  one  obtains  a 3rd  order  polynomial  in 

X1  xi  < *■> 
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-3 

1 W X1  1 -2  -2 
6 El  + 2 C3X1  + C2X1  + ci  “ 0 


(6.4-5) 


The  zeroes  of  the  above  cubic  equation  can  easily  be  obtained  as  a 
closed  form  solution  [61].  Further  treatment  of  this  constraint  is 
exactly  the  same  as  in  Section  6.3. 


6.4.2  Displacement  Constraint  for 
Space  Frame 

Let  S2(x;l)  and  ^3^xi^  be  displacements  in  the  x2  and  directions, 
respectively.  The  total  displacement  £(x^)  can  be  expressed  as 

«<*!>  - U^)  + ^(Xj)}172  (6.4-6) 


Using  the  same  approach  as  in  previous  section,  the  general  expressions 
for  elastic  curves  in  the  x2  and  x^  directions  can  be  expressed  as 
(Fig.  6.1) 


’2V  1 


6 3 1 


"2  1 


1 1 


(6.4-7) 


and 


^3^X1^  " 6 d3xl  + 2 d2Xl  + dlXl  + d0 


(6.4-8) 


where  c^  and  di  (i  - 0 to  3)  are  constants  of  integration  and  are 
given  in  Appendix  3.  Using  the  necessary  conditions  for  stationary 
points  of  the  displacement  ^(Xj),  one  obtains  a 5th  order  polynomial 
in  x^ 


~r  ■ 


V. 


I,.,-- — 
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Vi  + 


+ a^ 


ao  = ° 


(6.4-9) 


where  the  coefficients  *^1  = 0 to  5)  are  given  in  Appendix  3.  One 
can  solve  this  5th  order  polynomial  for  stationary  values,  in  order 
to  implement  the  constraint. 


6.5  Design  Variable  Constraint 

In  this  research,  mean  radius  (R)  and  thickness  (t)  of  a hollow 
circular  member  are  treated  as  design  variables.  Thus,  the  design 
variable  constraints  of  Eq.  (2.2-5)  are  expressed  as 


RL  < R < RU  (6.5-1) 

and 

tL  < t < tU  (6.5-2) 

where  RU,  RL,  tU,  and  tL  are  the  upper  and  lower  bounds  on  R and  t, 
respectively.  These  bounds  on  design  variables  are  to  be  furnished 
by  the  designer.  In  the  present  work,  bounds  on  mean  radius  are 
provided  by  the  designer,  whereas  bounds  on  thickness  are  selected  as 


for  a thin  circular  section  | 62 ] and 


to  avoid  local  buckling  of  a member  [56]. 


(6.5-3) 


(6.5-4) 
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6.6  Design  Sensitivity  Analysis 

The  cost  function  J(b)  given  by  Eq.  (6.2-1)  depends  only  on  design 

variables  and  is  nonlinear.  Thus,  the  design  sensitivity  vector  of 

3j 

J(b)  is  easily  calculated  as  5J  = -ft  6b. 

The  stress  constraints  for  Cases  (1)  and  (2)  of  Section  6.3  at 
a stationary  point  are  continuous  functions  of  state  variables  z and 

D 

Zj.  and  design  variable  b.  The  violated  stress  constraints  at  a 
desired  point  are  first  expressed  in  terms  of  nodal  forces  by  using 
Eqs.  (6.3-13)  to  (6.3-15)  for  plane  frames,  and  Eqs.  (6.3-31)  and 
(6.3-32)  for  space  frames.  These  equations  are  then  expressed  in 
terms  of  nodal  displacements  in  a global  coordinate  system  by  using 
Eq.  (2.3-14).  The  sensitivity  coefficients  ~ , and  are 

dZg  dZj  db 

then  routinely  evaluated.  The  computations  are  somewhat  lengthy 
but  are  straightforward.  In  the  present  work,  the  effective  length 
factor  k,  defined  by  Eq.  (6.3-3)  or  (6.3-4),  is  treated  as  constant 
in  a particular  design  cycle.  This  assumption  is  reasonable,  since 
the  value  of  k does  not  change  appreciably  from  one  design  cycle  to 
another,  which  was  also  observed  in  Ref.  28.  However,  at  the  start 
of  each  design  cycle,  k values  for  members  subjected  to  axial  com- 
pression are  recomputed,  or  read  in. 

Lastly,  the  design  sensitivity  analysis  for  displacement  and 
design  variable  constraints  is  straightforward  and  can  easily  be 
performed.  Further,  bounds  on  design  variables  are  calculated  at  the 
beginning  of  each  design  cycle.  At  this  stage,  all  necessary  infor- 
mation is  available  to  assemble  the  matrix  A of  Eq.  (2.6-6). 




■ tnCSjftM*  * ' 


■I'S 
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6.7  Example  Problems 

A computer  program  based  on  the  algorithm  of  Chapter  2 and  the 
flow  diagram  of  Chapter  4 is  developed  for  optimal  design  of  plane  and 
space  frames.  The  following  structures  are  optimized  as  examples: 

(i)  Example  6.1:  One-Bay  Two-Story  Frame 

(ii)  Example  6.2:  Two-Bay  Six-Story  Frame 

(iii)  Example  6.3:  Helicopter  Tail  Boom  Structure 

The  first  two  examples  are  plane  frames  that  were  solved  in 
Refs.  28  and  41  by  treating  the  moment  of  inertia  of  each  member  as 
a design  variable.  In  the  present  formulation, ' the  mean  radius  (R) 
and  thickness  (t)  of  each  member  are  considered  as  design  variables. 
Further,  the  constraints  imposed  here  are  different  than  those  imposed 
in  Refs.  28  and  41.  Therefore,  direct  comparisons  of  final  results 
are  not  possible. 

The  last  example  is  a space  frame.  This  example  is  new  and  has 
not  been  solved  earlier.  A truss  idealization  of  this  structure  was 
optimized  in  Ref.-  63.  Here,  a frame  idealization  of  the  same  structure 
is  optimized  and  the  final  solution  analyzed. 

In  the  following  subsections  each  example  is  discussed  separately 
and  the  final  results  are  presented.  The  computation  times  reported 

I 

herein  are  for  an  IBM  360/65(H)  computer  (double  precision).  The 
value  of  e for  imposing  e-active  constraints  varies  from  1 to  5%  for 
all  example  problems.  The  weighting  matrix  is  selected  as  explained 
in  Chapter  5.  For  all  examples,  the  value  of  weighting  multiplier 
(w^)  for  mean  radius  (R)  varies  from  2 to  100,  whereas  for  thickness 
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(t)  it  is  taken  as  unity.  The  factor  of  safety  (FS)  associated  with 

Eq.  (6.3-1),  for  evaluating  allowable  stress  (oa) , is  selected  based 

on  the  recommendations  of  Ref.  64.  Further,  the  bound  (£a)  used  in 

imposing  the  displacement  constraint  in  the  interior  of  an  element 

l 

(see  Eq.  (6.4-1))  is  taken  as  -r^r  , where  £ is  the  length  of  the  member 

JoU  m 

to  which  the  element  belongs.  Finally,  the  allowable  nodal  displace- 
ments are  taken  in  accordance  with  Refs.  65  and  66,  such  that  the 
deflection  index  (ratio  of  deflection  to  height)  is  restricted  to 
0.002. 

6.7.1  Example  6.1:  One-Bay  Two-Story  Frame 
This  small  scale  frame  serves  as  a good  example  for  testing  the 
computer  program  based  on  the  present  algorithm.  Figure  6.3  shows 
the  dimensions  and  the  loading  conditions  for  this  structure.  In  the 
present  formulation  the  frame  is  divided  into  two  substructures  by 
partitioning  it  at  nodes  2,  6,  and  7.  Thus,  the  first  substructure 
consists  of  elements  1,  5,  6,  and  8 and  has  9 boundary  degrees  of 
freedom.  The  remaining  elements  are  in  the  second  substructure,  which 
has  6 boundary  and  9 interior  degrees  of  freedom.  Table  6.1  gives 
the  design  data  for  this  example.  It  may  be  noted  that  in  Refs.  28 
and  41  symmetry  in  the  structural  design  was  guaranteed  by  imposing 
an  additional  loading  condition.  In  the  present  work  design  variable 
linking  is  used  in  order  to  provide  symmetry  in  the  frame  (see  Table 
6.2).  This  feature,  as  pointed  out  in  Chapter  3,  is  more  efficient 
since  it  reduces  the  number  of  loading  conditions  as  well  as  design 


variables. 


3 


4 


/-6K/ft 
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TABLE  6.1.  DESIGN  DATA  FOR  ONE-BAY 
TWO- STORY  FRAME 


= 3 x 10^  ksi 
= 0.2836  lb/in.3 
= 36.0  ksi 
= 1.51 

= 16.610  in. 

= 0.450  in. 

= 1.0  in. 

= None 

!±0 . 7 2 in.  for  nodes  3 to  5 
±0.36  in.  for  nodes  2,  6,  & 7 
±0.67  in.  for  elements  3 to  6 

= 0.037 

Starting  weighting  multiplier,  = 0.04 

Number  of  loading  conditions  = 2 
(see  Fig.  6.3) 

The  above  data  is  same  for  both  the  cases. 


Modulus  of  elasticity 
Material  density 
Yield  stress 
Factor  of  safety 
Initial  value  of  mean  rad.,  R 
Initial  value  of  thickness,  t 
Lower  limit  on  mean  rad.,R 
Upper  limit  on  mean  rad.,R 

Displacement  limits  in 
and  directions 

Starting  reduction  ratio,  r 
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The  framed  structure  Is  optimized  under  two  separate  constraint 
conditions:  stress  constraint  only  (Case  I);  and  stress  and  displace- 

ment constraints  (Case  II).  Table  6.2  gives  the  final  designs,  values 
of  cost  function,  number  of  active  constraints,  maximum  constraint 
violations,  | | 6b^  | | ^ at  the  optimum,  the  maximum  value  of  | | 6b^ | | 2 » and 
the  total  CPU  time.  The  set  of  active  constraints  at  the  optimum  design 
includes : 

Case  I:  Buckling  constraint  for  elements  3,  5,  7,  and  9 under  loading 

condition  2,  and  minimum  element  thickness  for  all  elements. 
Case  II.  Buckling  constraint  for  element  8 under  loading  condition  2, 
displacement  constraint  at  node  3 in  x^-direction  under 
loading  condition  2,  and  minimum  element  thickness  for  all 
elements. 

The  total  CPU  time  reported  in  Table  6.2  includes  initial  set-up  time 
(0.707  sec  for  both  cases),  and  15  and  21  iterations  of  design  time 
(10.017  sec  and  14.135  sec)  for  Case  I and  II,  respectively. 

Table  6.3  gives  the  cost  function  history  for  both  cases.  As 
pointed  out  in  Section  6.6,  the  effective  length  factor  (k)  does  not 
change  appreciably  from  one  design  cycle  to  another.  This  fact  is 
also  observed  numerically  by  evaluating  k at  each  design  iteration. 

A side  computation  was  made  by  keeping  k fixed  after  2 design  iterations 
and  essentially  the  same  optimum  design  was  achieved.  The  CPU  time  per 
design  iteration  was  reduced  by  12%  for  this  run. 

The  Von  Mises  equivalent  stress  constraint  is  imposed  at  the 
relative  maximum  points  (x^.x^)  obtained  from  solutions  of  Eqs.  (6.3-25) 


' 1 


i 


TABLE  6.3.  COST  FUNCTION  HISTORY  FOR  ONE-BAY 
TWO- STORY  FRAME 


DI 

Number 


Weight  in  lbs. 


Case  I 


Case  II 


m 


W— •*-? 

*'*  • i *• 


15805.09 

8895.28* 

8054.98 
8432.11 
8137.61 
8369.07 
8642.53 
8741.71 
8863.67* 
8924.57 

8954.21 

8968.43 

8975.21 

8978.44 

8979.99 
8980.73+ 


15805.09 
8895.28* 

10534.66 

11383.10 
11213.80 
10820.25 
10420.79 
10391.42 

9902.82* 

10335.22 
8015.05 
8270.03 
9142.02 

10141.47 

10162.23 
10162.35* 
10175.35 
10176.06 
10169.92 
10166.22 

10166.48 
10166.63* 


■underlined  iteration  is  not  counted  as  it  is  used 
for  checking  the  solution  only. 

* 

Means  change  in  step  size. 

+ Design  iteration. 
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and  (6.3-24),  the  ends  of  the  element,  and  the  point  of  maximum  bending 
moment.  The  Von  Mises  equivalent  stress  constraint  was  not  active  at 
any  of  the  relative  maximum  points  obtained  from  the  necessary  condi- 
tions of  Eq.  (6.3-25),  during  the  entire  design  process.  Thus,  the 
computing  times  reported  in  Table  6.2  are  for  the  case  in  which  stress 
constraints  are  imposed  only  at  the  end  points  and  at  the  point  of 
maximum  bending  moment.  When  critical  points  arise  in  the  necessary 
condition  (solutions  of  Eqs.  (6.3-25)  and  (6.3-24))  and  the  Von  Mises 
equivalent  stress  constraint  is  checked  at  these  points,  the  computa- 
tional times  are  increased  by  30%. 

Based  on  the  foregoing,  the  polynomial  defined  by  Eq.  (6.3-25)  is 
not  solved  in  the  subsequent  examples  and  stress  constraints  defined 
by  Eq.  (6.3-18)  are  imposed  along  the  x^-axis,  at  the  ends  of  the 
elements  and  at  the  point  of  maximum  bending  moment.  The  remainder 
of  the  procedure  is  the  same  as  explained  in  Section  6.3.2.  Also,  the 
effective  length  factor  (k)  is  not  evaluated  in  each  design  iteration. 
Rather,  it  is  kept  constant  after  the  first  few  design  iterations. 

These  simplifications  are  reasonable  from  a practical  standpoint  and 
considerable  efficiencies  result  from  them,  especially  for  large  scale 
problems. 

6.7.2  Example  6.2:  Two-Bay  Six-Story  Frame 

Figure  6.4  shows  the  geometry  and  dimensions  of  a two-bay  six- 
story  structure.  In  the  present  formulation,  the  frame  is  divided 
into  three  substructures  by  partitioning  it  at  joints  7 to  9 and 
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13  to  15.  Joints  1 to  3 are  also  treated  as  boundary  nodes,  in  order 
to  provide  a minimum  upper  bandwidth  for  . Members  1 to  10  are  in 
substructure  1 and  members  11  to  20  in  substructure  2.  The  remaining 
10  members  (21  to  30)  are  in  substructure  3.  The  first  two  substruc- 
tures have  18  boundary  degrees  of  freedom  and  the  third  has  only  9. 
Further,  all  substructures  have  9 interior  degrees  of  freedom. 

Table  6.4  gives  design  data  for  this  example.  Design  variable 
linking  is  used  in  order  to  provide  symmetry  in  the  frame  (see  Table 
6.5).  This  reduces  the  number  of  loading  conditions  from  4 (used  in 
Refs.  28  and  41)  to  2,  as  well  as  design  variables  from  60  to  36.  As 
in  the  previous  example,  this  framed  structure  is  optimized  under  two 
separate  constraint  conditions:  stress  constraint  only  (Case  I);  and 
stress  and  displacement  constraints  (Case  II). 

Table  6.5  gives  the  final  designs,  values  of  cost  functions, 
number  of  active  constraints,  maximum  constraint  violations,  ||fb^||2 
at  the  optimum,  the  maximum  value  of  IJfib^ljj  and  the  total  CPU  time. 
The  set  of  active  constraints  at  the  optimum  design  includes: 

Case  I:  Buckling  constraint  for  elements  1,  3,  10,  11,  13,  and  20 

under  loading  condition  1 and  for  elements  9,  16,  19,  21, 

24  to  26,  and  29  under  loading  condition  2,  Von  Mises 
equivalent  stress  constraint  for  element  7 under  loading 
condition  1,  and  minimum  element  thickness  for  all  elements. 
Case  II:  Buckling  constraint  for  elements  1,  3,  4 under  loading 
condition  1,  Von  Mises  equivalent  stress  constraint  for 
element  7 under  loading  condition  1,  displacement  constraint 
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TABLE  6.4  DESIGN  DATA  FOR  TWO-BAY 
SIX-STORY  FRAME 


Modulus  of  elasticity 
Material  density 
Yield  stress 
Factor  of  safety 


= 3 * 10^  ksi 
= 0.2836  lb/in.3 
■ 36.0  k.si 
- 1.51 


Elements 

Initial  values  of  design  variables 

for  (in.) 

Case  I 

Case  II 

R 

t 

R 

t 

1-10 

9.0 

0.20 

9.0 

0.20  | 

11-20 

10.0 

0.22 

11.0 

0.24 

21-30 

12.0 

0.26 

14.0 



0.31 

Lower  limit  on  R = 1.0  in. 
Upper  limit  on  R = None 


Displacement  limits 
in  and  directions = 


Starting  reduction  ratio. 


1.728 

in. 

for 

nodes 

1-3 

1.440 

in. 

for 

nodes 

4-6 

1.152 

in. 

for 

nodes 

7-9 

0.864 

in. 

for 

nodes 

10-12 

0.576 

in. 

for 

nodes 

13-15 

0.288 

in. 

for 

nodes 

16-18 

10. 670 

in. 

for 

all  beam  elements 

-r 

081 

Case  I 

1.0. 

054 

Case  II 

Starting  weighting  multiplier,  w^  = 0.05  for  both  the  cases 
Number  of  loading  conditions  ■ 2 


Jo 


»?r-j 


r - 


TABLE  6.4.  (cont'd) 


Load  data:  " “ " ~~ — 

Loading  Condition  1.  Uniformly  distributed  load  (in  negative 

^“direction)  of  4.0  kip/ft.  on  elements,  1,  7, 
11,  17,  21  and  27;  and  1.0  kip/ft.  on  elements 
2,  6,  12,  16,  22  and  26. 

Loading  Condition  2.  Uniformly  distributed  load  (in  negative 

Xj-direction)  of  1.0  kip/ft.  on  elements  1,  2, 
7,  11,  12,  16,  17,  21,  22,  26,  and  27;  and 
loads  of  9.0  kips  each  at  nodes  1,  4,  7,  10,  13 
and  16  in  x^-direction. 


Group 

Number 

Element 

Numbers 

Optimum  Design  in  inches 

Case 

I 

Case 

II 

Mean  Rad. 

Thickness 

Mean  Rad. 

Thickness 

4* 

1.2 

10.235 

0.2233 

10.009 

0.2183 

3.5 

10.571 

0.2306 

10.361 

0.2261  ' 

4 

5.547 

0.1210 

6.322 

0.1379 

' 

6.7 

9.808 

0.2140 

9.682 

0.2112 

8,10 

7.208 

0.1573 

7.825 

0.1707 

; 

9 

9.140 

0.1994 

10.953 

0.2390 

11,12 

10.082 

0.2200 

10.639 

0.2321 

8 

13,15 

10.535 

0.2298 

9.891 

0.2158 

9 

14 

9.233 

0.2014 

10.763 

0.2348 

10 

16,17 

10.395 

0.2268 

11.406 

0.2488 

11 

18,20 

8.973 

0.1958 

8.687 

0.1895 

12 

19 

11.336 

0.2473 

13.104 

0.2859 

13 

21,22 

11.076 

0.2417 

12.114 

0.2643 

14 

23,25 

10.619 

0.2317 

10.690 

0.2332 

15 

24 

11.236 

0.2452 

11.919 

0.2601 

16 

26,27 

10.986 

0.2397 

11.707 

0.2554 

17 

28,30 

10.687 

0.2332 

11.413 

0.2490 

18 

29 

15.021 

0.3277 

13.915 

0.3036 

Weight  in  lbs 

22530.52 

24405.38 

g 

No . 

of  Active 

3 

e 

•H 

Constraints 

5 5 

15 

cu 

Max 

Constraint 

..-2 

. --2 

° 

u 

c 

Violation 

U. jbxlU 

0.52x10 

1 Ub1 1 |2 

2764.0 

4223.0 

Max.  | | 6b1 1 | 2 

34120.0 

43020.0 

Total  CPU 
Time  In  sec 


87.766 


87.077 
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at  node  7 in  the  x^-direction  under  loading  condition  2, 
and  minimum  element  thickness  for  all  elements. 

The  total  CPU  time  reported  in  Table  6.5  includes  initial  set-up  time 
(2.100  sec  for  both  cases),  and  22  and  26  iterations  of  design  time 
(85.666  sec  and  84.977  sec)  for  Case  I and  II,  respectively.  The  total 
computing  time  reported  herein  also  includes  the  computation  time  re- 
quired for  calculation  of  the  effective  length  factors  (k's)  at  each 
design  iteration.  This  calculation  takes  10%  of  total  computing  time. 

Table  6.6  gives  the  co9t  function  history  for  this  example.  For 
Case  I,  the  final  design  is  achieved  in  22  design  iterations  and  for 
Case  II  in  26  iterations.  As  pointed  out  in  chapters  4 and  5,  the  rate 
of  convergence  is  highly  dependent  on  step  size  and  the  weighting 
matrix.  Lastly,  two  additional  runs  were  made  when  critical  points 
from  the  necessary  conditions  (solutions  of  Eqs.  (6.3-25)  and  (6.3-24)) 
were  also  calculated.  The  Von  Mises  equivalent  stress  constraint  was 
not  active  at  any  of  the  above  critical  points.  However,  the  computa- 
tional time  per  design  iteration  was  increased  by  30%  when  the 
critical  points  were  evaluated. 

6.7.3  Example  6.3:  Helicopter  Tail 
Boom  Structure 

The  helicopter  tail  boom  structure  shown  in  Fig.  6.5  is  idealized 
using  48  frame  elements  (see  Fig.  6.6).  The  total  number  of  degrees 
of  freedom  is  144.  The  element  numbering  system  for  a typical  panel 
is  given  in  Fig.  6.7.  For  convenience,  the  joint  coordinates  and 
element  locations  are  given  in  Tables  6.7  and  6.8,  respectively. 
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TABLE  6.6.  COST  FUNCTION  HISTORY  FOR  TWO-BAY 
SIX-STORY  FRAME 


DI+ 

Number 

Weight  in  lbs. 

Case  I 

Case  1 1 

1 

23336.48 

28536.81 

2 

22723.48 

26286.30 

3 

21441.48 

25000.70 

4 

21819.70* 

24901.38 

5 

22216.96 

24423.88 

6 

22153.72 

23176.56 

7 

22253.80* 

24313.63 

8 

22321.84 

22661.77* 

9 

22416.93 

24325.53 

10 

22414.44* 

24420.39 

11 

22474.35 

23644.30 

12 

22499.20 

24397.14 

13 

22517.74 

23833.90 

14 

22529.27 

24532.80 

15 

22540.89 

23797.65* 

16 

22541.22 

24429.36* 

17 

22546.73 

24547.51 

18 

22545.92* 

24214.25 

19 

22540.41 

24458.22 

20 

22540.69 

24280.10 

21 

22524.64 

24505.66 

22 

22530.48 

24178.72 

23 

22530. 52* 

24459.04* 

24 

24335.72 

25 

24416.48 

26 

24305.77 

27 

24405. 38* 

^Underlined  iteration  is  not  counted  as  it  is  used 
for  checking  the  solution  only. 

* 

Means  change  in  step  size. 

^Design  iteration. 
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TABLE  6.7.  NODAL  COORDINATES  FOR  HELICOPTER  TAIL  BOOM 

STRUCTURE 


Node  Number 


Coordinate  in 

in. 

x 1 

1 

X2 

X3 

0.0 

-11.950 

12.375 

0.0 

-11.750 

-12.375 

0.0 

11.950 

12.375 

0.0 

11.750 

-12.375 

33.50 

-10.666 

11.105 

33.50 

-10.485 

-11.105 

33.50 

10.666 

11.105 

33.50 

10.485 

-11.105 

66.50 

-9.401 

9.855 

66.50 

-9.239 

-9.855 

66.50 

9.401 

9.855 

66.50 

9.239 

-9.855 

99.50 

-8.136 

8.604 

99.50 

-7.994 

-8.604 

99.50 

8.136 

8.604 

99.50 

7.994 

-8.604 

127.50 

-7.063 

7.543 

127.50 

-6.937 

-7.543 

127.50 

7.063 

7.543 

127.50 

6.937 

-7.543 

151.50 

-6.143 

6.634 

151.50 

-6.030 

-6.634 

151.50 

6.143 

6.634 

151.50 

6.030 

-6.634 

173.50 

-5.30 

5.80 

173.50 

-5.20 

-5.80 

173.50 

5.30 

5.80 

173.50 

5.20 

-5.80 
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In  the  present  formulation,  the  idealized  tail  boom  is  divided 
into  three  substructures  by  partitioning  it  at  nodes  9 to  12  and  17  to 
20.  For  the  same  reasons  as  in  the  previous  example,  nodes  25  to  28 
are  also  treated  as  boundary  nodes.  Each  substructure  consists  of  16 
elements;  elements  1 to  16  in  the  first  substructure,  17  to  32  in  the 
second,  and  the  remaining  elements  33  to  48  in  the  third.  The  first 
two  substructures  have  48  boundary  degrees  of  freedom  and  the  third 
has  24.  Further,  all  substructures  have  24  interior  degrees  of 
freedom.  Design  variable  linking  is  employed  in  this  example  (see 
Table  6.10).  Table  6.9  gives  design  data  for  this  problem. 

Based  on  the  conclusions  of  the  previous  two  examples,  transcen- 
dental Eqs.  (6.3-35)  and  (6.3-37)  are  not  solved  for  relative  maximum 
points  of  the  Von  Mises  equivalent  stress.  Rather,  constraints  defined 
by  Eqs.  (6.3-33)  and  (6.3-36)  are  imposed  at  points  of  maximum  bending 
moment  and  maximum  shear  force. 

The  idealized  tail  boom  is  optimized  under  stress  and  displace- 
ment constraints.  Table  6.10  gives  the  final  design,  the  value  of  the 
cost  function,  the  number  of  active  constraints,  the  maximum  constraint 
violation,  | | 6b^  | | at  the  optimum,  and  the  maximum  value  of  | | 6b^  ! ! 2 * 
It  also  gives  the  total  computing  time,  which  includes  problem  set-up 
time  (3.44  sec)  and  12  design  iterations  time  (108.58  sec).  The  set 
of  active  constraints  at  the  optimum  includes:  displacement  con- 
straints at  nodes  25  and  27  in  the  X2~direction  and  minimum  element 
thickness  for  elements  5 to  8,  13  to  16,  21  to  24,  29  to  32,  and 
33  to  48. 
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TABLE  6.9.  DESIGN  DATA  FOR  HELICOPTER 
TAIL  BOOM  STRUCTURE 


Modulus  of  elasticity  * 

(Aluminum) 

Material  density  =* 

Poisson's  ratio 

Modulus  of  rigidity 

Yield  stress  = 

Factor  of  safety  = 

Initial  values  for  mean  rad.,  R 

Initial  values  for  thickness,  t - 

Lower  limit  on  mean  rad.,  R = 

Upper  limit  on  mean  rad.,  R =* 

Displacement  limits  in 

x^,  x^,  and  x^  directions  * 

Starting  reduction  ratio,  r = 

Starting  weighting  multiplier,  wi  = 

Number  of  loading  conditions  = 

Load  data: 


10.5  x 103  ksi 

0.1  lb/in.3 
0.3 

4038.46  ksi 

42.0  ksi 
1.25 

2.0  in . 

0.051  in. 

0.25  in. 

None 

±0.50  in. 

0.081 

0.1 

1 


Node  Number 

Load  Component  (kips)  in  direction 

X1 

x2 

X3 

13 

0.0 

0.0 

-0.140 

14 

0.0 

0.0 

-0.140 

15 

0.0 

0.0 

-0.140 

16 

0.0 

0.0 

-0.140 

25 

1.4903 

1.6918 

0.0 

26 

1.4903 

-1.3658 

0.0 

27 

-1.4903 

1.6918 

0.0 

. 28 

-1.4903 

-1.3658 

0.0 

TABLE  6.10.  RESULTS  FOR  HELICOPTER  TAIL  BOOM 
STRUCTURE 


Group 

Element 

Optimum  Design  in  inches 

Number 



Numbers 

Mean  Radius 

Thickness 

L 

1-4 

2.6695 

0.0880 

2 

5-8 

1.9152 

0.0487 

} 

9-12 

2.6530 

0.0829 

* 

13-16 

2.1035 

0.0535 

5 

17-20 

2.6784 

0.0753 

5 

21-24 

2.1488 

0.0547 

7 

25-28 

2.6238 

0.0673 

8 

29-32 

2.1569 

0.0549 

) 

33-36 

2.5038 

0.0637 

10 

37-40 

2.1730 

0.0553 

11 

41-44 

2.3748 

0.0604 

12 

45-48 

1.9707 

0.0502 

Weight  in  lbs. 

111.20 

B 

No.  of  Active 

a 

Constraints 

10 

•H 

U 

CL 

o 

Max. 

Constraint 

-4 

Violation 

0.92x10 

< 

116b1 

II, 

44.42 

Max 

116b1 

ll2 

618.30 

Total  CPU  Time  in  sec. 

112.02 

1 


~T' 

O'* 
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The  cost  function  history  (in  lbs)  for  this  problem  Is:  69.11, 

86.23,  105.78,  112.36,  112.40,  100.32,  110.49,  111.46,  111.23,  111.16, 

* 

105.23  , 110.96,  111.20  (*  means  change  in  step  size).  A side  computa- 
tion was  made,  in  which  16  more  design  iterations  were  performed  beyond 
the  result  reported  in  Table  6.10.  The  cost  function  was  reduced  by 
0.6%  to  110.64.  There  was  some  redistribution  of  the  material  and 
I | 6b^ | J ^ was  reduced  from  44.42  to  36.75.  This  slow  convergence 
behavior  was  observed  earlier  in  some  structural  design  problems. 

Lastly,  the  value  of  effective  length  factor  (k)  required  for 
evaluating  a ^ (see  Eq.  (6.3-2))  is  taken  as  2.0.  This  value  of  k 
is  held  constant  for  all  design  iterations.  In  the  last  design  itera- 
tion, the  value  of  k for  each  member  was  computed.  Its  maximum  value 
was  1.7.  Thus,  the  effective  length  factor  of  2 is  slightly  on  the 
conservative  side. 

A direct  comparison  of  results  cannot  be  made  with  those  of  Ref. 
63,  as  the  two  idealizations  are  entirely  different. 


6.8  Discussion  and  Conclusions  Based 
on  Results  of  Chapter  6 


In  this  section  results  of  all  example  problems  of  this  chapter 


are  briefly  discussed.  The  one-bay  two-story  frame  is  a small  scale 
structure  and  is  used  for  developing  the  computer  program  for  optimal 
design  of  framed  structures.  The  second  example  is  of  a medium  scale 
plane  frame,  whereas  the  third  example  is  a space  frame  and  is  fairly 


large. 
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The  optimum  weight  obtained  by  the  present  method  for  the  first 
two  examples  cannot  be  compared  with  available  results  in  the  literature 
[28,41],  but  it  is  noted  that  material  distribution  follows  practically 
the  same  pattern  as  in  the  references  cited  above. 

It  is  interesting  to  note  from  the  results  of  these  two  examp Les 
that  during  the  entire  design  process  for  both  problems,  the  Von  Mises 
equivalent  stress  constraint  was  not  active  at  any  of  the  relative 
maximum  Doints  obtained  from  the  necessary  conditions  of  Eq.  (6.3-25). 
Also,  the  effective  length  factor  (k)  does  not  change  appreciably  after 
the  first  few  design  iterations.  Therefore,  some  simplifications  in 
imposing  the  Von  Mises  equivalent  stress  constraint  and  in  evaluating 
the  effective  length  factors  (k's)  are  reasonable  and  considerable 
efficiencies  result,  especially  for  large  structures. 

The  last  example  is  a frame  idealization  of  a helicopter  tail  boom 
structure  that  has  not  been  treated  earlier.  It  is  noted  here,  contrary 
to  the  idealized  wing  structures  of  Chapter  5,  that  the  displacement 
limit  for  this  example  is  quite  severe.  Thus  apart  from  the  lower 
bound  on  thickness  t,  only  the  displacement  constraint  was  active,  which 
was  also  observed  in  Ref.  63. 

It  is  also  noted  here  that  all  design  variables  associated  with  the 
thickness  of  members  for  the  first  two  examples  are  at  their  lower 


bounds  at  the  optimum.  For  the  third  example,  many  thickness  variables 
also  attain  their  lower  bounds.  This  indicates  that  the  lower  bound 


on  the  thickness,  selected  on  the  basis  of  local  buckling  consideration 


CHAPTER  7 


CONCLUSIONS  AND  EXTENSIONS  OF  THE  METHOD 
TO  FAIL-SAFE  STRUCTURAL  DESIGN 

A number  of  key  features  of  the  method  and  problem  solutions  have 
already  been  noted  in  Chapters  A,  5,  and  6.  Sections  A. 5,  5. A,  and  6.8 
are  devoted  to  discussions  of  results  obtained  in  their  respective 
chapters.  This  final  chapter  summarizes  some  of  these  points  and  cites 
areas  that  need  further  study. 

The  general  conclusions,  based  on  the  present  study  are  that  the 
sensitivity  analysis  and  optimal  design  method  developed  with  sub- 
structuring, in  the  state  space  setting,  is: 

(a)  general ; in  the  sense  that  no  restrictions  are  imposed 
on  the  type(s)  of  finite  element (s)  employed,  on  design 
variable  linking  within  or  across  the  substructure 
boundaries,  and  on  the  type  of  behavioral  and  design 
constraints  included.  This  aspect  has  been  demonstrated 
by  considering  a large  class  of  structural  optimization 
problems,  which  consist  of  both  two  and  three  dimensional 
trusses,  idealized  wings,  and  framed  structures. 

(b)  computationally  efficient;  as  compared  with  a similar 
algorithm  without  substructuring  [l]  and  other  available 
methods  in  the  literature  [ 2— A ] , in  the  sense  that  com- 
puting time  per  design  iteration  is  considerably  reduced. 
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Also,  a comparable  number  of  design  iterations  [2-4]  are 
needed  with  the  present  formulation  in  an  interactive  mode 
of  step  size  and  weighting  matrix  selection. 

The  example  problems  of  Chapter  4 serve  to  measure  the  performance 
of  the  algorithm,  by  employing  the  substructuring  concept.  Results  are 
obtained  for  optimization  problems  involving  truss  idealization  of 
structures  under  static  loadings,  with  constraints  on  stresses  and 
nodal  displacements  under  multiple  loading  conditions,  on  natural  fre- 
quency, and  on  design  variables.  Comparing  the  present  results  with  a 
similar  formulation  without  substructuring  [1],  reductions  of  up  to 
23%  for  the  10  member  cantilever  truss  and  up  to  66%  for  the  200  member 
plane  truss  in  computing  time  per  design  iteration  are  achieved. 

Also  comparing  the  present  results  with  other  available  results 
in  the  literature,  it  is  observed  that  for  the  10  member  cantilever 
truss  reductions  in  computing  time  by  factors  of  up  to  5.6  compared 
with  Schmit  and  Miura  [2]  and  2.4  compared  with  Rizzl  [3,4]  ate 
achieved.  For  the  third  example  of  Chapter  4,  the  63  member  wing- 
carry-through  structure,  a reduction  in  computing  time  per  design 
iteration  by  factors  of  up  to  19  is  achieved.  Based  on  the  comparison 
of  the  present  results  with  Refs.  1 to  4,  it  is  concluded  that  for 
large  scale  structures  still  greater  efficiency  in  computing  time  can 
be  achieved.  Further,  in  all  examples  the  optimum  weights  obtained 
with  the  present  method  are  the  same  as  obtained  in  Refs.  1 to  4, 
whereas  they  are  better  than  those  reported  in  Refs.  5 and  6. 

The  results  obtained  in  Chapter  5 for  optimal  design  of  idealized 
wing  structures  serve  to  confirm  the  general  purpose  finite  element 
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structural  synthesis  capability  of  the  present  algorithm.  Constraints 
on  Von  Mises  equivalent  stress  and  on  nodal  displacements  under  multiple 
loading  conditions  and  on  design  variables  are  imposed.  The  optimum 
weights  for  all  examples  obtained  by  the  present  method  are  the  same 
as  obtained  in  Refs.  2 to  4.  As  compared  to  Refs.  7 and  8,  there  is  a 
reduction  of  about  9%  in  the  final  weight  for  the  18  element  wing  box 
beam.  The  computing  times  with  the  present  method  are  significantly 
better  than  those  given  in  Refs.  2 to  4.  For  example  problems  of  this 
chapter,  reductions  in  computing  time  per  iteration  by  factors  up  to 
6.0  compared  with  Schmit  and  Miura  [2],  and  up  to  9.4  compared  with 
Rizzi  [3,4]  are  achieved. 

The  results  obtained  in  Chapter  6 for  optimal  design  of  plane  and 
space  frames  present  a practical  application  of  the  present  formulation. 
The  constraints  imposed  are  same  as  in  Chapter  5,  with  the  addition  of 
a buckling  constraint  for  elements  under  direct  compressive  force 
and  a deflection  constraint  in  the  interior  of  elements.  No  direct 
comparison  of  final  results  is  possible,  since  in  the  present  formula- 
tion constraints  imposed  are  different  than  those  in  Refs.  28  and  41. 
Also,  the  frame  idealization  of  a helicopter  tail  boom  structure  is  new. 

As  mentioned  in  Sections  4.5,  5.4,  and  6.8  the  number  of  analyses 
needed  in  the  present  algorithm  depends  upon  the  step  size  and  weighting 
matrix.  Selection  of  these  parameters  requires  some  experience.  In 
the  present  example  problems,  these  parameters  were  monitored  and  ad- 
justed to  obtain  a final  solution.  More  work  is  needed  in  the  area  of 
step  size  and  weighting  matrix  selection. 
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Extensions  to  Fall-Safe  Structural  Design 

So  far,  the  objectives  of  the  first  phase  of  the  ongoing  research 
project  have  been  accomplished.  An  optimal  design  technique  with  sub- 
structuring has  been  developed,  with  the  associated  design  sensitivity 
analysis  for  various  types  of  structures.  In  the  second  phase  of  the 
project,  extensions  of  the  method  developed  will  be  applied  for  fail- 
safe design  of  large  structures.  In  doing  so,  one  needs  to  modify  the 
statement  of  the  problem  given  in  Section  2.4.  Here,  one  needs  to 
impose  all  the  constraints  of  Eqs.  (2.4-2)  and  (2.2-6),  for  all  damage 
and  loading  conditions.  Also,  one  needs  to  solve  the  state  equations 
(2.3-2),  (2.3-7)  and  (2.2-3)  for  all  damaged  structures.  The  real  advantage  of 
substructuring  becomes  apparent  now,  because  in  reanalyzing  the  structure 
under  damage  conditions,  only  the  matrices  related  to  the  substructure 
where  damage  occurs  are  changed  and  need  to  be  recomputed.  The  design 
sensitivity  analysis  is  also  accomplished  in  a similar  manner.  This 
feature  will  enhance  computational  efficiency  of  the  fail-safe  optimal 
design  algorithm. 

In  imposing  various  constraints,  the  worst  violated  constraint  concepts, 
explained  in  Section  3.4,  can  be  extended  to  include  constraints  for 
damaged  structures.  The  design  sensitivity  analysis  for  the  damaged 
structure  proceeds  as  for  the  case  of  an  undamaged  structure.  Thus,  in 
the  second  phase  of  the  project,  the  fail-safe  optimal  design  problem 
will  be  formulated  with  substructuring.  A computer  code  will  be  developed, 
based  on  this  formulation  and  several  problems  will  be  solved  to  demonstrate 


its  efficiency. 
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where  6^  is  the  same  as  in  Eq.  (2.3-12),  and 


M^r)  = Q(r)  + Q(r)  m£}  Q(r)  4-  Q(r) 


(Al-5) 


is  the  boundary  mass  matrix  for  the  rth  substructure.  Then  the 
approximate  eigenvalue  problem  is  defined  by  the  Eq.  (2.3-16). 
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APPENDIX  2 

FINITE  ELEMENTS  EMPLOYED 

The  element  stiffness  matrices  for  truss,  frame  and  constant 
strain  triangular  (CST)  elements  are  readily  available  in  the  literature 
[34],  and  are  therefore  not  given  here.  However,  element  stiffness 
matrices  for  symmetric  shear  panels (SSP)  and  symmetric  pure  shear  panels 
(SPSP)  cannot  be  easily  found.  Their  derivation  is  briefly  summarized 
here  [2]. 


A2.1  Notations 

a = length  of  SSP  or  SPSP  element 
b = height  of  SSP  or  SPSP  element 
E * modulus  of  elasticity 

&l»Ol  = direction  cosines  of  the  local  x^  axis  in  the  datum 
coordinate  system 


uCx^i^) ' 
v(xlfx2V 


1 * 2 * 3 

X1’X2,X3 

eil,e22,e12 

°ll’°22,a12 

t 

0 


W - * ■ 


displacement  in  local  coordinate  system 

datum  coordinate  system 
local  coordinate  system 

strain  components  in  local  coordinate  system 
stress  components  in  local  coordinate  system 
thickness  of  SSP  or  SPSP  element 
aspect  ratio  of  SSP  or  SPSP  element  (0  - 


> 1 
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r * vector  of  nodal  displacement  In  local  coordinate 
system 

r * vector  of  nodal  displacements  In  datum  coordinate 
system 

B « strain-displacement  relation  matrix 
C,C  ■ stress-displacement  relation  matrices  in  datum  and 
local  coordinate  system,  respectively 
D * stress-strain  relation  matrix 
k,k  *■  element  stiffness  matrices  in  datum  and  local 
coordinate  systems,  respectively 
A = local  to  datum  coordinate  transformation  matrix 


mi 
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u(0,l/2  b)  - 
u(a,l/2  b)  = 
v(0,l/2  b)  ■ 
v(a, 1/2  b)  = 


and 


\ (A2.2-2) 


uCx^O)  = 0.0 


The  local  to  datum  coordinate  transformation  for  nodal  displacements 
is  expressed  as 

r = A r (A2.2-3) 


where 


and 


{il  i2  ?3  f4}T 


r 


{rl  r2  r3  r4 


(A2.2-4) 


(A2.2-5) 


(A2.2-6) 


~r — 
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From  the  assumed  stress  state,  the  strain-displacement  and  the 
boundary  conditions,  the  displacement  state  can  be  obtained  as: 


u(xltx2) 


v(xltx2) 


A11 

0 

A13 

0 

T—i 

CM 

< 

A22 

A23 

A24 

r i. 


(A2.2-7) 


where 


and 


2 X1  X2 


13 


21 


ab 


~ 2 . ~2 

X1  t X1  + VX2  _ vb 
b ab  4a 


X1 

A,,  - 1 - — 
22  a 


A23  “ ” A21 


A24  a 


) (A2.2-8) 


The  strain  displacement  relations  are  obtained  using  Eq.  (A2-5)  as: 

(A2.2-9: 


c = B r 


where 
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e - (en  e22  e12) 


(A2.2-10) 


The  stress-strain  relation  for  plane  stress  is: 


where 


a = D e 


G = {all  °22  °12 ' 


(A2.2-1I ) 


(A2.2-12) 


(A2.2-13) 


1 V 0 


v 1 0 


(A2.2-14) 


0 0 


Substituting  in  Eq.  (A2-12)  the  values  of  strains  in  terms  of 
displacements,  one  obtains  the  stress-displacement  relation  as: 


o » C r 


(A2.2-15) 


r 
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where 


^2 

ab 

0 

2h 

ab 

0 

0 

0 

0 

0 

(A2.2-16) 

1 

1 

1 

1 

2 (l+v)b 

2(l+v)a 

2(l+v)b 

2(l+v)a 

In  datum  coordinate  system,  the  stress-displacement  relation  Is: 


a = C r 


(A2.2-17) 


where 


C = C A 


(A2.2-18) 


The  element  stiffness  matrix  In  the  local  coordinate  system  Is 


B D B dv  = t 


B D B ds 


(A2.2-19) 


Thus,  one  obtains 


k - 


Et 


12(l+v) 


+ 30  -3  2(1^v)  +30  3 


m T1*  30 


3 

0 


3 

0 


3 

0 


1^+30 


3 

0 


(A2.2-20) 
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Thus,  the  element  stiffness  matrix  in  the  datum  coordinate  system  can 
be  expressed  as: 


k = AT  k A 


(A2.2-21) 


Finally,  the  Von  Mises  equivalent  stress  a for  this  element  (see 
Eq.  (5.2-3))  is  given  as 


oC  « (a2n  + 3o22) 


1/2 


(A2.2-22) 


For  calculating  the  maximum  value  of  o from  Eq.  (A2.2-22)  the  following 
expressions  for  o^  and  a ^ are  used  (from  Eq.  (A2.2-15)): 


”11  ' I (r3  - rl> 


(A2.2-23) 


and 


'12  2 (1+v) 


{ 7 <f4  - f2>  + I (i3  + *!>}  (A2.2-24) 


A 2. 3 Symmetric  Pure  Shear  Panel  (SPSP) 

The  element  stiffnes  matrix  for  this  pure  shear  element  (Fig. 
A2.1)  is  also  obtained  by  following  the  previous  procedure  and  by 
assuming  the  stress  state  to  be  as  follows:  = 0,  a 22  " 0,  = 

where  3^  is  a constant.  The  element  stiffness  matrix  is  then  given 


as: 


o< 


The  stress  state  is 

o = C r 


(A2.3-1) 


(A2.3-2) 


(A2.3-3) 


[-**-%*?*■? 
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APPENDIX  3 

SUPPLEMENTAL  EQUATIONS  FOR  OPTIMAL  DESICN  OF  FRAMES 


A3.1  Critical  Stresses  for  Eg.  (6.3-2) 

The  values  of  a and  t for  a member  are  obtained  by  considering 
cr  cr  j n 

a single  force  acting  alone  and  are  briefly  summarized  here.  For 
complete  details,  the  reader  is  referred  to  Ref.  56. 


A3. 1.1  Circular  Tube  Subjected  to  Axial 
Compression  Alone 

A nondimens ional  slenderness  ratio  X for  a member  with  mean 
radius  R and  length  f.  is  given  as  [56]: 


0.450  ^ ^ 

R " E 


(A3. 1-1) 


In  terms  of  X,  the  critical  axial  stress  (a  ) of  a member  can 

ca 

be  expressed  as  follows  [56]: 


for  elastic  behavior,  X < Jl 


(A3. 1-2) 


°ca  * (1  - °*25  X > °y 


(A3. 1-2) 


for  elastic  behavior,  X > /2 


o - a /X 
ca  y 


(A3. 1-3) 


A factor  of  safety  that  is  a function  of  the  slenderness  ratio  is 
often  used  and  is  given  by  the  following  equation  [56]: 


e;  ' \ '-V  V <■ 


'.t/  V ' 

* / 


Jr ' ■ * Ifa.rfW vi  s T 


V 
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For  A < /2, 


FS  = 1.67  + 0.265A  - 0.044A' 


(A3. 1-4) 


and  for  A > , 


FS  - 1.92 


(A3. 1-5) 


2 VAV-  • tfVi  ■ • 

ft 

.... 


1 ' ■ . ■ f ' ■ « • 


T" 


A3. 1.2  Circular  Tubes  Subjected  to  Bending  Alone 
The  buckling  behavior  of  circular  tubes  in  bending  is  similar  to 
that  of  axially  compressed  cylinders  except  that  bent  tubes  have  a 
stress  gradient  [56].  Donnell  [67]  found  the  elastic  buckling  stress  in 
bending  to  be  somewhat  higher  than  the  critical  stress  for  axial  com- 
pression. Fltigge  [68],  and  Timoshenko  and  Gere  [60]  have  reached  the 
same  conclusion.  Other  Investigators  [69-71]  have  indicated  that 
there  is  not  much  difference  between  the  critical  stress  in  bending 
and  in  axial  compression  [56].  In  the  present  work,  the  critical 
stress  in  bending  is  taken  to  be  the  sane  as  for  axial  compression. 

This  is  also  a current  design  practice. 

A3. 1.3  Circular  Tubes  Subjected  to  Transverse 
Shear  Alone 

There  is  very  little  information  available  regarding  tubes 
subjected  to  transverse  shear.  Schilling  [48]  suggests  that  for 
manufactured  tubes  the  critical  stress  in  transverse  shear  be  taken 
as  1.25  times  the  critical  stress  in  torsion  when  buckling  is  elastic, 
and  for  inelastic  buckling  it  should  be  taken  the  same  as  for  torsion. 


7T--  - 
■ ■ 


■ • m 


i 
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A3. 1.4  Circular  Tubes  Subjected  to 
Torsion  Alone 

Using  the  Alcoa  Structural  Handbook  [72],  the  shear  buckling 
stress  (xct)  In  elastic  range,  for  a tubular  member  subjected  to 
torsion  is  given  by  the  following  equation: 


(A3. 1-6) 


in  which  A,  the  equivalent  slenderness  ratio,  is  approximated  by 


A 


3.73  (D/t) 


0.75  w 0.5 


where 


* 

W 


Min 


*,\0.5 


0.5611 


GD 


0.25 


(A3. 1-7) 


(A3. 1-8) 


and  is  the  clear  length  between  circumferential  stiffeners.  In 

the  present  work,  no  circumferential  stiffeners  are  used,  therefore, 

= l.  Further  there  is  not  much  information  available  regarding  the 

inelastic  torsional  buckling.  Hence  in  this  work,  t defined  by  Eq. 

ct 

(A3. 1-6)  is  used  for  both  elastic  as  well  as  inelastic  behavior  which 
is  slightly  conservative. 

Now  the  expressions  for  a and  t , after  simplification,  can 

cr  cr 

be  summarized  as  follows: 


i 
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a 

cr 


4-538’il(£)2  • for*>/2 

{l  - 0.0506  (^)2(^)}^  , for  X < ^2 


(A3. 1-9) 


2 

T 

cr 


0.0983 


2 3 
E t: 

3*2 
R W 


(A3. 1-10) 


for  transverse  shear  acting  along,  and 


2 3 

2 F r J 

t = 0.0629  — rvr 
cr  3*2 

R W 


(A3. 1-11) 


for  combined  torsion  and  transverse  shear,  where  FS  is  defined  by  Eqs. 

* 

(A3. 1-4)  and  (A3. 1-5);  A and  W are  defined  by  Eqs.  (A3. 1-1)  and 
(A3. 1-8),  respectively. 


A3. 2 Computation  of  Q,  Eq.  (6.3-20) 


Figure 


A3.1 


Cross-Section  of  Tube 


V 


From  the  definition  of  Q, 


dQ  = CdA 


where 


C = R sin  0,  and 


dA  “ Rtd6 


Hence  the  total  value  of  Q is 


rir-0 


Q (0 ) “ dQ  =•  2R  t cos  0 

'0 


In  terms  of  x2» 


Q(x2)  = 2R2t  U - (x2/R)2}1/2 


A3. 3 Coefficients  of  Polynomials 


A3. 3.1  Coefficients  in  Eq.  (6.3-25) 


a8  “ 16  r2w5 


a7  “ 2 r2?2w4 


»6  “ \ R2w3(3P2  - P3w  - 3R2w2) 


.2  2 


a - R‘w‘P2(2P‘  - 3P  w - 9R2w2) 


2v-*  2 ji3¥ 
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(A3. 2-1) 


(A3. 2-2) 


(6.3-20R) 


I 

I 

l 


aA  - R w |9R  w - 6Rw2  (^3  P2  - P3w)  + P4  - 6P3P2w 


*1  *’♦($£)} 

a3  = R2w  |36R4P2w3  - 6R2P2w(3P2  - 4P3w)  + 6P2P2w 

-4P3P2+(¥)  V} 

a2  = r2w  |54R4P2w2  - 6R2(P4  - 5P3P2w  + P2w2) 

+ 6P3P2  - 2P3“  + (t5t)  <1  P2>  } 

- R2  |36R4P2w2  - 12R2UP3P2 (PjU  - p2)  - 4P2P2v 

+(!il)2pu 

' AR  ' 2 / 

a0  ■ {9r4p2U  - 6R2p2P3”  + p3"  - {-£;)  (p2  + P3")P3}  r2 


A3.3.2  Coefficients  in  Eq.  (6.3-35) 


a3  = R(F5F6  + F2F3R  } 


a2  - - (I/A)  F1F5  - 1.5  FjF  R* 


a - (I/A)  F F + 1.5  F F.R^ 
1 16  3 4 


a0  * f {r2(F3  • p2)  + F6  “ F5 


mm  .-..jjasrsBa 
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